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^ ; Abstract 

, We present approximate non-perturbative solutions for the propagators as well as the 
^ . running coupling of QCD. We solve a coupled system of renormalised, truncated Dyson- 

■ Schwinger equations for the ghost, gluon and quark propagators in flat Euclidean space- 
^ ' time. We employ ansdtze for the dressed vertices such that the running coupling and the 
^ ' quark mass function are independent of the renormalisation point. At large momenta we 

' obtain the correct one-loop anomalous dimensions for all propagators. Our solutions are 
^ . in good agreement with the results of recent lattice calculations. 

^ ■ In the Yang-Mills sector of Landau gauge QCD we find a weakly vanishing gluon 

g propagator in the infrared and a ghost propagator more singular than a simple pole. This 
is in accordance with Zwanziger's horizon condition and the Kugo-Ojima confinement 
criterion. The running coupling possesses an infrared fixed point at a(0) = 8.92/Nc. To 
P-i^ investigate the influence of boundary conditions on the propagators we solved the ghost 
^ and gluon DSEs also on a four-torus. Our results show typical finite volume effects but 

■ are still close to the continuum solutions for sufficiently large volumes. 

^ ■ In general ghost- antighost symmetric gauges we study the infrared behaviour of the 

c3 I ghost and gluon propagators. No power-like solutions exist when replacing all dressed 
verticed with bare ones. The results of the Landau gauge limit are recovered from a 
different direction in gauge parameter space. 

For the quark propagator we find dynamically generated quark masses that agree 
well with phenomenological values and corresponding results from lattice calculations. 
The effects of unquenching the system are found to be small. In particular the infrared 
behaviour of the ghost and gluon dressing functions found in pure Yang-Mills theory is 
almost unchanged as long as the number of light flavors is smaller than four. 
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Chapter 1 
Introduction 



Quantum Chromo Dynamics (QCD) is the quantum field theory which describes the strong 
interaction of the fundamental building blocks of matter, the quarks and gluons PQ El El 
HlElin]. In contrast to Abelian gauge theories like Quantum Electro Dynamics (QED), 
the non-Abelian nature of the gauge symmetry of QCD not only induces interactions 
between quarks and gluons but also among gluons themselves. This last effect is expected 
to generate the phenomenon of confinement, i.e. the permanent inclusion of all colour 
charges in colour neutral objects, the hadrons. 

'Quantum chromo dynamics is a Lagrangian field theory in search of a solution. ' This 
statement, quoted from the classical review of Marciano and Pagels on QCD [Zj, has not 
lost its relevance since it has been written down in 1977. Although in the meantime a lot 
of progress has been made it is still not clear how the plethora of observed bound state 
objects, the hadrons, can arise from the fundamental quark and gluon fields of QCD. In 
the last thirty years a lot of different strategies have been employed to explore both the 
large and small momentum properties of hadrons. The physical phenomena encountered 
at large momentum transfers are very well described by perturbation theory. Asymptotic 
freedom means that the interaction strength of QCD tends to zero at small distances. High 
energy probes therefore picture hadrons as quark and gluon lumps with definite quantum 
numbers described by so called structure functions. This picture, however, starts to break 
down at energies around 1-2 GeV and is surely inadequate at length scales corresponding to 
the size of the nucleon. At such scales the strong interaction is strong enough to invalidate 
perturbation theory and one has to employ completely different methods to deal with what 
is called Strong QCD. 

There are two phenomena of QCD which are important for this work: the mechanism of 
confinement and that of dynamical chiral symmetry breaking, i. e. the generation of quark 
masses via interactions. Neither of these phenomena can be accounted for in perturbation 
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theory, thus they are genuine effects of Strong QCD. Interestingly, both phenomena appear 
to be connected. From finite temperature studies of QCD we infer that both effects seem 
to disappear at roughly the same temperature but the reasons for this are yet unclear jH]. 

The framework chosen in this work to investigate the small momentum regime of QCD 
are the Dyson- Schwinger equations of motion for correlation functions of the fields. Cer- 
tainly a great step forward in understanding QCD would be the detailed knowledge of the 
basic correlation functions, the propagators. Information on certain confinement mecha- 
nisms is encoded in these two-point functions. Furthermore the mechanism of dynamical 
chiral symmetry breaking can be studied directly in the Dyson- Schwinger equation for 
the quark propagator, which is the gap equation of QCD. Besides being related to the 
fundamental issues of QCD, the quark and gluon propagators are vital ingredients for 
phenomenological models describing low and medium energy hadron physics. Bound state 
calculations based on the Bethe-Salpeter equations for mesons or on the Faddeev equations 
for baryons (for reviews see [01 EI] ) might one day be capable to bridge the gap between 
the fundamental theory, QCD, and phenomenology. 

Throughout this work we will compare the results obtained from Dyson-Schwinger 
equations with those of lattice Monte Carlo simulations (see e.g. IH]). Combining the 
strengths and weaknesses of both approaches allows one to make quite definite statements 
for the propagators in a large momentum range. Lattice Monte Carlo simulations include 
all non-perturbative physics of Yang-Mills theories and are therefore the only ab initio 
calculation method available so far. However, the simulations suffer from limitations at 
small momenta due to finite volume effects. One has to rely on extrapolation methods 
to obtain the infinite volume limit. Furthermore calculations including quarks are subtle 
on the lattice, as it is very difficult to implement fermions with small bare masses. On 
the other hand Dyson-Schwinger equations can be solved analytically in the infrared and 
are the proper tool to assess the effects of dynamical quarks. However, in order to obtain 
a closed system of equations one has to employ ansdtze for higher correlations functions. 
The quality of these truncations can be ascertained by comparison with lattice results. 

We will recall in the next chapter some basic aspects of Strong QCD. Based on the 
symmetries of the (generalised) QCD Lagrangian certain aspects of dynamical chiral sym- 
metry breaking as well as confinement are reviewed. In particular we will recall how 
information on the so called Kugo-Ojima confinement mechanism, the notion of positiv- 
ity and Zwanziger's horizon condition are encoded in the propagators of QCD. From the 
generalised QCD Lagrangian of Baulieu and Thierry-Mieg |E] we will derive a modified 
Dyson-Schwinger equation for the ghost propagator, which together with the correspond- 
ing equations for the gluon and quark propagators are the basic tools of our investigation. 
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The third chapter is devoted to Landau gauge, where a novel truncation scheme is 
introduced that allows to solve the coupled ghost and gluon Dyson-Schwinger equations 
of pure Yang-Mills theory. Contrary to earlier attempts this is done without any angular 
approximations in the loop integrals of the equations. Besides the numerical solutions for 
general momenta we obtain analytical results for the ultraviolet and infrared region of 
momentum. We arc thus able to show that the anomalous dimensions of one-loop pertur- 
bation theory arc reproduced by our solutions for the full ghost and gluon propagators. 
For small momenta the ghost and gluon dressing functions follow power laws, which are 
in accordance with the Kugo-Ojima confinement criterion as well as Zwanziger's horizon 
condition. We are able to show that only one out of two infrared solutions already found in 
earlier investigations is connected to the numerical solutions for general momenta. The re- 
sulting running coupling possesses an infrared fixed point. Our results for the propagators 
are in nice agreement to recently obtained lattice calculations. 

Landau gauge is a special gauge in the sense that it allows for surprisingly simple vertex 
ansdtze in the truncation of the Dyson-Schwingcr equations. In chapter four we will explore 
whether a truncation employing bare vertices can be extended to general gauges. We solve 
the corresponding Dyson-Schwingcr equations analytically in the infrared employing power 
laws for the dressing functions. The main results of Landau gauge, an infrared vanishing 
gluon propagator and an infrared diverging ghost turn out to be persistent for the class 
of linear covariant gauges. However, no power law solutions for general ghost-antighost 
symmetric gauges can be found. Numerical solutions for the limit of Landau gauge from 
different directions in gauge parameter space turn out to be stable. 

In chapter five wc change the base manifold and investigate the Dyson-Schwingcr equa- 
tions of Landau gauge Yang-Mills theory on a four-torus. There arc three ideas motivating 
such a change. The first idea stems from the observation that torus calculations share the 
finite volume problem with lattice Monte Carlo simulations. However, as the continuum 
limit is known in the framework of Dyson-Schwinger equations, one is able to judge extrap- 
olation methods on a torus. Our solution for the gluon propagator on a torus resembles 
closely the one found in the continuum whereas the ghost dressing function deviates in 
the very infrared. The second idea is a technical one. Although the treatment is inverted 
in this work, we first have been able to obtain results without angular approximations 
on the torus and only subsequently in the continuum. This is due to the torus acting as 
an effective regulator in the infrared. The third idea, which remains for future work, is 
to include topological obstructions like twisted boundary conditions in Dyson-Schwinger 
equations on a torus. 

The last chapter of this work focuses on the quark propagator. We present solutions 
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for the quenched system of quark Dyson-Schwinger equations using the results of chapter 
three as input, and compare to lattice calculations. We then go one step further and solve 
the unquenched coupled system of equations for the ghost, gluon and quark propagators 
of QCD. This is done without any angular approximations. Our results again reproduce 
the anomalous dimensions of one- loop perturbation theory in the ultraviolet. The effects 
of unquenching the system, i.e. including the quark loop in the gluon equation, is found 
to be small if the number of light flavours is small, A^j*^'^* < 4. In particular the infrared 
behaviour of the ghost and gluon dressing functions remain unchanged. The quark masses 
generated by dynamical chiral symmetry breaking are found to be close to phenomenolog- 
ical values. 



Chapter 2 

Aspects of Strong QCD 



2.1 The generating functional of QCD 

Working in Euclidean space-time^ the generating functional of the quantum field theory 
of quarks and gluons is given by 



where we have introduced the Grassmann valued sources fj and t] for the quark fields \i/ 
and ^ and the source for the gauge field ^4°. Furthermore we used the abbreviation 
p= "^i^Dfi with Euclidean 7-matrices^ and the covariant derivative given in eq. ()2.4j) . 
The quark fields are spin- 1/2 fermions which transform according to a fundamental rep- 
resentation of the gauge group SU{Nc). The central focus of this work is QCD, i.e. the 
gauge group SU{3). However in the course of our investigations we will run across some 
results, that are valid for general gauge group SU{Nc). Some comparisons with lattice 
calculations will be done for SU{2). 

The non-Abelian gluon fields, A^, transform according to the adjoint representation of 
the gauge group. The corresponding field strength tensor and the covariant derivative in 

^We will adopt Euclidean metric throughout this work. A justification of this choice will be given in 
subsection 12.3.21 

^We use hermitian 7-matrices defined in appendix lA. II 




(2.1) 
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2.1. The generating functional of QCD 



the adjoint representation are given by 

F;, = d.At-d.A^-gr'^AlAt, (2.2) 
Df = d^S'^' + gr'^A^. (2.3) 

Here g is the (unrenormahsed) couphng constant of the theory and /"^^ are the structure 
constants of the gauge group. With the help of the generators t"' of SU{Nc) we can rewrite 
the covariant derivative in the fundamental representation 

Df, = d^ + igA^ , (2.4) 

with A^ = A^^r and = if'H^. The Lagrangian C = ^ {-p+ m) ^ + ^F^^ of our 

theory is invariant under local gauge transformations. 

One of the most intricate tasks in the quantisation of a field theory is the separation 
into physical and non-physical degrees of freedom, which is a prerequisite for the definition 
of the physical state space of the theory. The integration over all possible gauge field con- 
figurations A in the generating functional ()2.1|) includes the ones that are gauge equivalent. 
Therefore the integration generates an infinite constant, the volume of the gauge group 
Q, which has to be absorbed in the normalisation. More important, the gauge freedom 
implies that the quadratic part of the gauge field Lagrangian has zero eigenvalues and 
therefore cannot be inverted^. This prevents the definition of a perturbative gauge field 
propagator [T8] . 

In order to single out one representative configuration from each gauge orbit 

[A^] ■= {A^ = UAU^ + UdU^ : U{x) G SU{N^)] (2.5) 

one has to impose the gauge fixing condition f°'{A) = on the generating functional. This 
is conveniently done by inserting the identity 

1 = A[A] j Vg6ir{A)) (2.6) 

into the generating functional 1)2.11) and absorbing the group integration J Dg in a suitable 
normalisation jT3]. We will discuss problems with this gauge fixing prescription in more 
detail in subsection 12.3.31 

•^If one avoids the generating functional and employs canonical quantisation this problem manifests 
itself on the level of commutation relations of the fields. These are usually fixed at time zero and should 
then determine the commutators for all times. However, this cannot be the whole story, since one can 
always gauge transform to a field that vanishes at time zero. So one has to remove the freedom of gauge 
transformations here as well 
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In linear covariant gauges the Faddeev- Popov determinant A[A] reads explicitly 

A[A] = Det {-d,D;') (2.7) 

and can be written as a functional integral over two new Grassmann valued fields c and c. 
Furthermore the gauge fixing condition f"'{A) = d^A^ — iXB"" = employing the auxiliary 
field B can be represented by a Gaussian integral centred around iXB°- = (see e.g. jHI 
for details). We then have 

A[A]6{r{A)) = j V[c, c, B] exp |- y d'x (^-id^cD^c + iB'^d^A'^^ + ^B'^B'^^ | 

= J V[c, c] exp |- y d'x (^-td,cD,c - } ■ (2.8) 

Introducing the sources a and a for the antighost and ghost field respectively we arrive 
at the gauge fixed generating functional 

Z[J,a,a,r],f]] = j ©[A^f^cc] 

exp | - y d'^xCeff + j d^x [A^J^ + f]^ + ^r] + ac + ca)^ , (2.9) 

with the effective (unrenormalised) Lagrangian 



\2 



C,ff = Cm + Ca + Cgf = ^ {-P+ m)^ + -^F^ + - id^cD.c. (2.10) 



Note that a factor of i appears in front of the ghost terms as we have used real ghost and 
antighost fields. We will see the importance of this choice in section 12. HI 

In effect we have modified the Lagrangian of our theory with a term containing the 
(unphysical) ghost and antighost fields c and c which transform according to the adjoint 
representation of the gauge group and the gauge fixing part which may or may not be 
written with the help of the Nakanishi-Lautrup auxiliary field -B". 

Apart from the convention of real ghost fields the Lagrangian ()2.1()|1 is the usual one 
employed in perturbation theory. It has some important properties: 

(i) it is of dimension 4, 

(ii) it is Lorenz invariant and globally gauge invariant, 

(iii) it is BRS and anti-BRS invariant^, 



^The explicit definitions of the BRS and anti-BRS transformations used in this work are given in 
subsection 12.2.31 
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2.1. The generating functional of QCD 



(iv) and hermitian. 

The number of dimensions and Lorenz invariance are certainly dictated by experiment. 
The further symmetries, global gauge invariance and BRS symmetry, are discussed in 
section 1^7^ Hermiticity is necessary to define the physical S-matrix of the theory. We will 
say a little more about this issue in subsection 12.3.11 

The Lagrangian ()2.10|) arises from a specific gauge fixing procedure, the Faddeev-Popov 
method. This, however, is not the only gauge fixing procedure that has been employed 
so far. Indeed, as will be discussed in more detail in subsection I2.3.3| the Faddeev-Popov 
method is not capable to fix the gauge completely. Although it is currently not known to 
what extend this poses a problem for strong QCD, it is desirable to develop alternatives 
which do not suffer from such a deficiency. Examples for different gauge fixing procedures 
are topological gauge fixing [IHllini, where the aim is to represent the partition function 
of QCD by a topological invariant, or stochastic gauge fixing [T2J [THj, which employs a 
Fokker-Planck equation for the probability distribution in gauge field space^. 

On the other hand, one could reverse arguments and claim the properties (i)-(iv) to be 
crucial for the quantum field theory of strong interaction. Without bothering about the 
details of the gauge fixing procedure one could then search for the most general Lagrangian 
satisfying (i)-(iv). This view has been adopted in [I2, 19j. It has been shown that, omitting 
topological terms, the most general polynomial in the fields A^, c, c, \E' and ^ satisfying 
(i)-(iv) can be written 



Here the abbreviation (c x c)" = gf"'^'^ c^c'^ is used. Again both ghost fields, c and c, are 
chosen to be real, which is necessary here to maintain the hermiticity of the Lagrangian 
for all values of the gauge parameters A and a, see e.g. [20] and references therein. 

We easily see, that the new Lagrangian p.llj) is a generalisation of the Faddeev-Popov 
Lagrangian ()2.1()|1 with a new, second gauge parameter a. This gauge parameter controls 
the symmetry properties of the ghost content of the Lagrangian. For the cases a = 
and a = 2 one recovers the usual Faddeev-Popov Lagrangian ()2.10|) and its mirror image, 
respectively, where the role of ghost and antighost have been interchanged. For the value 
a = 1 the Lagrangian 1)2.111) is completely symmetric in the ghost and antighost fields. 
Compared to the Faddeev-Popov case the four ghost interaction is an additional term in 

^ A pedagogical treatment of this topic can be found in |2] . 



c 




c. 



(2.11) 
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the theory. Note that such a term is e.g. found in topological gauge fixing scenarios ^Hj 
or as a result of partial gauge fixing in maximal Abelian gauges \Z\\ |22l IIB] ■ 

In ref. ^2] it has been shown that the S-matrix of the theory ()2.11|) is invariant under 
variation of the gauge parameters A and a. Therefore gauge invariance of physical observ- 
ables is ensured. One-loop calculations confirm in particular the independence of the first 
nontrivial coefficient of the /3-function from the gauge parameters. 

Furthermore, the existence of a renormalised BRS-algebra has been proven ^2]; thus 
the theory given by ()2.1H) is multiplicatively renormalisable. From one-loop calculations 
one finds that the Faddeev- Popov values of the gauge parameters, a = and a = 2, are 
fixed points under the renormalisation procedure. The same is true for the ghost- ant ighost 
symmetric case a = 1. The case of Landau gauge, A = 0, corresponds to a fixed point as 
well, because the constraint d^A^ = is not affected by a rescaling of the gluon field. 

The correspondence between the bare Lagrangian ()2.11|) and its renormalised version 
including counterterms is given by the following rescaling transformations 

Al y%A1, Z^ec\ ^ Z2^^, (2.12) 

9 ZgQ, a^Z^a, A ZaA, (2.13) 

where six independent renormalisation constants Z3, Z3, Z2, .^g, Z^ and Z\ have been in- 
troduced. Furthermore five additional (vertex-) renormalisation constants are related to 
these via Slavnov-Taylor identities, 

Zi = ZgZ^ , Zx = ZgZ^Zr^^ , Zip = Z gZ ^ ^2, Z = Z^Z^, Z^ = Z^Z^. (2.14) 

Most of the calculations performed in this work are done in Landau gauge. By partial 
integration it is easy to see that in Landau gauge the additional gauge parameter a drops 
out of the Lagrangian 1)2.111) due to the condition S^A^ = 0. Our results are therefore 
completely insensitive to a decision between the usual Faddeev-Popov version of QCD 
or the generalised version. The only exception occurs in chapter |3J where we investigate 
the infrared behaviour of the ghost and gluon propagators for general values of the gauge 
parameters A and a. 

2.2 Symmetries of QCD 

In the following we will recall some of the symmetries of the generalised Lagrangian 1)2.111) . 
While ghost number symmetry and BRS symmetry are vital ingredients for the Kugo- 
Ojima confinement scenario discussed in the next section, chiral symmetry and its dynam- 
ical breaking will be important in chapter IHl when we solve the Dyson-Schwinger equation 



10 



2.2. Symmetries of QCD 



of the quark propagator. All considerations in this section are independent of the specific 
value of the gauge parameters A and a of our general Lagrangian ()2.1H1 . 

2.2.1 Chiral symmetry 

Let us first discuss the quark sector of the Lagrangian ()2.1H1 . The approximate chiral 
symmetry of the quark terms in the Lagrangian proves to be very fruitful to generate low 
energy expansions of QCD (for reviews see e.g. 123 123 I2Z|)- A lot of qualitative and 
quantitative properties of hadrons have been inferred from the concept of approximate 
chiral symmetry. Effective models like the global colour model PB| or the NJL model (jSHl 
OU] . see also [211 ^ review) try to capture the basic properties of QCD by approximating 
the gluon sector with a simple effective interaction but maintaining the chiral symmetry 
aspects of the quarks. 

From eq. (j2.1H) we have the quark part of the QCD Lagrangian 

= t (-|)+m)^ , (2.15) 

where m is a diagonal matrix containing the masses of six different flavours of quarks, 
generated in the electroweak sector of the standard model (see e.g. jHS])- The chiral 
limit, m = 0, is appropriate for the three light quarks up, down and strange, which are 
considered only in the following. In the chiral case the quark fields can be split in left- 
and right-handed Weyl-spinors 

^ii = ' = (2-16) 

The resulting Lagrangian is symmetric under the global unitary transformation SU(3) x 
SU(3) X U(l) X U(l), which generates the currents 

= 'S'Tm^, (2.17) 

jj = ^7m75^, (2.18) 

3i = ^i,t^^, (2.19) 

fj = ^7/.75t^^, (2.20) 

where t-f = ^ denote the generators of SU(3) flavour transformations given by the Gell- 
Mann matrices A-^. These currents are conserved on the classical level of the theory. 
Quantum corrections, however, spoil the conservation law for the axial current, eq. ()2.18p . 
This effect is known as Adler-Bell-Jackiw anomaly and has the observable consequence of 
allowing the otherwise forbidden decay of the uncharged pion 7r° into two photons. 
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Under the presence of a non-vanishing mass matrix m and including the effects of the 



anomaly we have the divergences of these currents given by 

d^j, = 0, (2.21) 

d% = 2z^m^,^--£^e^'"^''F^,F^^, (2.22) 

d^^jl = ^[tf,m]^, (2.23) 

df'jf/ = ^{t^^m}^. (2.24) 



Thus only one current, eq. (jSSJ, is conserved and describes baryon number conservation 
in strong interaction processes. The vector current, eq. ()2.23|) . is conserved in the case of 
identical quark masses and thus describes the approximate flavour symmetry in the light 
quark sector of QCD. 

The axial vector current, eq. ()2.24|) . is broken if we have a non- vanishing quark mass 
matrix in the Lagrangian of QCD. This situation is called explicit chiral symmetry breaking. 
Since the current quark masses of at least the up and down quark are very small, we still 
expect approximately degenerate parity partners of the lowest lying hadron spectra, if the 
current masses were the only reason for broken chiral symmetry. However, such parity 
partners are not observed in nature. 

The solution of this puzzle is the effect of dynamical chiral symmetry breaking described 
by the Dyson-Schwinger equation for the quark propagator. We will see in detail in chapter 
ini how the strong interaction in the quark equation generates physical quark masses of the 
order of several hundred MeV even in the chiral limit of zero bare masses in the Lagrangian 
of our theory. 

2.2.2 Ghost number symmetry 

The conserved Faddeev-Popov ghost number Npp is generated by the scale transformation 

e e-®c" , (2.25) 

with a real parameter G. Via the Noether-theorem this symmetry leads to a conserved 
current and a conserved charge Qc, the FP ghost charge. The charge leaves all other fields 
invariant and acts on the ghost fields as 

[2Q„c'^(x)] = c'^(x), 

[iQc,c\x)] = -c\x). (2.26) 
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The ghost number Npp is then identified with the eigenvalue of the operator Qc multiphed 
by i. Note that the appearance of the hermitian operator Qc with purely imaginary 
eigenvalues is perfectly consistent in the presence of indefinite metric, which seems to be 
the case in quantum field theories [201 • Note further, that only the scale transformation 
(12 .2511 and not the corresponding phase transformation is compatible with the choice of real 
ghost and antighost fields. Ghost number conservation and the BRS-symmetry discussed 
in the next subsection play an important role in the Kugo-Ojima confinement scenario 
summarised in subsection I2.3.11 

2.2.3 Global gauge symmetry and BRS-symmetry 

Although the Lagrangian ()2.11|) is gauge fixed such that local gauge invariance is not 
present any more, there are two gauge symmetries left: The global gauge symmetry and 
the so called BRS-symmetry which has been found by Becchi, Rouet and Stora^ jIMj . 
The global gauge transformations of the gauge field and the quark field are given by 



with space-time independent parameters A'^ and the generators t"" of the gauge group. 
Although the global gauge transformation is a symmetry of the Lagrangian it is not clear 
whether a corresponding well defined charge exists, i.e. whether global gauge symmetry 
is spontaneously broken or not. This will play an important role in the discussion of the 
Kugo-Ojima confinement criterion in subsection I2.H.1I Note that there are no fundamen- 
tal reasons why the global gauge symmetry should be unbroken, as Elitzurs theorem of 
unbroken gauge symmetry only applies to local transformations 

The most efficient way to introduce BRS-symmetry is by means of the nilpotent BRS- 
operator s. The BRS-transformation of the gluon, ghost and quark fields as well as the 
auxiliary field B are given by 




(2.27) 
(2.28) 



s^ 



-igt"-c"-^ , 



(2.29) 




(2.30) 
(2.31) 
(2.32) 
(2.33) 



sB' 



•a 



. 



^At the same time the symmetry has been discovered independently by Tyutin, see 



Chapter 2. Aspects of Strong QCD 



13 



Note that the Nakanishi-Lautrup auxihary field B can be ehminated from the BRS- 
transformations by using its equation of motion, dfj,A^ — iXB"' = 0. Note furthermore 
that the apphcation of the BRS-operator s on a field increases the ghost number by +1, 
thus we can assign the value Npp = +1 to the BRS-operator itself. 

The BRS-transformations ()2.33|) can be seen as (local) gauge transformations with the 
ghost field c{x) as parameter. Thus the transformations describe a global symmetry, since 
one is not free to treat different space-time points independently. Similar to global gauge 
symmetry, it is not clear whether the BRS-symmetry generates a well defined BRS-charge 
Qb- Indeed, it has been argued [33113 that BRS-symmetry is broken as a consequence 
of the presence of Gribov copies^. The Kugo-Ojima confinement scenario, discussed in the 
next section, assumes a well defined, i.e. unbroken, BRS-charge Qb- 

With the help of the BRS-operator s the Lagrangian 1)2.111) can be written as 

C = Cm + Ca + Cgf (2.34) 
= Cm + CA-ts (^d^c'^A; - h-B- + t^gp^'^c'^c^c^^ . (2.35) 

Thus the BRS-invariance of the Lagrangian can be inferred from the local gauge invariance 
of the part Cm + Ca and the nilpotency of the BRS-operator, = 0. 

Under the assumption of its existence the BRS-charge Qb together with the ghost 
charge Qc constitute a simple algebraic structure characterised by the following relations 

{Qb.Qb} = 2(Qb)2 = 0, (2.36) 
[iQc.Qb] = Qb, (2.37) 
[^Qc,Qc] = 0. (2.38) 

This algebra is called BRS-algebra in what follows. 

It is interesting to note that the Lagrangian 1)2.111) is also invariant under the anti-BRS 
transformations defined by 

= -igt^e^i , (2.39) 
-sAl = Dfc", (2.40) 
sc" = —r^^'c^c" , (2.41) 

sc" = -iB'' - p^^^c'^c" , (2.42) 
= -gT^^c^B^ . (2.43) 



^The problem of Gribov copies is discussed in some more detail in subsection 12. 3. 31 
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Similar to the BRS-operator s the anti-BRS operator s is nilpotent and both operators 
are related by 

ss + ss = (2.44) 

It has been argued, however, that although adding structure to the mathematical frame- 
work of the theory the presence of anti-BRS symmetry has no influence on the physical 
content [201 • Consequently this symmetry will play no role in the following discussions. 

2.3 Aspects of confinement 

In this section we will summarise some aspects of confinement. All three topics discussed, 
the Kugo-Ojima confinement scenario, the notion of positivity and Zwanziger's horizon 
condition generate testable predictions for the behaviour of the propagators of QCD. 

2.3.1 The Kugo-Ojima confinement scenario 

It has been stated already above that one of the most intricate problems in quantum field 
theories is the separation of physical and unphysical degrees of freedom. In QCD this 
problem is directly connected with the issue of confinement, since we are searching for the 
mechanism which eliminates the coloured degrees of freedom from the physical state space 
Vphys of the theory, which is supposed to contain the colourless hadronic states observed 
in experiment. 

From a theoretical point of view to be able to define a physical S-matrix between the 
physical states of the theory three conditions should be satisfied fIU\ : 

• The Hamiltonian H corresponding to the Lagrangian of the theory should be her- 
mit ian. 

• The physical subspace Vphys of the state space of the theory should be invariant 
under time evolution, i.e. HVphys ^ Vphys- 

• The physical subspace should be positive semidefinite, i.e. {"^l"^) > 0, if I^E') G Vphys- 

Certainly, the first of these criteria holds, because the general Lagrangian 1)2.111) is her- 
mitian due to our choice of real ghost fields^ c and c. The second criterion suggests 
the definition of the physical subspace via a conserved charge, which commutes with the 
Hamiltonian of the theory and thus guarantees the invariance of the subspace under time 

*Note that for general values of the gauge parameters A and a this is not the case in the original version 
of the Lagrangian in ref. ^2], where complex ghost fields have been chosen. 
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evolution. The third criterion is necessary to allow for the usual probabilistic interpreta- 
tion of the quantum theory. In general the complete state space V has indefinite metric 
and one therefore has to prove explicitly, that the third criterion holds. 

To proceed we recall briefly on an intuitive level what is meant by the notion of asymp- 
totic states^. What is really observed in particle physics are not fields but particles. Such 
particles are present long before and after scattering processes and are described by so 
called asymptotic states. These states |$as) are created by asymptotic fields which 
are defined by the weak operator limit of the corresponding field operators $ for large 
absolute times [HH] : 

(mix) - <^Ux)\g) if x^^ioo, (2.45) 

for any two states \f) ,\g) E V. 

The asymptotic states constitute two asymptotic state spaces, Vm for x° — — cxo and 
Vout for — i> +00, which can be shown to be isomorphic to the Fock space of free fields. 
One of the crucial postulates of axiomatic field theory, called asymptotic completeness, 
states the equivalence between the asymptotic and the complete state spaces of the theory: 

Vin = V = Vout (2.46) 

This is most important to define the S-matrix as unitary transformation between the in- 
and out-states of the theory. Whereas the complete S-matrix acts on the whole state space 
V, the physical S-matrix acts on the space Vphys C V of physical states only. 

Once one has succeeded to define Vphys, it is necessary to show that it only contains 
colourless states according to the confinement hypothesis. Based on symmetries described 
in the last section the Kugo-Ojima confinement scenario describes a mechanism, by which 
such a positive (semi-) definite state space containing only colourless states is generated 

One of the basic assumptions of the scenario is the existence of the conserved BRS- 
charge Qb, which is used to postulate the physical subspace Vphys of the state space V 
byio 

Vphys = {\phys) : Qnlphys) = 0}. (2.47) 

This space can be shown to have a positive semidefinite metric Recalling the BRS- 
transformations given in the last section the space Vphys contains two different sorts of 
states. The first ones are the so called BRS-daughter states, Each of these states can 

^See |37) for a mathematical rigorous introduction into the concept of asymptotic states and the 
problems related with asymptotic bound states as well as asymptotic massless particles. 

^°The corresponding construction d^A^\phys) = in QED is known as Gupta-Bleuler condition. 
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be generated by applying the BRS-operator to a corresponding parent state |7r), which is 
not element of Vphys- We thus have |0) = s|7r). The second ones are BRS-singlet states for 
which no such parent states exist^^. 

Let us first discuss the BRS-daughter states. The annihilation of these states by the 
BRS-charge is a trivial consequence of the nilpotency of the BRS-operator, which 
implies that QbQb = 0. Due to the BRS-algebra ()2.38|) of the ghost charge Qc and the 
BRS-charge Qb daughter states |0) and their parents |7r) always occur in pairs, i.e. two 
daughters and two parents form a so called BRS-quartet. Denoting the ghost number 
by Npp we have the quartet (|7r, Npp), |0, Npp + 1), |vr, —Npp — 1), \4>, —Npp)) related by 
the BRS- and ghost number-transformations: 



Qb\t!-,Nfp) 
Qb\7i,-Nfp-1) 

Q,\71,Nfp) 
Qek, -Npp - 1) 



^, Npp + 1) 
-Nfp) 
-Npp) 

^, Npp + 1) 



(2.48) 
(2.49) 
(2.50) 
(2.51) 



It is easy to show, that the BRS-daughter states |0) are orthogonal to all states |\E') of 

^phys '^jS 

(^|0) = (^|gB|7r) = O. (2.52) 



As a consequence the daughter states do not contribute to the physical S-matrix, i. e. the 
corresponding asymptotic states are not part of the physical spectrum of the theory. This 
is also true for the asymptotic states of parent states ISHj- We therefore have the result, 
that the asymptotic states of all members of BRS-quartets do not correspond to physical 
particles. This confinement of quartet states is known as quartet mechanism. 

As an example the so called elementary quartet can be constructed, which consists of the 
parent states \A"j), \c"') and the daughters \D"^c^), \B"'). From these states corresponding 



^^In geometrical language the space of BRS-singlets, H{Qb, V), is called a cohomology, the physical state 
space Vphys = Z{Qb,V) is denoted as cocycle space and the space of BRS-daughter states, B{Qb,V), 
is called a coboundary space. The cocycle space Vphys contains the closed forms with respect to the 
BRS-charge, whereas the coboundary space contains the exact forms (see e.g. 00]). 
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asymptotic states of massless particles can be inferred [SH] 



c^''-{x) = 7"(x) + ... (2.53) 

^a«,a^^^) = + (2.54) 

(D^cr'^ix) = d,rix) + ... (2.55) 

B''''''{x) = (3%x) + ... (2.56) 



It turns out, that these asymptotic states describe ghosts, antighosts and longitudinally 
polarised gluons, which are therefore confined by the quartet mechanism. 

We now return to the remaining states in Vphys, the BRS-singlets |$). The asymptotic 
states of the BRS-singlets are candidates to describe the physical particles of the theory, 
namely baryons and mesons. It can be argued jHS], that all BRS-singlet states have 
vanishing ghost number Npp, as it is expected for physical states. In order to have 
confinement one has to show that 

• there is a well defined, i.e. unbroken, global colour charge with 

($|g''|$') = (2.57) 

for all BRS-singlet states |$), |$') G Vphys, 

• the cluster decomposition property is violated for these states. 

As the cluster decomposition property is not investigated in this work we briefly explain 
what it means and put it aside afterwards. Speaking intuitively, cluster decomposition 
means the possibility to divide each given lump of particles into subsets which can be torn 
apart. For QCD this is not what is observed in experiment, as the cluster decomposition 
property would imply the possibility to split an observable colourless object into observ- 
able coloured objects. A rigorous mathematical formulation of the cluster decomposition 
property is given in |1T]. Here we just mention that the cluster decomposition property 
may fail only for a quantum field theory with an indefinite metric and without a mass 
gap in the whole state space V |^. This certainly not excludes a mass gap in Vphys as is 
expected in the case of confinement in QCDjl2j. 

Let us now come back to the global colour charge Q"^. This charge can be defined by 

Q'' = j d^xJ^{x), (2.58) 

^^In this is done by a thorough analysis of the correlation functions {D'^^'c^c"') and The 
intuitive argument given in jSj is based on the free field limit g = 0, which is certainly not appropriate in 
a strong coupling gauge theory. 
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where the divergence free current is generated by the unfixed global gauge symmetry 
c.f. subsection 12.2. HI With the help of the equations of motion the current can be 
written as"^^ 

9j;i = d,F;:, + {QB,Dfc^}. (2.59) 

For obvious reasons this equation is called quantum Maxwell equation. One can thus 
separate the global colour charge Q°' into two different charges and A^'* corresponding 
to the two terms of eq. ()2.59|) : 

(2.60) 



Q" = J d^x^ {diFS, + {Qb, Dfc'}) =G^ + N\ 



One crucial point with a Noether current is that it is only defined up to an arbitrary term 
of the form of a total derivative. Thus if eq. ()2.60|) is well defined, one can redefine the 
global colour charge, 



d^x^{QB,Dfc^} = N". (2.61) 

If the second equation is well defined, we immediately satisfy eq. ()2.57j) due to the nilpo- 
tency of the BRS-charge, {QbY = 0. 

The problem is, however, that there is the possibility that the three dimensional integral 
in eq. 1)2.611) does not converge i. e. global gauge symmetry is spontaneously broken. One 
version of the Goldstone theorem states the equivalence of the following conditions on a 
conserved current and its global charge Q jSHj: 

(a) Q = J d^xJo does not suffer from spontaneous symmetry breakdown; 

(b) contains no discrete massless spectrum: (0| J^|^(p^ = 0)) = 0. 

However, we already know from eq. ()2.5Hj) that the antighost c is a member of the 
elementary quartet and contains a one-particle contribution from the massless asymptotic 
field 7. We therefore have 

D^(^ = d^c'' + gr'''A''^(^ {l + u)d^r ifxo^ioo. (2.62) 

Here the proportionality factor 1 stems from the first term of the covariant derivative and 
the proportionality factor u from the term containing the gauge field. As the global colour 
charge is proportional to -D^c", see eq. ()2.6ip . we arrive at the condition 

l + u = (2.63) 
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Our treatment follows closely ref. |4H| . 
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for the global colour charge to be well defined [HH] . 

This condition has been derived for general linear covariant gauges^^, a = 0. A special 
case is Landau gauge, which is our preferred choice for most of the investigations in the 
next chapters of this thesis. In Landau gauge it has been shown, that the condition ()2.63|) 
is connected to the ghost dressing function G{p'^). With the definition 

Aghast = -2 (2-64) 



for the ghost propagator one obtains the relation 

(2.65) 




We are thus in the position to check the Kugo-Ojima confinement criterion, eq. (j2J 
by calculating the ghost dressing function G{p'^) in the infrared. This is most conve- 
niently done in the framework of Dyson-Schwinger equations, as analytic expressions for 
the infrared behaviour of the propagators can be obtained. 

The condition eq. ()2.65p has also been subject to various investigations on the lattice. 
There one aims either at the direct determination^^ of the parameter u EHl , or 
determines the ghost dressing function in the lattice analogue to Landau gauge [HI HHl EH! • 
Whereas exploratory calculations of the first class give a result of u ~ —0.8, the ghost 
dressing functions obtained in the second class of investigations indicate a divergence when 
extrapolated to zero momentum and thus are in agreement with the criterion eq. ()2.f)5j] . 
However, one has to keep in mind that infinite volume extrapolations from lattice results 
might suffer from fundamental problems. 



2.3.2 Positivity in an Euclidean quantum field theory 

In the Kugo-Ojima confinement scenario the conserved BRS-charge is employed to 
define a subspace Vphys from the state space of QCD, which can be shown to be positive 
semidefinite. However, this is only one particular mechanism to ensure the probabilistic 

^''As the quantum Maxwell equation 12.59|l can be written down for the general case a 7^ as well, we 
believe it to hold for general values of the gauge parameters A and a. 
^^This is done using the function U{p'^) which is defined as 

J d^xe'P^O\TD^c''{x)D,c\0)\0) = 

Here the symbol T denotes time ordering. It can be shown that the parameter u is given hy u — U {p^ = 0) 



(Sp.i' 



2 -)uip') - 2 



Sab 



(2.66) 
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interpretation of the quantum theory. Even if the Kugo and Ojima scenario eventually 
will turn out not to apply for QCD, there has to be some mechanism which singles out 
a physical, positive semidefinite subspace in QCD due to the arguments given at the 
beginning of subsection 12.3.11 

This suggests another, quite general criterion for confinement, namely violation of 
positivity. If a certain degree of freedom has negative norm contributions in its propagator 
it cannot be part of the physical content of the theory^^. In the following we will briefly 
explain what positivity means in the context of a Euclidian quantum field theory. In 
chapter ini we will search for negative norm contributions in our solutions for the quark and 
gluon propagators. 

Quantum field theories are completely described in terms of correlation functions, 
which can be ordered by an infinite hierarchy with respect to the number of contributing 
space-time points. Basic examples of such correlation functions are propagators (two point 
correlations) and vertices (three point correlations). Both are the central objects inves- 
tigated in this work. These correlation functions are subject to mathematical properties 
described by the axioms of axiomatic quantum field theory^''. A set of such axioms has 
been given first in 1957 by Wightman [20] for field theories formulated in Minkowski space. 
The Euclidean counterpart of the Wightman axioms have been found by Osterwalder and 
Schrader [HI E2] in 1973. 

This second set of axioms will be relevant for us, as we will work in Euclidean space- 
time throughout this thesis. There are two reasons for this choice: one technical and one 
physical. The technical aspect is the absence of poles on the real positive p^-axis for the 
propagators of Euclidean field theory. This is a huge practical simplification when dealing 
with Dyson- Schwinger equations^^. Maybe even more important is the physical aspect. In 
the course of this work we will run across several results which are most advantageously 
compared to lattice Monte Carlo simulations, which are performed in Euclidean space- 
time. Whereas both methods alone suffer from specific problems detailed in later chapters, 
the interplay between lattice simulations and Dyson-Schwinger calculations leads to well 
justified statements on the infrared behaviour of QCD correlation functions. 

The Euclidean counterpart to the notion of positivity in Minkowski space is the Oster- 
walder-Schrader axiom of reflection positivity. For a thorough mathematical formulation 
of the axiom the reader is referred to refs. [HZIEIIESI- We are interested in the special case 

^®If the Kugo-Ojima construction is appropriate, this imphes that the corresponding particle is a member 
of an asymptotic BRS-quartet. 

^^An introductory overview on axiomatic quantum field theory is given in the book of Haag |37|. 

^^Only some exploratory calculations of Dyson-Schwinger equations in Minkowski space can be found 
in the literature, see e.g. [^IKlj . 
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of a two point correlation function, i.e. a propagator S{x — y), for which the condition of 
reflection positivity can be written as 

J d\ d'y fix, -xo) S{x - y) f{y, y,) > . (2.67) 

where / is a complex valued test function with support in {(x, xq) : xq > 0} and time 
ordered arguments. After three-dimensional Fourier transformation this condition implies 

oo 

J dtdt' f{t\p)S{-{t + t'),p) f{tj) >0 (2.68) 



where S'(xo,p) := / d^x S{x,Xo) e*^^'. The momentum dependence of the corresponding 
Fourier transform of the test function / has been chosen to provide a suitable smearing 
around the three-momentum p. This condition will be tested for the quark and gluon 
propagator in chapter IHl 



2.3.3 The horizon condition 

The last aspect of confinement summarised in this section is the horizon condition formu- 
lated by Zwanziger [HSl IHE] • This brings us back to the gauge fixing procedure, described 
in section 12.11 Recall that the problem of fixing a gauge is equivalent to singling one 
representative configuration from each gauge orbit 

[A^] ■= [A^ = UAU^ + UdU^ : U{x) G SU{N^)] . (2.69) 

It has been shown by Gribov ^l^, that the simple Faddeev-Popov procedure generates a 
hyperplane F in gauge field configuration space which still contains gauge field configura- 
tions connected by a gauge transformation. These multiple intersection points of a gauge 
orbit with F are called Gribov copies. An almost unique representative of each gauge orbit 
is obtained, if one restricts the hyperplane F to the so called Gribov region Q. This is 
conveniently done by minimising the following L^-norm of the vector potential along the 
gauge orbit 

Fa{U) ■= \ \A^\ p = I \A\\^ -2i j d^x ti{udA) + j d^x ti{ujFP{A)uj) + 0(cj=^) , (2.70) 

with the gauge transformation U = exp{iuj{x)) and the Faddeev-Popov operator 

FP{A) = -dD{A)''^ = -8^6"^ - gT^^d^A^ . (2.71) 

Any local minimum thus implements strictly the Landau gauge condition dA = 0, and 
the Faddeev-Popov operator has to be a positive operator. The Gribov region fl defined 
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Figure 2.1: Sketch of the hyperplane F in gauge field configuration space obtained by the 
Faddeev-Popov gauge fixing method. Shown are the first Gribov region fi, and the fundamental 
modular region A containing the trivial configuration A = 0. The arrows pointing to the 
boundary OA of the fundamental modular region indicate that different boundary points have 
to be identified thus generating topological obstructions. 

by this prescription can be shown to be convex, to contain at least one intersection point 
with each gauge orbit and to be bounded in every direction of the hyperplane F [HH] . 
Furthermore the lowest eigenvalue of the Faddeev-Popov operator approaches zero at the 
boundary dQ, the first Gribov horizon. 

The absolute minimum of the Fa defines the fundamental modular region A. It can 
be shown that A is convex as well, bounded in Q and contains the origin A = at a 
finite distance of the boundary dA. Furthermore points on the boundary dA have to be 
identified to completely fix the gauge [HDl EHI- The set dfl fl dA contains the so called 
singular boundary points which are related by infinitesimal gauge transformations. This 
whole situation is sketched in figure I2.ll 

The horizon condition has been formulated in an attempt to restrict the generating 
functional of the gauge fixed theory to the fundamental modular region A. It was shown 
for lattice gauge theory in the thermodynamical limit [HSl IM] as well as for continuum 
theory [SEl, that such a restriction might be possible, if the probability distribution inside 
the fundamental modular region is concentrated at the Gribov horizon, i.e. at the region 
dfindA. Entropy arguments have been employed to argue for this condition. Furthermore 
it has been argued, that this implies the quantum field theory to be in the nontrivial, 
confining phase 
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Interestingly enough, the horizon condition originally formulated in different terms can 
be connected to the ghost dressing function G{p^) IHT ] IH^. which has been defined in 
the last subsection, eq. ()2.64j) . Due to the proximity of infrared {i.e. almost constant) 
gauge field configurations to the Gribov horizon jJT] the horizon condition is equivalent to 

\wi\G{p^)]-^ = . (2.72) 

This is the same condition as the Kugo-Ojima confinement criterion, eq. ()2.65jl . 

Furthermore, the same entropy arguments have been been employed to argue for a 
vanishing gluon propagator -D(p^) in the infrared [U^ 135]: 

\im\D{p')] = . (2.73) 

p— >o 

Both of these conditions can be checked by solving the corresponding Dyson-Schwinger 
equations (DSEs) for the ghost and gluon propagators in Landau gauge. Turning the 
argument round, eqs. ()2.72|) . ()2.73|) are appropriate boundary conditions for the DSEs to 
generate solutions corresponding to a restriction of the generating functional ()2.1j) to the 
Gribov region Q jn2]- In the remaining chapters of this thesis we will demonstrate the 
agreement of our solutions of the DSEs with the conditions ()2.72|) . ()2.73|) . 

2.4 The Dyson-Schwinger equations for the QCD prop- 
agators 

Having summarised some general aspects of strong QCD which are relevant for the dis- 
cussion of the propagators of the theory we focus our attention on the Dyson-Schwinger 
equations of motion for the ghost, gluon and quark propagator. These are a coupled set 
of integral equations which are derived from the generating functional ()2.9j) together with 
the generalised Lagrangian ()2.11|) . We concentrate on the derivation of the ghost DSE, as 
the four-ghost interaction term of our generalised Lagrangian generates new loops different 
from those of standard Faddeev- Popov gauges. As the derivation is rather lengthy we just 
give a summary in this chapter, deferring details to appendix El The formal structures of 
the gluon and quark DSEs remain unchanged compared to the Faddeev-Popov case and 
we therefore just give the results at the end of the section. 

With the conventions defined in appendix lA . 31 the Dyson-Schwinger equation of motion 
for the ghost propagator reads 

(^c''(i/)^ = <5(.-y)4.. (2.74) 
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Here the brackets (.) indicate the expectation value of the enclosed field operators. The 
derivative of the action S[J, c, c] = / d'^x C is given by 

= + f - f ) ^9'r''f^''-A^)c^i^)c'i^) 

+^ (l - f ) 9r'%. {Al{z)Az)) + z^gr'^^A;{z)d,c'{z). (2.75) 

We now decompose the full four-ghost correlation function into connected parts and use 
the relation 

S{y - = j d^z [D'^(^{z - y)] D^J{x - z) 
for the ghost propagator Dq to arrive at 



1-1 



- f - f ) \3'r''f^'' j [D-{x - z)r X 

{{<^{y)A^)cf{zY{z)) + {<^{yy{z)){A^)cf{^)) 

-{<*{y)cf{z)){c\zy{z))] 

- ^ (l - f ) ar'^ j d'z [Z^-(x - z)]-'{c\y)d, {Al^izyiz))) 

- ^^9r''J d'z [D^Jix - z)]-\c\y)A';,{z)d,c'{z)). (2.76) 

The remaining task is to decompose the connected Green's functions into one-particle 
irreducible ones. Plugging in the definitions of the bare ghost-gluon and the bare four- 
ghost vertex defined in appendix IA.3I we arrive at the ghost Dyson-Schwinger equation in 
coordinate space: 

[D^'{x-y)]-' = [DP^\x-yr' 

- J d'[zuvz^z,z,] rW^^^^d/, u, v) D^iiv - z,) rf ^(^1, Z3, x) D^Jiu - Zs) 

- I d'luv] rl°f ^"(x, u, V, y) D{f{v - u) 

~ \ j d'^[zuvuiU2UsU4] TfJ^'^'^^{y, z, v, u) D^^^u - U4) x 

D^civ - U2) Hg^l^s, X, U4, U2) D^(^{u3 - z) 

- \ j d'^[zuvuiU2UsU4U5] Tf^l^'^'^^{y,z,v,u) Df^{ui - u^) x 

Dfi{u - Ms) Tl-\u,,x,u^) D'^iv - U2) rf{uuU3,U2) 0^,%^ - z). 

(2.77) 
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Here we used the abbreviation d^[xyz\ := d'^x d'^y d'^z. The ghost-gluon vertex is denoted 
by Vy and the four-ghost vertex by T^^gh. The superscript (0) indicates a bare vertex or 
propagator. 

We now perform a Fourier transformation of eq. ()2.77|) and finally introduce renormal- 
isation factors at the appropriate places: 

- Zi^, I d'qT'^;\p,q)D,,{p-q)T,{q,p)DG{q) 

+ J d^[<li(l2]'r4gh Dciqi) Dg{p -qi- ^2) ^Agh{p, gi, 92) ^0(^2) 

Ty{p, qi) DG{q2) ^/.{-p + Qi + ^2, ^2) Dg{p - gi - 52) • 

(2.78) 



The colour traces have already been carried out and the reduced vertices defined in ap- 
pendix IA.3I have been used. The four-ghost interaction generates three new diagrams in 
the ghost equation, a tadpole contribution and two two-loop diagrams. Furthermore the 
bare ghost-gluon vertex depends on the gauge parameter a, 



gr'%2T,f5\k + q-p)V^;^\p,q) 



a\ a ' 
1 - 2 j + 2^^ 



(2.79) 



Note the symmetry between the ghost momentum p^ and the antighost momentum g^, 
when the gauge parameter a is set to one. 

The respective equation for the gluon propagator is formally the same as in the 
Faddeev- Popov case, where the equations have been derived in ref . [HI] • Differences occur 
in the explicit form of the bare ghost-gluon vertex and the dressed vertices in general 
depend on the gauge parameters. As the derivation of the DSE from the generating func- 
tional brings nothing new we refrain from displaying it explicitly and merely give the 
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result: 



1 g'^N^ f 

- ^4-^^ y d^lqm] rJIJJ^A ^PP'(?2) /^<xa'(P - 92 - ^l) X 

^p'ux'a'ip, qi,q2) D\\'{qi) 

- ^42^;^ J r(°U - - q,) D^^^q,) x 

rp'c<x'(p - gi - 92, 92) ^cc'(p - 9i) rcvA'(p - gi, gi) Dxy{qi). 

(2.80) 

The definitions and conventions for the gluon propagator D^^^, the three-gluon vertex F^^^p 
and the four-gluon vertex F^j^p^ are given in appendix IA.31 

The Dyson-Schwinger equation for the quark propagator is derived in a similar way 
from the generating functional and reads explicitly 

[S{p)]-' = [S^'\p)]''-Z,F^^ J d'qVf^p^q) D,^{p~ q) S{q)Tl{p,q) , (2.81) 

where S denotes the quark propagator and F^ the quark-gluon vertex. The factor Cf = 
[N^ — l)/{2Nc) stems from the colour trace which has already been carried out. 

All three Dyson-Schwinger equations are shown diagrammatically in figure 12.21 One 
clearly sees the striking similarity between the ghost and the gluon equation once a four- 
ghost interaction has been introduced. Both equations have bare and one loop parts, a 
tadpole contribution, a sunset and a squint diagram. 
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Figure 2.2: The coupled Dyson-Schwinger equations for the gluon, ghost and quark propaga- 
tors from a BRS and Anti-BRS symmetric Lagrangian. Each of the equations for the gluon and 
ghost propagators contains genuine one-loop diagrams, a tadpole contribution and a sunset 
and a squint diagram. All signs and weight factors have been absorbed in the diagrams. 



Chapter 3 



Propagators of Landau gauge 
Yang-Mills theory 

In this chapter we investigate the Dyson-Schwinger equations for the propagators of Yang- 
Mills theory. The knowledge of the two point functions of Yang-Mills theory, the ghost and 
gluon propagator, might shed light on those fundamental properties of QCD, which are 
generated in the gauge sector (for a recent review see JHl)- This is clearly the case for the 
phenomenon of confinement, as can be inferred from lattice calculations^. Furthermore the 
knowledge of the interaction strength in the gauge sector of QCD provides the basis for a 
successful description of hadronic physics [111110! • Based on the idea of infrared slavery older 
works on this subject assumed a gluon propagator that is strongly singular in the infrared. 
Recent studies based either on Dyson-Schwinger equations jHSl EH IHHl 1101 1121 HOI UHj or 
Monte-Carlo lattice calculations [ZH [Z21 IZHl jEEl IZll IZSl HO] in Landau gauge indicate quite 
the opposite: an infrared finite or even infrared vanishing gluon propagator. 

Lattice simulations and the Dyson-Schwinger approach are complementary in the fol- 
lowing sense: On the one hand, lattice calculations include all non-perturbative physics of 
Yang-Mills theories but cannot make definite statements about the far infrared due to the 
finite lattice volume. On the other hand, Dyson-Schwinger equations (DSEs) allow one 
to extract the leading infrared behaviour analytically and the general non-perturbative 
behaviour with moderate numerical effort. However, the infinite tower of coupled non- 

^In these simulations the string tension between a pair of static quarks has been calculated and found to 
be linearly rising as is expected for confined quarks. No dynamical quarks are involved in the calculations, 
therefore one concludes that the string tension is generated by the gauge field only. There is even evidence 
that only very particular gauge field configurations, namely center vortices, are responsible for the linearly 
rising potential |65| . This might shed new light on the origin of the violation of the cluster decomposition 
principle discussed in subsection 12.3. II 
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linear integral DSEs has to be truncated in order to be manageable. As we will see in 
the course of this chapter, the propagators of SU(2) and SU(3) Landau gauge Yang-Mills 
theory coincide for these two different approaches reasonably well. Thus we are confident 
that our results for the qualitative features of these propagators are trustworthy. 

Throughout this chapter we stay in the framework of ordinary Faddeev-Popov quan- 
tisation. Interesting enough, some of our results can be directly compared with recent 
calculations obtained in a framework employing stochastic quantisation JH] • We will thus 
be able to check for systematic errors connected to the appearance of Gribov copies^ [Hj , 
c.f. the discussion in section 1^. H.HI 

Landau gauge, which has been chosen for all DSE studies of Yang-Mills theory so far, 
is special for a number of reasons. First, it is a fixed point under the renormalisation 
procedure. This means that the gauge parameter A is not renormalised when A = 0, a 
fact which simplifies the renormalisation of the DSEs considerably. Second, as we saw in 
the last chapter. Landau gauge is a ghost-antighost symmetric gauge. This is of principal 
interest as we are then allowed to interpret ghost and antighost as (unphysical) particle 
and antiparticle. On the other hand it simplifies matters if one attempts to construct a 
non-perturbative dressed ghost-gluon vertex, as one is guided by a symmetry. This has 
been exploited in references 1^21113 ED]- Third, the ghost-gluon vertex does not suffer from 
ultraviolet divergences in Landau gauge, as has been shown by Taylor |HI]. Again, this 
simplifies the search for a suitable ansatz for this vertex considerably. Indeed, one is even 
allowed to use the bare ghost-gluon vertex, as we will see in the course of this chapter. 

This chapter is organised as follows: We first give a brief summary of previously em- 
ployed truncation and approximation schemes for the coupled gluon and ghost DSEs in 
Landau gauge. We discuss the key role of the ghost-gluon vertex in these truncations and 
show how a non-perturbative definition of the running coupling can be inferred. One of the 
obstacles encountered in providing numerical solutions are the angular integrals inherent 
to these equations. Therefore approximated treatments of the angular integrals have been 
applied so far IHTj . In general these angular approximations proved to be good for 
high momenta but less trustworthy in the infrared. Recent studies HOI therefore 

concentrated on the infrared analysis, where exact results have been gained for the limit of 
vanishing momentum. However, as we will see in the course of this chapter, not every ex- 
tracted infrared solution is connected to a numerical solution for finite values of momenta. 
The main part of this chapter is therefore devoted to the construction of a novel trunca- 
tion scheme, which allows to overcome the angular approximation for the whole range of 

corresponding investigation of the influence of Gribov copies on the gluon propagator in lattice 
simulations can be found in ref. |Z3 EHI 
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momenta [HHl IH^- Thus we are able to single out the physical infrared solutions of the 
schemes used in |H21 1021 HOI ■ In the ultraviolet region of momentum we obtain the correct 
one-loop anomalous dimensions of the propagators known from resummed perturbation 
theory. 

3.1 Gluon and ghost Dyson— Schwinger equations in 
flat Euclidean space-time 

The coupled set of gluon and ghost Dyson-Schwinger equations, which has been given 
diagrammatically in Fig. 12. 21 in the last chapter, loses a considerable amount of complexity 
in Landau gauge. There, the four-ghost vertex vanishes and we are left with one dressing 
loop in the ghost equation. As we will only be concerned with pure Yang-Mills theory 
in this chapter the quark loop in the gluon equation disappears as well. Furthermore 
in Landau gauge the tadpole term provides an (ultraviolet divergent) constant only and 
will drop out during renormalisation. Thus we will neglect this contribution from the 
very beginning. All these simplifications are due to our choice of Landau gauge. The 
resulting system of equations is still very complicated as it contains full two-loop diagrams 
in the gluon equation. The first assumption of all truncation schemes up to now is, that 
contributions from these two- loop diagrams may safely be neglected^. We will join in this 
assumption and provide some arguments for its validity in subsection 13.4.21 

Thus, we will effectively study the coupled system of equations as depicted in Fig. 13.11 
The corresponding equations are given by 

[Dcip)]-' = Z,[DP{p)]-'-Z,0^^Jd\T(^\p,q)D,Up-q)TAq,p)DG{q), (3.1) 
[D{p)];l = Z,[D^'\p)]].l + Z,^^jd'qTf{p,q)Da{p-q)TM,p)Da{q) 

-^i2^ J A r(°l(p, 9) ^pp'(p- 9) rpw'(g,p)i^..'(?)- (3.2) 

Here the ghost-gluon vertex is denoted by the symbol T^{q,p), whereas the three-gluon 
vertex is given by Tpi,„{q,p). Furthermore we have the coupling g and the number of 
colours Nc stemming from the colour trace of the respective loops. Suppressing colour 
indices the explicit expressions for the ghost and gluon propagators as well as the inverse 



■^The only exception is the scheme discussed in [HHl, which has not been solved yet. 
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Figure 3.1: Diagrammatic representation of the truncated Landau gauge gluon and ghost 
DSEs studied in this chapter. In the gluon Dyson-Schwinger equation terms with four-gluon 
vertices and quarks have been dismissed. 



of the gluon propagator are given by 



Dcip) = (3.3) 
n f \ fx Pt^PA ^iP^) , xPf^P-^ A\ 

d,Ap) = y^^--—r)—^ + ^-^^ (3-4) 

[D^]-^ = - + jp^p, . (3.5) 

For all linear covariant gauges the longitudinal parts of the full and bare inverse gluon 
propagators cancel each other in the gluon equation ()3.2|) . Furthermore in Landau gauge 
we have A = 0. 

At this stage of treating the equations we spot a problem: The left hand side of the 
gluon equation ()3.2j) is transverse to the gluon momentum therefore the right hand side of 
this equation should be transverse as well. This is certainly the case in the exact theory. 
However in an approximate treatment the gluon polarisation on the right hand side may 
acquire spurious longitudinal terms due to breaking gauge invariance. In general there are 
two possible sources for this violation: The first one is the use of vertices which violate 
the corresponding Slavnov- Taylor identity. The second one is the use of a regularisation 
scheme which breaks gauge invariance, such as a cutoff in the radial momentum integral. 

We postpone the problem of the gauge invariance of the vertex ansatz to the next 
subsection and discuss the regularisation problem first. In the following we outline some 
rather abstract arguments which will become more transparent in section 13. 3| where we 
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investigate the infrared and ultraviolet properties of our new truncation scheme in detail. 

In general a cutoff in the loop integrals can lead to quadratically ultraviolet divergent 
terms in the gluon equation. Such terms are scheme dependent and therefore unphysical. 
Furthermore they are highly ambiguous because they depend on the momentum routing 
in the loop integral. Unfortunately, a gauge invariant regularisation scheme avoiding these 
terms is hard to implement in Dyson-Schwinger studies^. 

It has been argued though [1211201, that quadratic divergences can occur only in that 
part of the right hand side of the equation which is proportional to the metric 6^^. There- 
fore an alternative procedure to avoid quadratic divergences is to contract the equation 
with the tensor {HHl 

n,,{p)=S,.-4^, (3.6) 

which is constructed such that TZ^uip) ^'^^ = 0. However, as has become obvious recently 
[THj . the use of the tensor (jH.fij) interferes with the infrared analysis of the coupled gluon- 
ghost system (see also ref. jHI] for a corresponding discussion in a much simpler truncation 
scheme) . 

In order to study this problem more carefully we will contract the Lorentz indices of 
eq. ()3.2|) with the one-parameter family of tensors 

^i9(p) = 5..-C^. (3.7) 

This allows us to interpolate continuously from the tensor ()3.6|) to the transversal one 
(with = 1). We will then encounter quadratic divergences proportional to the factor 
(4 — (^), which can be identified unambiguously and removed by hand. We will be able 
to show that this procedure restores the correct perturbative behaviour of the equations 
even with a finite cutoff A. 

Having removed all quadratic divergences we are then in a position to evaluate the 
remaining degree of breaking gauge invariance. As a completely transversal right hand 
side would be independent of C after contraction with the projector ()3.7|) . the variation of 
our solutions with is a measure for the influence of the artificial longitudinal terms on 
the right hand side of the equation. 



"'For the corresponding use of dimensional regularisation see e.g. refs. jHHl EZI EHl ) • 



Chapter 3. Propagators of Landau gauge Yang-Mills theory 



33 




Figure 3.2: A diagrammatical representation of the ghost-gluon vertex. 

3.2 The ghost-gluon vertex in DSE studies and the 
running coupling 

We now focus on two issues that are connected with the ghost-gluon vertex of Landau 
gauge. We will argue for the surprising fact, that one can safely use a bare ghost-gluon 
vertex even in the non-perturbative momentum region of the DSEs. Furthermore we 
describe, how one is able to relate the running coupling of the strong interaction to the 
ghost and gluon dressing functions by the renormalisation properties of the ghost-gluon 
vertex. 



3.2.1 Dressing the ghost-gluon vertex 

For our further treatment of the ghost and gluon system of equations, ()3.ip and ()3.2|) . we 
have to specify explicit forms for the dressed ghost-gluon vertex T^{q,p) and the dressed 
three-gluon vertex Tpi^f^{q,p). As has already been mentioned in the introduction to this 
chapter, the ghost-gluon vertex does not attribute an independent ultraviolet divergence 
in Landau gauge, i.e. one has Z\ = \ ^HT. Therefore a truncation based on the tree-level 
form for the ghost-gluon vertex function, 

^M^V)=np. (3.8) 

is compatible with the desired short distance behaviour of the solutions. Here the momen- 
tum is the momentum of the outgoing ghost, see Fig. 13.21 Thus we obtain the correct 
ultraviolet behaviour of the ghost loop in the gluon equation ()3.2p and the dressing loop 
in the ghost equation ()3.1|) . as will be shown explicitly below. 

However, as the effects of non-perturbative vertex dressing are supposed to be most 
pronounced in the infrared, one might wonder whether the tree-level form of the ghost- 
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gluon vertex leads to a sensible infrared behaviour of these equations. Furthermore the 
bare vertex violates the Slavnov- Taylor identity (STI), which restricts that part of the 
ghost-gluon vertex which is longitudinal in the gluon momentum. Such an identity is 
the manifestation of gauge invariance and can be derived using the BRS-invariance of the 
gauge-fixed Lagrangian of Yang-Mills theory. 

Considering this, obviously the best way to obtain a properly dressed ansatz for the 
ghost-gluon vertex is to solve the corresponding STI. This strategy has been followed by 
von Smekal, Hauck and Alkofer in |M1 ES] • On the level of connected Green's functions 
the STI for the ghost-gluon vertex in general linear covariant gauges reads [7^1 

j{c^{z)c\y)dA-{x)) - j{c^{z)c'^ix)dA\y)) = -^-f"^ {c\z)c\z)-c\x)c\y)) . (3.9) 

Here we have the spatial coordinates x, y and z, the bare coupling the gauge parameter 
A and the real structure constant /'^'''^ of the gauge group SU(3). The left hand side of this 
equation can be decomposed into the full ghost-gluon vertex and respective propagators. 
However, the right hand side contains the connected ghost-ghost scattering kernel, which 
is completely unknown. 

By neglecting this irreducible correlation von Smekal, Hauck and Alkofer were able to 
construct a vertex ansatz which solves the resulting approximate STI. Together with a 
similar construction for the three-gluon vertex they obtained a closed system of equations 
which has been solved numerically using an angular approximation (c.f. subsection 15. l.ip . 
The results, an infrared vanishing gluon propagator and an highly singular ghost in the 
infrared have been confirmed since in other DSE-calculations as well as lattice Monte-Carlo 
simulations [7TJ [721 112 ESI- Thus the old idea of infrared slavery based on the notion of 
an infrared divergent gluon propagator has been abandoned since. 

However, as became clear later, the vertex construction of refs. |66t 179] is somewhat 
problematic. It has been shown in ref. that this vertex causes inconsistencies in the 
infrared behaviour of the ghost equation once the angular integrals of the dressing loop 
are treated exactly. Furthermore it has been shown in refs. 1^^] that the neglection of 
the irreducible ghost-ghost scattering kernel in the identity ()3.9p is at odds with pertur- 
bation theory. On the other hand it is hard to see, how one could include this irreducible 
correlation and thus improve the construction of refs. I79j. 

Therefore Atkinson and Bloch chose a different strategy and employed a bare ghost- 
gluon vertex in their truncation scheme jHTl |H2] , which keeps only the ghost loop in the 
gluon equation^. The numerical calculations in this scheme are also obtained using angular 

^Although there is an attempt in ref. [S7] to include the gluon loop as well, the authors themselves 
note an inconsistency between their construction and perturbation theory. Thus the focus was mainly on 
the 'ghost-loop only' situation. 
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approximations in the integrals. Amazingly, though, the bare vertex scheme and the 
one from EH] provide results with qualitatively similar infrared behaviour: the gluon 
propagator vanishes in the infrared and the ghost propagator is highly singular there. 

The surprising conclusion from the comparison of these two truncation schemes is, 
that a bare ghost-gluon vertex ()3.8|) is not only capable of providing the correct ultraviolet 
behaviour of the ghost loop in the gluon equation ()3.2|) and the dressing loop in the ghost 
equation ()3.H1 . but in addition leads to a satisfactory infrared behaviour of the equations in 
accordance with lattice Monte-Carlo simulations. Consequently recent analytical infrared 
investigations concentrated on the bare vertex [H21 IBS] • The influence of multiplicative 
corrections to the bare vertex has been assessed in refs. [701 MM found to be 

irrelevant for the qualitative behaviour in the infrared. 

For our novel truncation scheme, detailed in section 13. 3| we will therefore use the 
bare ghost-gluon vertex ()3.8p . keeping in mind that we have to check for the influence of 
artificial longitudinal terms due to the violation of the Slavnov- Taylor identity ()3.9j) . As 
a major improvement to previous calculations we will be able to overcome the angular 
approximation and give solutions for the ghost and gluon dressing functions which include 
the full angular dependence of the loops in the equations. 



3.2.2 The running coupling of the strong interaction 

Before we give the details of our truncation scheme there is still more to be learned from 
the ghost-gluon vertex in Landau gauge. In the following we consider the renormalisation 



of the vertex term in the unrenormalised Lagrangian, eq. (j2]TTj). In Landau gauge we have 

Cghost-gluon = gr'^C^A^c' , (3.10) 

which is identical to the respective term in the Faddeev- Popov Lagrangian, eq. ()2.1()j) . 
Multiplicative renormalisability means that it is possible to render all Green's functions 
finite by renormalising the fields and parameters of the Lagrangian without changing its 
form. Recall that the coupling g, the gluon field and the ghost field are renormalised 
according to (c./. eq. ()2.13j) ) 

A« ^ v^a;:, (3.11) 

c" \fKc\ (3.12) 

c'^ ^ (3.13) 
9 ^ Zgg, (3.14) 
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where we have unrenormahsed objects on the left hand side and renormahsed ones of the 
right hand side of the relations. Furthermore, by definition, the vertex renormalisation 
constant is related to the renormalisation constants of the constituent fields of the 
vertex by 

Zi = Zg~Z^Z'J''. (3.15) 
We then have the renormalised ghost-gluon vertex part of the Lagrangian 

^fhost-gluon = ^iQ f°'^'^dfjC°' A'j^^C^ , (3.16) 

which has the same form as the respective term ()3.10|) in the bare Lagrangian. 

We are now able to exploit the fact again that Zi = 1 in Landau gauge [SI]. As the 
strong running couphng is defined by a = 5f^/(47r), we obtain the relation 

«(A2) = ^ ?^ (3.17) 

zl{^^\K^)z,{^^\^^) 

from the renormalisation of the coupling, eq. ()3.14|1 . Here we gave the explicit arguments 
of the renormalised coupling a(/i^), evaluated at the renormalisation point /x, and the bare 
coupling a(A^), which depends on the cutoff A of our theory. From the relations ()3.11|) 
and ()3.13|1 we furthermore infer that the ghost and gluon dressing functions, GijP') and 
ZIp^), are renormalised according to 

G{p',A') = G(p^^2)Z3(/i^A2), 

Z{p',A^) = Z{p\fi^)Zs{fi',A'), (3.18) 

where the unrenormahsed quantities are on the left hand side of the equations. Certainly 
the renormalisation point /i is completely free. We are thus allowed to substitute these 
relations into eq. ()3.17p and renormalise once at an arbitrary point /i and once at the 
specific value p. We obtain 

a{A') G\p\A^) Z{p\A^) = a{i?) G\p',n') Z{p\^?) = a{p') G\p\p') Z{p\p^) 

(3.19) 

It suffices now to impose the non-perturbative renormalisation condition 

G\p\p')Z{p\p') = G\^^\^?)Z{^^\^^) = 1 (3.20) 
on equation ()3.19|1 to arrive at 



a{p')=aif^')G\p^fi')Zip\fi') 



(3.21) 



This is a defining equation for the non-perturbative running coupling a(p^) of Landau 
gauge QCD [^11 ESI- Note that the running coupling a(p^) defined this way does not 
depend on the arbitrary renormalisation point /x. This is equivalent to saying that the 
right hand side of this equation is a renormalisation group invariant [02] ■ 
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3.3 The truncation scheme 

Having spent some time with the discussion of Landau gauge we now return to the coupled 
set of eqs. fl3.1|) and ()3.2|) and specify exphcit expressions for the ghost-gluon vertex T^{q, p) 
and the three-gluon vertex Tp^(j{q,p). In the following we will show, that the bare ghost- 
gluon vertex 

T^.iq,p) = ^q>. (3-22) 

and the construction 

P , , _ 1 G{{q - pf)^^-^l^-^) ^ 

^pua{q,P) z^{ii\K^) 'p-^^'^^ ^'^'^'^^ 

with the bare three-gluon vertex vf^a given in eq. ()A.3H) and the new parameters a and h 
lead to the correct one-loop anomalous dimensions of the dressing functions in the ultra- 
violet, provided the quadratic divergences have been removed. This is true for arbitrary 
values of the parameters a and 6, although we will later argue for the specific values 
a = b = 36, where 6 is the anomalous dimension of the ghost. 

For simplicity we introduce the abbreviations x '■= p"^, y '■= and z := k"^ = {q — p)"^ 
for the squared momenta appearing as arguments of the dressing functions. Furthermore 
s := /i^ and L := denote the squared renormalisation point and the squared momentum 
cutoff of the theory. Substituting the two vertices in the ghost and gluon system ()3.1|1 and 
fl3.2|l we then arrive at 

' - ' g^Nj ^,^^i^Giy)Ziz) , (3.24) 



G{x) ^ 7 (27r)4 xy 

1 ry , 2^c [ A M{x,y,z) 

- Z^ + g ^ TTTTZ G[y)G[z) 



Z{x) S J (27r)4 xy 



^N, f d\ Q{x,y,z) GiyY'-''/^-^''^ G{zY'-'/^-^''^ 



+^ 3 J (27r)4 xy Z{yY Z{zf ' ^^"^^^ 



The kernels ordered with respect to powers of z have the form: 
K{x,y,z) -- ' 



z^ \ 4 J z \ 2 J 4' 

1 y Cy^\ 1 Cy Cz 

I flx^ 2 19-C 5-C 2 Cl/^ 

Q{x,y,z) = — xy + ——y + - — 

z^ \8 y 8 4 8x 

1 fx' 15 + C 17-C ^,y' 

H — — X — y + C— 

z \y 4 4 X 

— 7, \ A ^ — + -2 - H — + -2 7, — • (3.26) 

8 y 4 4 X J \x Ay J 8xy 
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3.3.1 Ultraviolet analysis 



In order to identify the quadratically divergent terms in the kernels K, M and Q we now 
analyse eqs. ()3.24|) and ()3.25|) in the limit of large momenta x. It is known from resummed 
perturbation theory (see e.g. [H]) that the behaviour of the dressing functions for large 
Euclidean momenta can be described as 



Z{x) 
G{x) 



Z{s) 
G{s) 



UJ log 
u log 



+ 1 
+ 1 



1 7 



(3.27) 
(3.28) 



Here Z{s) and G{s) denote the value of the dressing functions at some renormalisation 
point s := and 7 and 6 are the respective anomalous dimensions. To one loop one has 
6 = —9/44 and 7 = —1 — 26 = —13/22 for arbitrary number of colours and no quarks, 
Nf = 0. Furthermore, uj = llNca{s)/127i'. 

The slowly varying logarithmic behaviour of the dressing functions in the ultraviolet 
justifies the angular approximation 



Ziz),G{z)^Ziy),Giy) for y>x 



(3.29) 



in the ultraviolet. The angular integrals in eqs. ()3.24|1 and ()3.25j] can then be trivially 
calculated using the angular integration formulae of appendix IC.ll Furthermore, as the 
cutoff L = h? can be chosen arbitrary large, the integrals will be dominated by the part 
from X to L. We then obtain 
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(3.30) 



. (3.31) 



We are now able to identify the quadratic divergences in the gluon equation, which are 
the two terms independent of the integration momentum. Both, the one from the ghost 
loop and the one from the gluon loop, are proportional to (4 — Q and therefore vanish 
when we use the Brown- Pennington projector, eq. ()3.(i|l . For general values of C, we have 
to subtract these terms by hand. However, this cannot be done straightforwardly at the 
level of integrands: Such a procedure would disturb the infrared properties of the Dyson- 
Schwinger equations. Since we anticipate from previous studies and analytic work jH21 HOI 
that the ghost loop is the leading contribution in the infrared the natural place to subtract 
the quadratically ultraviolet divergent constant is the gluon loop. We do this by employing 
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the substitution 



Q{x,y,z) Q{x,y,z) = Q{x,y,z) + -{4- C) 



(3.32) 



in eq. ()3.25p . Only the problematic terms in eq. ()3.31|) then disappear. 

Having subtracted the quadratic divergences from the ghost and gluon system we now 
check for the logarithmic divergences, which should match to perturbation theory. We 
choose the perturbative renormalisation condition G{s) = Z{s) = 1 at a large Euclidean 
renormalisation point s = /i^ and plug the perturbative expressions ()3.27|) and ()3.28p in 
eqs. ()3.3()j) . ()3.31|1 . Thus we arrive at 
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(3.33) 



Note that these equations are completely independent of the parameters a and b in the 
three-gluon vertex, eq. ()3.23p . Note also, that the ultraviolet behaviour of the equations 
is independent of the parameter ( in the projector ()3.7p . We are therefore left with a 
transversal structure in the gluon equation for ultraviolet momenta. The renormalisation 
constants Z^^s, L) and Z^{s, L) cancel the cutoff dependence, i.e. the respective first terms 
in the brackets. Thus, the power and the prefactor of the second term have to match with 
the left hand side of the equations. This leads to three conditions: 



7 + 2(5 + 1 
3 N,g^ 
4(7 + (5 + l)167r2c<; 
13 N,g^ 
6 {25 + 1) IGtt^cu 



0, 
1, 

1 . 



(3.34) 
(3.35) 

(3.36) 



Eq. ()3.34|) is of course nothing else but consistency of the ghost equation with one-loop 
scaling. All three equations together result in the correct anomalous dimensions 7 = 
— 13/22 and 5 = —9/44 for an arbitrary number of colours and zero flavours. 

Having established this result let us pause for a moment and reflect our construction 
fl3.23|l for the three-gluon vertex. There are two possible ways to obtain this construction 
and we leave it to the taste of the reader which philosophy to prefer. The first way to 
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see eq. (I3.23|l is to take it at face value as a minimally dressed vertex ansatz, constructed 
in such a way as to obtain the correct one-loop scaling of the gluon loop. Then the 
parameters a and b are completely free. Note that the choice a = b = corresponds to 
the truncation scheme of [531 whereas a = 36,b = together with the appropriate vertex 
dressings reproduces case c) of ref. [66] . 

The second point of view, which is to be found in reference jnEj, is to employ a bare 
three-gluon vertex, r°^^(g,p), and then ask the question how the renormalisation con- 
stant Zi has to behave, given that the theory should have the correct perturbative limit. 
Certainly, as both vertices violate their respective Slavnov- Taylor identities, there is no 
reason why the vertex renormalisation constant Zi should obey the corresponding iden- 
tity Zi = Z3/Z3. For the present truncation scheme the answer is, that the constant Zi 
acquires a momentum dependence according to 

Z, Z,ix,y,z;s) = ^^^^^ " ^^'^^^ 

This is precisely the form required to transform our vertex ()3.23|1 to the bare one and 
nevertheless obtain eqs. ()3.24|) and ()3.25j) . 

A reasonable choice of parameters is then one which keeps Zi as weakly dependent as 
possible on the momenta y and z, c.f. Fig. I(y.2l in appendix 1(^.51 The infrared behaviour 
of the gluon loop in the gluon equation depends strongly on a and b. Setting 6 = one 
can distinguish three cases: For a < the gluon loop is subleading in the infrared, for 
a = as in ref. jHS] the gluon loop produces the same power as the ghost loop, for a > 
the gluon loop becomes the leading term in the infrared. In the latter case we did not 
find a solution to the coupled gluon-ghost system. In Appendix l( 1 51 we demonstrate that 
a = b = 35 minimises the momentum dependence of Zi. Thus we use these values except 
stated otherwise explicitly. 

3.3.2 Infrared analysis 

The leading infrared behaviour of the propagator functions in this truncation scheme 
for the special case of the transverse projector {( = 1) has been determined recently 
[HS IZOl- Our analysis in this subsection is valid for general values of the parameter 
( and furthermore includes subleading contributions The general assumption at 

the beginning of all analytic infrared investigations is, that the ghost and gluon dressing 
functions, G and Z, behave like power laws in the infrared: 



Z{x) = Ax"', 
G{x) = Bx''\ 



(3.38) 
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This assumption is justified, if we are able to sfiow tliat tlie ansatz ()3.38j] solves tlie ghost 
and gluon system, eqs. ()3.24j) . ()3.25|1 . self-consistently in the infrared. As all loop integrals 
in the equations are dominated by contributions around the external momentum x = p"^, 
we are allowed to substitute the power laws for the whole momentum range up to the cutoff 
L = A^. Errors due to this approximation are subleading in the infrared. Furthermore, 
as has been shown in detail in reference [ZHI, the renormalisation constants and 
can be dropped for very small momenta x: They are either subleading in the infrared 
(gluon equation) or have to be zero when the renormalisation takes place at /i = (ghost 
equation). We will explain this point more precisely later on. Plugging the power laws 
into eqs. ()3.24j) . ()3.25|1 we then arrive at 

-/^^cA5 (3.39) 



5x^2 J (27r)4 xy 



Ax'^^ 3 J (27r)4 xy 

_^^2j^2J-4-12<5^-6<5 f C^^g Q{x, y, z) ^(_2_6^)k2-3^ki^(-2-65)k2-33ki ('349) 

3 J (27r)4 xy 

with the kernels M and Q given by eqs. ()3.26p . ()3.32|) and the parameters a = 6 = 35 
as motivated in the last subsection. 

We first investigate the ghost equation ()3.39j) . Shifting the cutoff L = to infinity 
we are able to use the formula [ZOI 

i', J,"/ = A-^-.^ r(2 + a)r(2 + 6)r(-a-t-2) 

T{-a)T{-h)T{A + a + h) ^ ' 

for the integration of the dressing loop on the right hand side of the equation. The 
straightforward but tedious algebra is done with the help of the algebraic manipulation 
program FORM [HE]- We obtain 



/ 



1 3 



Bx""-^ 167r2 2 (ki + K2) (-1 + Ki + K2) 

r(3 + /€i + /t2)r(2-fi:i)r(i-/€2)' ' 

Matching powers of x on both sides of the equation we arrive at the condition 

Ki = -2/t2 . (3.43) 
With the definition k := —K2 we thus have from eqs. ()3.38|) the power laws 

(3.44) 
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Our derivation of these power laws used the special form of a bare ghost-gluon vertex. It 
has been shown, however, that relation ()H.4Hj] holds under the general assumption that the 
ghost-gluon vertex can be expanded in a power series We thus have the interesting 
situation that the infrared divergence of one of our dressing functions, Z or G, is always 
connected to the vanishing of the other. This has interesting consequences for the running 
coupling a{x). Recall the definition 

= a(s) G'^(x, s) s) (3.45) 

from eq. ()3.21|) with x = and s = fi'^ and a{s) = g^/4:'7i. Substituting the power laws 
fimil we obtain 

a{0) = ^AB^. (3.46) 

4:71 

The coupling approaches a finite value in the infrared. Within the framework of Dyson- 
Schwinger studies such a behaviour of the running coupling has first been found in ref. [HHj . 
An infrared fixed point is also found in the flow equation study in ref. jHZI and analytic 
perturbation theory [HHl EHl 11001 IM] • 

Consider now the gluon equation ()3.40j) . After some algebra using formula ()3.4H1 we 
find, that the ghost loop is proportional to x^**^, whereas the gluon loop is subleading in 
the infrared. We thus again obtain relation ()3.43|) . ni = — 2/t2, as matching condition for 
the ghost loop and the leading power on the left hand side of the equation in the infrared. 
Together with eq. ()3.42|) we arrive at the two conditions 



1 {2 + k){1 + k) T\2-k)T{2k) g^N, 



(3 - 2k) r(4 - 2K)r2(l + k) 487r2 

Ak-2 r2(2 - k)T{2k) g^N^ 



AB^ , (3.47) 
AB^ . (3.48) 



ACk-Ak + Q- 3C r(4 - 2^)12(1 + k) 487r2 

from the ghost and gluon equations in the infrared. Equating both left and right hand 
sides we are able to determine k, see Fig. 13.31 For the Brown-Pennington projector, i.e. 
C = A, one then finds the known solution k = 1 [nZj. However, as can be seen immediately, 
the left hand side of the second equation possesses a zero for k = 1/2 which is cancelled 
by a pole only for ( = A. Lowering ( only slightly a further solution with k slightly larger 
than 0.5 exists. For the transverse projector, i.e. ( = 1, this latter solution becomes 



^ ^ 93 V120T ^ Q 595353 ^3 
98 

in accordance with refs. 170] . Also the solution k = 1 changes continuously when 
lowering (. The corresponding k are then all larger than 1 and contradict the masslessness 
condition, see section 5.2 of ref. for a discussion of this condition. One of the main 
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FaddeevPopovgauges: a=0,2 




Figure 3.3: Here the graphical solution to equations (13.4711 and (13.4811 is shown. The thick 
line represents the left hand side of equation (I3.47j) . whereas the other four curves depict the 
left hand side of equation p.48jl for different values of the parameter (. The ellipse marks the 
bulk of solutions between k = 0.5 and k = 0.6 for different (, whereas the circles show the 
movement of the solution k = 1.3 for a transverse projector to k = 1 for the Brown-Pennington 
case, C = 4. 

results of the numerical section of this chapter is that the infrared behaviour k, ^ 0.5 
matches to a numerical solution whereas no numerical solutions could be found with the 
infrared behaviour k > 1. This is in perfect accordance with the Kugo-Ojima criterion 
and Zwanziger's horizon condition discussed in section 

Now let us come back again to the running coupling, eq. ()3.45|) . As can be seen directly 
from eqs. ()3.47j) . ()3.48|1 . the product Ncq'^AB'^ is constant for given k. With a(0) = 
g'^AB^/AiT one concludes immediately that a{x) is proportional to A*""^. Furthermore, 
the ghost and gluon dressing functions Z{x) and G{x) are independent of the number of 
colours: enters the Dyson-Schwinger equations only in the combination g'^N^ at our 
level of truncation. From the solution k = 0.595 of the infrared analysis with the transverse 
projector ( = 1 one determines the infrared fixed point of the running coupling to be 

= 3iv^ n(2 - r(2K) " '-'^'/^^ ^'-'^^ 

for general numbers of colours Nc- 
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Having determined the leading infrared behaviour of the ghost and gluon system we 
now check for subleading contributions, and judge their importance for the numerical 
treatment of the coupled system of equations. An obvious extension of the power law 
fl3.44|) is the ansatz 

Z{x) = Ax^'^il + fxP) 

G{x) = Bx-^{l+gxP) (3.51) 

which is substituted in the eqs. (I3.24|l and ()3.25j) . After integration the conditions on the 
leading term remain unchanged. Matching subleading powers leads to the coupled set of 
homogeneous equations for /, g and p: 

'^(k + p - 2) (/« + p - 1) (k + p) r(4 - 2k) r(i + k - p) r(3 + K + p) " )^ 

3 (-2 + 2/t + p) r(2 - h) r(2^ + p + 1) r(3 - /t - p) \ ^ 

V ''2{k + p-2){k + p-1){k + p) r(3 -2k- p) T{1 + K)T{'i + k + p) ) ^ 



A^K - 4k + 2p - 2Cp + 6 - 3a r(2 - k) T{2k - p)T{2- k + p) 
' 2k- p-l r(l + k)T{1 + k- p) r(4 - 2fi; + p) 



9 + f 



(3.52) 



Here v = Ncq'^ AB"^ / ASn"^ = a(0)/47r. There is either the trivial solution f = g = ot one 
has to set the determinant of these linear equations to zero. We then obtain the results 
P(i) = 0,p(2) = 0.58377, p(3) = 1.20300 and several other solutions with higher values of 
the power. The solution p(i) = corresponds to the pure power solution. The lowest 
non- vanishing solution, p(2) = 0.58377, is sufficiently high that we safely may neglect it in 
the numerical treatment of the infrared part of the equations. This will be detailed in the 
next section. 

Finally, let us come back to the role of the renormalisation constants Z3 and in 
eq. ()3.39j) . ()3.4()|1 . There are two possible situations, which have been clarified in 7Uj. 
First, consider an infrared vanishing dressing function, as is the case in the gluon equation. 
Then the left hand side of eq. ()3.40p as well as the loop integral on the right hand side 
are diverging, and the x-independent constant is subleading and therefore negligible in 
the infrared. Second, in the ghost equation we have the situation that the ghost dressing 
function G{x) diverges in the infrared. Therefore the equation is not renormalisable at the 
point /i^ = 0. Correspondingly the behaviour of the renormalisation constant Z^{p'^,h?) 
is such that for a given cutoff the renormalisation constant tends to zero as /i^ ^ 0. 
At the point = we have Z3(0,A^) = 0, the renormalisation process breaks down 
and no scale can be generated. The power solution ()3.44j] is then not only an infrared 
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approximation but a solution for the whole momentum range up to infinity. This in turn 
implies that the power solution can be determined with Z3 set to zero. We conclude that 
both renormalisation constants, Z3 and Z3, play no role in the infrared analysis of the 
ghost and gluon system in accordance with our assumption above eqs. ()3.40p . ()3.39|) . 



3.4 Renormalisation and numerical results 

Corresponding to the analytical solution of the equations ()3.24|1 and ()3.25|1 in the infrared, 
our numerical treatment is done without the help of any angular approximations. Thus 
the present investigation is the first calculation of the ghost and gluon dressing functions 
which takes into account the full angular dependence in the loops for all momenta [S3]. 
The technical details of the necessary numerical procedures to solve eqs. ()3.24|1 . ()3.25|1 
are given in appendix 1(141 In the following subsections we describe the renormalisation 
scheme employed in our calculations and give the numerical solutions for the ghost and 
gluon propagators as well as the running coupling. 



3.4.1 The renormalisation scheme 

We apply a MOM regularisation scheme similar to the ones used previously in refs. jMll67j . 
In such a scheme the equations for the ghost and gluon dressing functions for the external 
momentum x and a fixed subtraction scale tghost,tgiue are subtracted from each other. If 
we write the equations ()3.24|) and ()3.25|) symbolically as 

1 

Z3 + Ughostix) , (3.53) 
Z, + UgiUx), (3.54) 



G{x) 
1 

Z{x) 

this procedure yields 

1 1 



G{x) G(tghost 

1 1 



Z{X) Z{tgU 



- + '^ghostix) — '^ghostit ghost) ; (3.55) 
+ ^glue{x) — ^ghostitglue)- (3.56) 



We now see that the unknown renormalisation constants Z^ and Z3 drop out and instead 
of them we have to specify the two input variables G(tghost) and Z(tgiue)- For numerical 
reasons it is favourable to subtract the ghost equation at a very small momentum, tghost 
and the gluon equation at a perturbative scale /. We thus have to specify values for Z{1) 
and the parameter B in the power law ()3.44j] which describes the infrared behaviour of 
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the ghost dressing function. In practice we encounter a smooth transition at the infrared 
matching point (c.f. appendix l(I4|l of our numerical integrals only if B and Z{1) are 
uniquely related. This corresponds to the fact that we are not able to implement the 
perturbative renormalisation condition = Z{^'^) = 1 for general renormalisation 

points n but only the weaker condition Z{^'^) G'^{fi'^) = 1. A similar observation has been 
made in [011 HO]- In the actual calculation we use the arbitrary value Z{1 = (174 GeV)^) = 
0.83. The value of the renormalisation point fi is given implicitly by specifying the coupling 
a{fi^) = ^7l67r2 = 0.97 entering eqs. (jT^ . (HT^ . 

3.4.2 Numerical results 

For our numerical results we have fixed the momentum scale by calculating the running 
coupling for the colour group SU(3), and requiring the experimental value a{x) = 0.118 
at a; = M| = (91.187 GeV)^ [TI]^ . 

In Fig. I3.4l we show the results for the transverse projector, C = 1- The gluon and ghost 
dressing functions behave power-like for low momenta with k = 0.595 and obey one-loop 
scaling in the ultraviolet as expected. We found no numerical solution for the second 
analytical infrared power from subsection 13.3.21 k ^ 1.3. We thus conclude, that only 
the infrared power k = 0.595 connects to the numerical solution for all momenta. This 
result will be corroborated in chapter O As has been mentioned in the introduction to 
this chapter first results have been obtained in stochastic quantisation, which is supposed 
to solve the Gribov problem pB: . The infrared analysis of a bare vertex truncation scheme 
in stochastic Landau gauge leads to the power n = 0.521 for the transversal part of the 
gluon propagator. This is surprisingly close to the value in Faddeev-Popov quantisation 
and suggests the influence of Gribov copies within the Gribov region to be small. 

According to the power solution of the dressing functions the running coupling has a 
fixed point in the infrared. Furthermore it decreases logarithmically in the perturbative 
regime above several GeV in accordance with perturbation theory. For intermediate mo- 
menta the behaviour of our running coupling induces a node in the /5-function around 
(lOOMeV)^. This, however, corresponds to a double valued /3-function, which is a some- 
what strange result not expected from renormalisation group analysis. Therefore we regard 
it as an artefact of our truncation scheme. Recall from section I^TD that we have omitted 
the two-loop diagrams in the gluon equation in all our calculations. Since the two-loop 
diagrams contain gluon dressing functions only they are subleading in the infrared in the 
present framework. Furthermore they are subleading in the ultraviolet as can be seen 
from perturbation theory. Thus the only region where the inclusion of the two-loop dia- 
grams could lead to qualitative corrections of our solutions is the intermediate momentum 
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Figure 3.4: Shown are the results for the gluon dressing function Z , the ghost dressing function 
G and the running coupling a using a transverse projector, C = 1- The two sets of fit functions 
are given in eqs. (13.5711 and (I3.58j) . The /^-function corresponding to our DSE-solution is 
compared to the one- and two-loop expressions as well as to a polynomial in a. 



regime. This is the region where the bump in the couphng appears. 

The /3-functions in the bottom right diagram are defined in appendix IA.2I Compared 
to the one- and two- loop /5-functions of perturbation theory the /5-function from our DSE- 
solution |lU3j resembles the scaling behaviour of the one-loop result in the ultraviolet, that 
is for small values of a. 

The asymptotic behaviour of the solutions can also be seen from the functional form 
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of our fits. We employ two different fit functions |l()4j for the running coupling a{x): 

Fit A: a(x) = -. — --^ — '^^^^^ — — — -p^ rr-r , (3.57) 

ln[e + a,{x/AlcDY' + ^i(^/Aqcd)'^] 



Fit B: a(x) 



1 




-2 7 I (x/^QCd) 



(3.58) 



The value a(0) = 8.915/Nc is known from the infrared analysis. In both fits the ultraviolet 
behaviour of the solution fixes the scale, A = 0.714 GeV. Note that we have employed 
a MOM scheme, and thus Aqcd has to be interpreted as a'^Iqmi ^'^^^ scale has 
the expected magnitude. Fit A employs the four additional parameters: ai = 1.106, 
aa = 2.324, bi = 0.004, 62 = 3.169. Fit B has only two free parameters: a = 1.020, 
b = 1.052. The dressing functions Z{x) and G{x) are then described by 



R{x) 



c{x/AlcDr + d{x/Alcj,r 
l + c{x/AlcDY + d{x/AlcDY^' 

1+25 



a{x) 



a{jj,) 



Z(x) = ( 1 R^(x 



-5 



G(x) = (^) R-\x), (3.59) 

where c,d are fitting parameters for the auxiliary function R{x). They are given by 
c = 1.269 and d = 2.105. Recall from subsection 13.3.1! that the anomalous dimension 7 
of the gluon is related to the anomalous dimension 6 of the ghost by 7 = —1 — 25 and 
6 = —9/44 for the number of flavours Nf = 0. 

Whereas Fit A is better in the region 0.3 GeV^ < X < 1 GeV^ where a is strongly rising. 
Fit B is slightly better in the region 1 GeV^ < x < 10 GeV^, where hadronic r-decay takes 
place jl05| ll()6j . As can be seen in Fig. 13.41 both fits works very well and can be used as 
input for phenomenological calculations in future work. 

Our results for different values of the parameter ( {c.f. eq. ()3.7|l ) are shown in Fig. 13.51 
In accordance with the infrared analysis the power k, changes from k = 0.5953 for C = 1 
to K = 0.4610 for C = 5. The perturbative properties of the solutions remain unchanged. 
The bump in the running coupling gets smaller but does not disappear even for ( = 5. It 
has already been stated above that the dressing functions would be independent of C in a 
complete treatment of the gluon equation. As all our solutions are very similar even on a 
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Differentprojectors 

— C=l = 
C=2 ^ 

— C=3 ; 

- - C=3.9 J 

- - C=5 : 



1 - 




-4 -'1 -? -1 n 1 ? 

10 10 10 10 10 10 10 

x=p^[GeV ^] 

Figure 3.5: Shown are the results for the gluon dressing function, the ghost dressing function 
and the running coupling, c.f. Fig. 13.41 for different projectors. 




quantitative level we conclude that transversality is lost only to a moderate extent. This 
is a somewhat surprising result in such a simple truncation scheme as the one at hand. 

The Brown-Pennington projector, = 4, is an exceptional case as can be seen from 
eqs. ()3.47p . ()3.48|) . Here the K-dependence of the second term cancels and only one solution, 
K = 1, can be found {c.f. ref. jH2j). We found no numerical solutions for this case. However, 
within the limit of numerical accuracy, solutions for ( slightly different from 4 can be found 
leading to a value for k slightly different from 1/2. E.g. in Fig. I3.5l the case C = 3.9 leading 
to K = 0.5038 is depicted. 

In Fig. 13.61 we discuss two technical issues. First, we demonstrate what happens if we 
choose two different renormalisation points by specifying two different values for a{s) = 
(7^/47r. Starting from eqs. ()3.18j) we can easily see that a change from the renormalisation 
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Figure 3.6: Here we display two technical issues. The upper panel shows the ghost and 
gluon dressing function as well as the running coupling for two different renormalisation points 
/i^ = s = 0.9 GeV^ and v"^ = t = 900 GeV^. The independence of the running coupling on 
the renormalisation point is clearly demonstrated. The lower panel shows the variation of the 
dressing functions with the parameters a and b from the construction in eq. (I3.23j) . Recall 
S = -9/44. 



point s = /i to the new value t is performed by 



Z,{t,L) 

Zix,t) = Zix,s)^^. (3.60) 



As the renormahsation constants and are independent of momentum, this results 
in the mere multiplication of the dressing functions by a constant number. Our numerical 



Chapter 3. Propagators of Landau gauge Yang-Mills theory 



51 




3 4 5 1 2 
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Figure 3.7: Solutions of the Dyson-Schwinger equations compared to recent lattice results for 
two colours [48, 49||. 

results obey exactly this behaviour, as can be seen in the upper right diagram of Fig. IH.fil 
The running coupling, however, is independent of the renormalisation point, c.f. eq. ()3.2H) . 
as is clearly demonstrated in the upper right diagram of Fig. 13.61 

In the lower panel of Fig. I3.(jl we show what happens when we vary the parameters a and 
h in our construction of the three-gluon vertex ()3.23|) . Clearly, the qualitative behaviour 
of the curves does not change. In particular the ultraviolet behaviour of the solutions 
is independent of the values of a and h in accordance with our analysis in subsection 
13.3.11 However, we recognise the rising bump in the running coupling when we lower the 
absolute values of a and h from a = 6 = 35 to a = 6 = 25. This corresponds to the 
equations becoming more and more unstable and no solutions are found for even lower 
values of a and h. The extreme case a = 6 = 0, where the gluon loop becomes leading 
in the infrared, has been investigated in reference jSSl with the same negative result: no 
solutions have been found. 

Finally, we compare our results to recent SU(2) lattice calculations^ jlHlllH]- As has 
already been stated above the ghost and gluon dressing functions from Dyson-Schwinger 
equations are independent of the numbers of colours at least to our level of truncation. 
The only caveat in comparing our results with the lattice ones is the adjustment of the 
momentum scale. We used the lattice result Cisu{2){x) = 0.68 at x = 10 GeV^ as input, 

^Note that the same set of lattice data have been analysed differently in refs. |65[I107| . A comparison 
of our results with the differently analysed data can be found in ref . [H3 ■ 
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which leads to a shghtly different scale than the one used in the case of SU(3). The two 
graphs in figure ()H.7|1 show that the main qualitative features, the infrared suppression of 
the gluon dressing function and the divergence of the ghost dressing function are common 
properties of both, the lattice solutions and the one from Dyson-Schwinger equations. 
Even the power k ^ 0.595 of the gluon dressing function from the DSEs is very close to 
the one that can be extracted from the lattice fit to be k 0.5. The main difference 
between the two approaches is in the medium energy region around one GeV, where 
the Dyson-Schwinger solutions suffer from the missing two-loop contributions that are 
certainly present in lattice Monte-Carlo simulations^. The combined evidence of the two 
methods points strongly towards an infrared vanishing or finite gluon propagator and an 
infrared singular ghost propagator in Landau gauge. 

3.5 Summary 

In this chapter we have presented approximate non-perturbative solutions for the gluon 
and the ghost propagators as well as the running coupling in Landau gauge. We obtained 
these solutions for the Dyson-Schwinger equations in a truncation scheme working with a 
bare ghost-gluon vertex and an ansatz for the three-gluon vertex such that we reproduce 
the correct one-loop scaling of the ghost and gluon dressing functions. We attempt to 
obtain two-loop scaling by the inclusion of the diagrams that involve four-gluon vertices 
in future work. 

An important improvement to previous treatments has been the explicit numerical 
calculation of all angular integrals thus overcoming the angular approximations that have 
been made so far. We could show that for a given projector only one out of two analytical 
solutions in the infrared can be connected to a numerical solution for finite momenta. 
For a transversal projector, ^ = 1, we found that the gluon propagator is only weakly 
infrared vanishing, Dgiuonip^) oc (p^)^""^, k = 0.595..., and the ghost propagator is 
highly infrared singular, Dghost{p^) oc (p^)^''^^. This is in accordance with the Kugo- 
Ojima confinement criterion and Zwanziger's horizon condition. The running coupling 
possesses an infrared fixed point with the value a(0) ~ 2.97 (or, for a general number 
of colours, a(0) ^ 8.92/N^). 

Despite the simplicity of the truncation our solutions agree remarkably well with recent 

^Despite this shortcome the DSE-solutions are surprisingly good also in this momentum regime. Of 
particular importance is the reproduction of the bump in the gluon propagator at p « 0.8 GeV. In lattice 
calculations center vortices have been identified as promising candidates for field configurations which are 
responsible for confinement [SS]. Recent investigations on the lattice suggest |108| that the bump of the 
gluon propagator is induced by such field configurations. 
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lattice calculations performed for two colours. Due to the finite lattice volume the lattice 
results cannot, of course, be extended into the far infrared. In this respect our results are 
complementary to the lattice ones: We do obtain the infrared behaviour analytically. On 
the other hand, lattice calculations include, at least in principle, all non-perturbative effects 
whereas we had to rely on truncations. E.g. the deviations for the gluon renormalisation 
functions at intermediate momenta depicted in Fig. IH. 7! might be due to the neglect of the 
four-gluon vertex function in our calculations. 



Chapter 4 



Towards general gauges in the DSEs 
of Yang— Mills theory 

In the last chapter we have exploited some properties of Landau gauge that simplify the 
Dyson-Schwinger equations of the propagators of Yang-Mills theory considerably. The 
fact that the ghost-gluon vertex is not renormalised in Landau gauge has turned out to be 
the key for the formulation of our truncation scheme and leads to a useful definition of the 
non-perturbative running coupling. In other gauges matters are much more complicated 
and no one has been able to solve the DSEs for general gauges yet. As a first step towards 
such a solution we present an exploratory study in this chapter. We investigate what 
happens, when we extend the simple truncation scheme employing bare vertices to the 
whole family of general gauges which are given by the gauge parameters a and A of the 
general Lagrangian, eq. ()2.11|) . 

Away from the Landau gauge limit no direct connection between the Kugo-Ojima 
confinement criterion and the infrared behaviour of the ghost dressing function can be seen. 
This opens the possibility that other degrees of freedom like the longitudinal gluon may 
take over the infrared dominant role of the ghost dressing function in Landau gauge. As a 
matter of fact, infrared dominance of longitudinal gluons is seen if stochastic quantisation 
is used instead of Faddeev-Popov quantisation 'IS'. 

This chapter is organised as follows: In the next section we project the ghost and gluon 
Dyson-Schwinger equations for general gauges onto the respective dressing functions. Then 
we perform an infrared analysis for the bare vertex truncation of these DSEs. We show 
that in general ghost- ant ighost symmetric gauges the genuine two-loop terms (generalised 
squint and sunset diagram) in the gluon and the ghost DSEs become important in the 
infrared. This is different from Landau gauge at least in the framework of the bare vertex 
truncation. In general ghost- ant ighost symmetric gauges the infrared behaviour of the two- 
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Figure 4.1: The coupled gluon and ghost Dyson-Schwinger equations from a BRS and Anti- 
BRS symmetric Lagrangian. Each equation contains one-loop diagrams, a tadpole contribution 
and a sunset and a squint diagram. 

loop terms exclude power-behaved solutions for the gluon and ghost propagators when bare 
vertices are employed. In section il^^ we will provide numerical solutions for the DSEs in the 
Landau gauge limit of the ghost-antighost symmetric case of the Lagrangian and recover 
the solutions found in the last chapter from a different direction in the two dimensional 
gauge parameter space |lU9j . 

4.1 Projection of the gluon equation 

In Fig. 14.1! we show again the gluon and ghost Dyson-Schwinger equations for general 
gauges in diagrammatical notation, c.f. section 12.41 The ghost and gluon equations 
are remarkably similar, both having tadpole and non-perturbative one-loop and two-loop 
contributions. 

In order to sort the various contributions of the gluon equation to the inverse of the 
gluon propagator on the left hand side we project the equation on its longitudinal and 
transverse parts. It is well known that for linear covariant gauges, i.e. a = 0, the 
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1 7 _L Mghost I T/ff'we , -trglue . ^rglue . -triad i T/sun i T/stm 

I T/.s'un I T/S'un , T/sgwint , Trsquint Tr squint ^rsquint ^rsquint 



T.L 




rsun I wsti'/i I TT/-sf/ttnit , ^^rsquint , tt rsquint . tt rsqiiint . tt Tsquint 




Figure 4.2: Various contributions from the respective diagrams in the transversal (T) and 
longitudinal (L) gluon equation and the equation for the ghost dressing function. 
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longitudinal part of the gluon propagator remains undressed ^Oj- However, away from 
linear covariant gauges this is not the case as can be seen from the Slavnov- Taylor identity 
derived in P^. We then have three dressing functions, Z,L,G, in the general case and 
the propagators are given by 

DM = P„„(p)lr + 

V-^)^ + ALO>y-^. (4.1) 

Da(j.) = (4.2) 

The transversal and longitudinal gluon dressing functions Z{p'^) and L(p^) can be ex- 
tracted by contracting the gluon equation with the transversal and longitudinal projector 
respectively. The results are given graphically in Fig. 14. 2[ where we also specify our no- 
tation for the different contributions being analysed in the next section. Contributions in 
the transversal part of the gluon equation are denoted by the symbol V , contributions in 
the longitudinal part by W and the ones in the ghost equation by U . The subscripts T 
and L indicate the respective parts of the gluon propagator running around in the loops of 
the diagrams and abbreviations for the names of the diagrams are used. For example the 
symbol denotes a contribution from the sunset diagram to the longitudinal gluon 

equation with two longitudinal and one transverse part of the gluon propagator running 
in the loop. To isolate the dressing functions the left hand sides of the equations have 
already been divided by factors of Sp^ and p^, respectively. 



4.2 Infrared analysis with bare vertices for arbitrary 
gauge parameters 

We are now ready to determine the behaviour of the two-point functions at small momenta 
p^, in a way very similar to our analysis in the last chapter for Landau gauge. In section 
13.2.11 we discussed the central observation in Landau gauge, that there is no qualitative 
difference of the solutions found with bare vertices or with vertices dressed by the use 
of Slavnov- Taylor identities. This has been shown recently for a range of possible vertex 
dressings in a truncation scheme without any angular approximations ^TOJ. The reason 
for this somewhat surprising result has been attributed to the non-renormalisation of the 
ghost-gluon vertex in Landau gauge, i.e. Zi = 1. It seems as if the violation of gauge 
invariance using a bare vertex is not that severe in Landau gauge such that the resulting 
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equations still provide meaningful results. In the following we will explore to what extent 
such a simple truncation idea is applicable in other gauges. 
Our power law ansatz for the dressing functions is 

G{x) = Bx^, Z{x) = Ax% L{x) = Cx^, (4.3) 

where x = has been used. Together with the expressions for the bare vertices given 
in appendix IA.3I we plug the power laws into the ghost and the gluon equation. The 
formulae for the various integrals are given in appendix IC.2I The straightforward but 
tedious algebra is done with the help of the algebraic manipulation program FORM [OE] . In 
section IH. 3. 21 we have argued that the renormalisation constants Z3 and Z3 can be dropped 
in the infrared analysis. Furthermore the tadpoles just give constant contributions to the 
respective propagators which vanish in the process of renormalisation. Thus we can safely 
omit them in the present investigation. 

For general gauges, a ^ and A 7^ 0, we obtain the following structure: 

B-'x-" 
A-^x-" 



{cxy^x~p 



Here the primed quantities f/', V' and W are momentum independent functions of the 
exponents /?, a and p. The corresponding momentum dependent quantities ?7, V and W 
have been introduced in Fig. 14.21 The pattern of the equation is such that each primed 
factor on the right hand side is accompanied by the squared momentum x to the power of 
the dressing function content of the respective diagram. In appendix IC.3I we demonstrate 
how such a pattern emerges for the example of the sunset diagram in the ghost equation, 
(jjyun_ ^Q^g ^^^^ ^j^g contributions (P^OlF and [W) 

are zero and therefore missing 



J-T^o'n/M'^'J"*"* -L T-Scr+P/'T/M*'?"*"* _|_ „'^(7+2pc<ri\squint 
-rx yv )rpj,rpj, -r J. {V Jj'TTL'^-^ \^ JTTLL 

+x'^^'\V%TlI + x'p{v')ZTl , (4.5) 

= w^iwy'"'"' + x^'^iWYr^T + x"+^(iy')Tr 

+x'"^^P{W'Y^^Y!l- (4-6) 
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in the longitudinal gluon equation ()4.fj|l as momentum conservation cannot hold with three 
longitudinal gluons in the three gluon vertex. 

For the following argument we focus on one particular contribution on each right hand 
side of the equations, keeping in mind that all other contributions have no explicit minus 
sign in the exponents: 





= + . . . , 


(4.7) 




— X yv j'j^XTT ~r • • • 5 


(4.8) 


[cxy^x-f 




(4.9) 



In general the coefficients (f/')*"", {V'Yx'ttt nonzero and explicitly given 

in appendix IC. 31 First, it is easy to see from equations (jUTj), ()4.8p and ()4.9p that neither 
f3 nor cr nor p can be negative. If one of these powers would be negative the limit x ^ 
would lead to a vanishing left hand side of the respective equation whereas the right hand 
side is singular in this limit. This is a contradiction as the power on the left hand side of 
the equation should match the leading power on the right hand side. Second, if one of the 
exponents j3, a or p would be positive, then the diverging left hand side of the respective 
equation would require a diverging term on the right hand side as well. However, we 
already saw that (3, a or p cannot be negative and therefore we have no diverging term 
on the right hand side of the equations. Thus none of the exponents (3, a and p can be 
positive. The remaining possibility is (3 = a = p = 0, but then we have logarithms on the 
right hand side of the equation which do not match the constant on the left hand side. 
Thus we arrive at the conclusion that there is no power solution for general gauges when 
bare vertices are used. 

There are two limits for the gauge parameters a and A in which the situation changes. 
The first one is a = 0, that are ordinary linear covariant gauges. Due to the corresponding 
Slavnov- Taylor identity the longitudinal part of the gluon propagator remains undressed, 
L(p^) = 1 ^W\. The longitudinal gluon equation becomes trivial as the Slavnov- Taylor 
identity forces all nontrivial contributions on the right hand side to cancel each other. 
In the ghost equation the squint as well as the sunset diagram disappear, i.e. {Jj'y^^ = 
(jjiyqumt _ g rjj^g term of the dressing loop which contains the longitudinal part of the 
gluon propagator, {U'Yl^'^^ , vanishes as well, c.f. appendix IC.3I We are left with the 
contribution U^^^^^ ~ x^^~^'^\ similar to Landau gauge. This is the reason why we escape 
the argument given below eq. ()4.9|1 . The explicit expression for the ghost equation is given 
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by 



167r2 2(/3 + a)(-l + /? + a) 

^ T{3 + p + a)T{2 - a)T{l - p)' 



(4.10) 



which can be compared with the Landau gauge expression, eq. ()3.42|) . Matching left and 
right hand sides we conclude that the Landau gauge result 

a = -2/3 (4.11) 

is valid for all linear covariant gauges when bare vertices are used. 

Compared to general gauges the gluon equation does not change in structure when 
a = 0. Therefore we still have a > 0, according to the argument given below eq. ()4.9j) . 
The identity ()4.11|) thus requires (3 < and we have a divergent ghost and a vanishing 
gluon dressing function in the infrared similar to Landau gauge. With bare vertices this 
implies that the ghost loop is the dominant contribution in the gluon equation for small 
momenta. Since the ghost loop is independent of the gauge parameter A, our Landau gauge 
results of section l!i.3.2l for the exponent k, = —(3 = (j/2 are valid for all linear covariant 
gauges. However, a word of caution is in order. In Landau gauge there are indications 
j66[ l68] that the general result ()4.11|1 does not change when the vertices are dressed. This 
has been confirmed recently for a range of possible vertex dressings jZOl ■ It is a completely 
open question whether this is the same for A 7^ 0. 

Having addressed the case of linear covariant gauges with a = we now turn to the 
other interesting limit, i.e. A = while a 7^ 0. It is easy to see that the a-dependence 
of the Lagrangian 1)2.111) can be eliminated in this case by partial integration using the 
constraint = 0. However, on the level of the DSEs with bare vertices there remain 

spurious a-dependent terms on the right hand side of the gluon equation. In the next 
section we will investigate the dependence of the Landau gauge solution on the gauge 
parameter a if these spurious terms are present. 



4.3 Solutions in Landau gauge 

To assess the influence of the spurious a-terms in Landau gauge we use the truncation 
scheme of section 13.31 Again we employ the general tensor 




(4.12) 
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to contract the Lorenz indices in the gluon equation. By varying the parameter ^ we are 
able to test additionally for spurious longitudinal terms. For A = and general values 
of the gauge parameter a the coupled set of equations for the ghost and gluon dressing 
functions are (c.f. eqs. (IT^ . (IT^ . TTM . THi^ . where a = 0) 

' - Z,-g^NJ^,^^i^Giy)Ziz), (4.13) 



G{x) J (27r)4 xy 

- Zs + g — —— G{y)G{z) 



Z{x) S J (27r)4 xy 



^N, f d^q Q{x,y,z)G{yY^-''/^-^^W{zY^-^/^-^^^ 



3 J (27r)4 xy Z{yY Z{zf ' ^ ' > 

The kernels ordered with respect to powers of z := p"^ = {k — qY have the form: 

r.1 . ^( {x-yf \ I [x + y\ 1 
K{x,y,z) = — [ ^ ) + - ' 



z^ \ 4 ) z\ 2 ) 4' 

M(x,y,z) = - : a; H H - — 

^ ^ z\ 4 2 4 a; y 2 2a; 4a; 

1 Aa;^ 2 19-C 5-C o Cz/^ 

Q{x,y,z) = — xy + ——y + -— 

z^ \8 y 8 4 8 X 



1 15 +c 17 -c 

+ -[ — X — y + C— 

z \ y 4 4 X 

19-Cx , 17-C^9Cl/ 



y 4 4 X 

c , 5-c^ , .2 c ,5, 



Va; 4?/ y 8xy 4^ ^ ^ ^ 

First we accomplish the infrared analysis. With k, = —(3 = a/2 we employ the ansatz 

Z{x) = Ax"^" , G{x) = Bx-"" (4.16) 

in the equations ()4.13p and ()4.14|) . After integration we match coefficients of equal powers 
on both sides of the equations and obtain 

1 {2 + k){1 + k) 
18 (3 - 2k) 

(4k-2)(-1 + k) 

~ (C - 1) [4 ^2 (a2 - 2 a + 1) + 8 K a (2 - a) + 3 a (a - 2)] + K (10 - 7 C) - 6 + 3 C ' 

(4.17) 
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FaddeevPopovgauges: a=0,2 Ghostantighostsymmetricgauge: a =1 




Figure 4.3: Here the graphical solution to equation ()4.17j) is shown. The thick line represents 
the left hand side of equation (I4.17p . whereas the other curves depict the right hand side for 
different values of the parameters The left figure shows results for a = and a = 2, 
whereas in the figure on the right we have a = 1. The ellipses mark the bulk of solutions 
between k = 0.5 and k = 0.6, whereas the circles show spurious solutions with k > 1. 

The values of k for different tensors V^'^^ can be read off Fig. 14.31 The curve given by 
the fully drawn line represents the term on the left hand side of eq. ()4.17|1 . whereas the 
other lines depict the right hand side for several values of the parameter (. The solutions 
between k = 0.5 and k, = 0.6 remain nearly unchanged when a is varied, whereas the 
solutions for k > 1 are only present when a = 0,2. In the last chapter we found that only 
the solutions k, ~ 0.5 connect to numerical results for finite momenta. It is satisfactory to 
observe that by varying the gauge parameter a such solutions can be identified as spurious 
already on the level of the infrared analysis. 

We now explore the impact of the spurious a term on the behaviour of the solutions for 
all momenta x. We solved the equations ()4.13|) and ()4.14|) numerically using the technique 
described in appendix IC.41 Compared to the usual Landau gauge, a = 0, we obtain the 
greatest deviations for the ghost-antighost symmetric case, a = 1. The results can be seen 
in Fig. ()4.4|1 . As the dependence of the kernel of the ghost loop on a vanishes in the case of 
the transverse projector, C = I5 this solution is identical with the one displayed in Fig. 13.41 
in subsection 13.4.21 For the other cases the power k changes from 0.5953 for = 1 to 
0.5020 for ( = 3.9 in accordance with the infrared analysis. The ultraviolet properties of 
the solutions are slightly disturbed compared to the case a = 0, 2. An analytical ultraviolet 
analysis similar to the one in subsection 13.3.11 reveals that the a-term in the ghost loop 
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Figure 4.4: Here we display the gluon dressing function, Z{x), the ghost dressing function, 
G{x) and the running coupling a{x) in the truncation scheme of chapter El for the gauge 
parameters a = 1 and A = and different tensors V^''^ employed in the gluon equation. 

induces a spurious dependence of the anomalous dimensions on the parameter (: 

-26- (C- l)a(2-a) 
^ ~ 44 + (C - 1) a (2 - a) 

' = 44 + (C-^M2-a) f^-^^* 

For general a only the transverse projector removes the alpha-term in the ghost equation 
and leads to the correct one loop scaling of the equations, that is 5 = —9/44 for the ghost 
and 7 = —13/22 for the gluon dressing function for an arbitrary number of colours and 
zero flavours. 
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4.4 Summary 

We have studied the infrared behaviour of the ghost and gluon propagators in covariant 
ghost-antighost symmetric gauges. We derived the corresponding Dyson-Schwinger equa- 
tions for these propagators including the ones of hnear covariant gauges as the hmit of a 
vanishing gauge parameter. Note that ghost-antighost symmetric gauges are particularly 
interesting as they allow an interpretation of the antighost field being the antiparticle of 
the ghost which includes also the possibihty of ghost-antighost condensates. Due to the 
emergence of a four-ghost interaction term in the Lagrangian for general values of the 
gauge parameters the Dyson-Schwinger equation of the ghost propagator displays a rich 
structure very similar to the one of the gluon equation. On the other hand, in the gluon 
equation we obtain the same structure as in linear covariant gauges apart from the fact 
that the gluon propagator acquires a nontrivial longitudinal part which appears in turn 
in all diagrams. The gluon and ghost equations depend therefore on three independent 
dressing functions, one for the ghost, one for the transversal part of the gluon propagator 
and one for the longitudinal one. 

We then employed a truncation scheme for the Dyson-Schwinger equations that uses 
bare vertices in place of the dressed ones. The success of this particular truncation scheme 
in Landau gauge has been attributed to the non-renormalisation of the ghost-gluon vertex, 
i.e. Zi = 1. We addressed the infrared behaviour of the ghost and gluon propagators for 
general gauges by employing power law ansdtze for the respective dressing functions. We 
then have been able to evaluate the infrared behaviour of the gluon and ghost equations 
analytically. 

For all linear covariant gauges we find a similar result as compared to the one in Landau 
gauge: An infrared suppressed gluon propagator and an infrared enhanced ghost. Whereas 
in Landau gauge there are indications that this generic result is not changed when the 
vertices are dressed [70], it remains an open question whether this is the case in linear 
covariant gauges in general. Away from linear covariant gauges, i.e. in the general case 
a and A 7^ 0, we do not find power solutions for the dressing functions. Again, this 
might be altered significantly by appropriate vertex dressings. Nevertheless, it remains to 
be emphasised that therefore also the occurrence of a ghost-antighost vacuum condensate 
is excluded in this specific truncation scheme within this class of gauges. 

A special case among all gauges considered here is Landau gauge. In the limit A = 
the general Lagrangian ()2.11|) becomes independent of the second gauge parameter a, thus 
Landau gauge is also a special case of ghost-antighost symmetric gauges. Although the 
Lagrangian of the theory is independent of the gauge parameter a, our simple truncation 
scheme breaks this invariance and spurious a-dependent terms arise in the ghost loop of the 
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gluon Dyson-Schwinger equation. Examining the case a = 1 we showed that the influence 
of these spurious terms is very smalL We determined solutions for the ghost and gluon 
dressing functions both analytically in the infrared and numerically for finite momenta and 
found solutions close to the ones of chapter El We thus recovered the results of Landau 
gauge from a different direction in the two dimensional space of gauge parameters. 



Chapter 5 



Landau gauge Yang-Mills theory on 
a four-torus 

There are a number of central aims connected to the investigation of the Dyson-Schwinger 
equations on a four-torus. The first one is purely technical: This allows to study finite 
volume effects also in the Dyson-Schwinger approach. Monte-Carlo simulations on a lattice 
necessarily have to be done in a finite volume. Therefore in the latter kind of approach 
infrared properties are only accessible by extrapolations to an infinite volume. However, 
the available data are gained on volumes which differ at best by one order of magnitude 
due to limitations in computer time. We will see in the present Dyson-Schwinger approach 
that available volumes cover several orders of magnitude. And more importantly, one can 
compare to the results obtained in an infinite volume for several truncation schemes. 

The second issue is the connection of analytical solutions obtained in the infrared limit 
of the (truncated) Dyson-Schwinger equations to the numerical solution obtained for finite 
momenta up to the ultraviolet. In chapterlHlwe have seen that not every analytical solution 
in the infrared connects to such a numerical solution. By imposing different infrared 
boundary conditions on our numerical equations we could check, whether a given infrared 
solution is connected to a numerical solution for finite momenta and which is not. In 
this chapter we will verify this procedure by an even stronger selection process: Due to 
the finite volume there is no input in the infrared and our numerical solutions on a four- 
torus choose their infrared behaviour by themselves. We will see, that the same infrared 
solutions are chosen than in the continuum. 

The third issue is that once we know how to treat Dyson-Schwinger equations on 
a torus several interesting possibilities for further investigations open up. Choosing an 
asymmetric four-torus might allow the introduction of a non-vanishing temperature 
in a relatively simple way. Furthermore there is the possibility of topological obstructions 
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Figure 5.1: Diagrammatic representation of the truncated gluon and ghost Dyson-Schwinger 
equations according to the truncation schemes studied in this chapter. 

on a compact manifold. It is well known e.g. that a four-torus allows for a non-vanishing 
Pontryagin index Moreover one could think of choosing twisted boundary conditions 

[TT2irmj . 

This chapter is organised as follows: In the first section we summarise the truncation 
schemes of jHll HHj and [HZl in more detail. They serve as simple test cases to check the 
feasibility of Dyson-Schwinger calculations on a compact manifold. In the second sec- 
tion we set up the Dyson-Schwinger equations on a four-torus and discuss their numerical 
treatment. In the last section of this chapter we display our numerical results for three 
different truncation schemes: the two summarised in the next section and the new trun- 
cation scheme introduced already in chapter El We discuss the finite volume effects of our 
solutions [SSI perform a preliminary infinite volume extrapolation. 



5.1 A summary of two truncation schemes employing 
angular approximations 

One central aim of the present investigation is to answer the question, how the ghost and 
gluon propagators are affected by finite volume effects on a torus. We will therefore com- 
pare solutions of the continuum DSEs with results obtained from the respective equations 
on the torus. To be able to distinguish genuine finite volume effects from effects arising 
only in certain truncation schemes we will study not only the scheme defined in chapter 
01 but also the two previous truncation schemes of refs. [Ml ITHI IFTTj on a torus. All three 
truncation schemes take into account the loops displayed in Fig. 15.11 However, contrary to 
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our calculation in chapter El the two previous schemes employ an angular approximation 
in the loop integrals, c.f. the discussion in section IT^ This angular approximation turned 
out to be good in the ultraviolet but less trustable for small momenta. We therefore con- 
sider the system given in chapter El to capture more of the essential physics of the coupled 
ghost and gluon system. Nevertheless when it comes to a study of finite volume effects it 
seems helpful to employ all three truncation systems. 

To proceed we will summarise the truncation schemes of refs. [^^11131111 in the next two 
subsections. The main difference between these two truncation schemes is the treatment 
of the three-point functions. Whereas in ref. jnEj approximate Slavnov-Taylor identities 
have been employed to construct an ansatz for the vertices, in ref. [HZ! bare three-point 
functions have been used. Amazingly, though, both schemes provide results with identical 
qualitative infrared behaviour: the gluon propagator vanishes in the infrared, the ghost 
propagator is highly singular there, and the strong running coupling has an infrared fixed 
point. 



5.1.1 The dressed vertex truncation including the gluon loop 

In section 13.21 we already discussed the full Slavnov-Taylor identity for the ghost-gluon 
vertex, c.f. eq. ()3.9|) . Neglecting the irreducible four-ghost term this identity and the 
corresponding one for the three-gluon vertex are solved by the vertices 



r.(p,g) = + , (5.1) 



G(g2) G{p^ 

+2 — \ {S^iuP-q - Piyqti)ip - q)p + cyclic permut. , (5.2) 



with 



J J ,.2> GIF) / G(,=) ^ G(p^ 



(5.3) 



Here k"^ = {p — qY denotes the gluon momentum in the ghost-gluon vertex F^. The 
three-gluon vertex T^^p is completely symmetric in the three momenta. 

The Lorenz indices in the gluon equation are contracted with the Brown-Pennington 
projector, eq. ()3.6|) . In order to obtain the correct scaling behaviour of the gluon loop in 
the ultraviolet the substitution 

/y/ 2A(-2-6<5) 
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for the vertex renormalisation constant Z\ has been introduced. Here the squared momen- 
tum g2 denotes the momentum integrated over in the gluon loop. Recall the anomalous 
dimension of the ghost, b = —9/44, for an arbitrary number of colours and zero flavours. 
Furthermore in ref. ^66j the so called modified angular approximations 

Z{z),G{z)c^Z{x),G{x) ^ Z{y),G{y) if x<y, 

Z{z),G{z) — ^ Z{x),G{x) if x>y, (5.5) 

have been employed. The approximation in the first equation is designed to eliminate a 
spurious term stemming from the vertex construction (c./. the discussion in subsection 
I3.2.1|) . Since the presence of this spurious term causes a fatal inconsistency in the ghost 
equation [HZj, this modified angular approximation is a central ingredient in this truncation 
scheme. 

Collecting all this together one arrives at the equations 

^ - Zs- N^j£y^ I / ^ ^(^) ^(^) + / A ^ G{y) Z{y) \ , (5.6) 




G(x) 

V. 

^ = ^3 + ^7^1 /a ^^^^^^(y)-^^G(y)(--^^) 
Z{x) 6 (27r)^ ' ^ 



( X 

G''(x) + G(x) Giy) 
xyz 

i^i^ + ^1 G-(y) ) . (5.7) 

XZ xyz 

X 

for the ghost and gluon dressing function For ease of notation we have used the 

abbreviations x := p^, y := and z := {p — qY for the squared momenta. Furthermore 
an 0(4)-invariant momentum cutoff L = h? has been introduced. The integral kernel 
N{x, y, z) is given by 

N{x, y, z) = {4:Z^ + 32z^y + 2z^x - 26z^xy - Idz'^x^ - 72z^y^ + 32zy^ 

-38zyx'^ + 8x^z - 2Qzy'^x + 2y^x + x^ + Ay'^ - Iby'^x'^ + 8x^y) , (5.8) 

and the angle 9 is defined hj z = {q — pY = x + y — 2^Jxycos 9. The equations ()5.6p and 
()5.7|) will be implemented on a four-torus in section 15.21 

In the continuum one is now able to carry out the angular integrals analytically and 
solve the equations along the lines described in ref. [011 IHH HHj. In the infrared the 
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solutions Z{x) and G{x) behave power-like, 

Z{x) = Ax^^ 



G{x) = Bx 



—ft 



(5.9) 



with uniquely related coefficients A and B. In this truncation scheme one obtains k ^ 0.92. 
For three colours the running coupling approaches the fixed point «(0) ~ 9.5 in the infrared 



The dressed vertex truncation scheme of Hauck, Smekal and Alkofer has been the 
first one to include both the gluon and the ghost dressing function. Compared to the 
old truncation of Mandelstam j^S], which completely neglects the effects of ghosts, this 
has been a major improvement. Whereas the Mandelstam equation is solved by a gluon 
dressing function which diverges in the infrared, the inclusion of ghosts leads to the qual- 
itatively different picture of an infrared vanishing gluon dressing function and a diverging 
ghost. This new picture, contradicting the old idea of infrared slavery, has been corrobo- 
rated since in other Dyson- Schwinger studies [HZl IHHl ED] as well as lattice calculations 



One of the original strengths of the dressed vertex truncation scheme, namely the con- 
struction of the vertices as solutions of approximate Slavnov- Taylor identities, is considered 
as a weakness in the meantime. The so constructed vertex ansdtze have been shown to be 
not in accordance with perturbation theory jHSl IHH] and, much worse, lead to inconsistent 
equations when no angular approximations are employed. This is the central reason why 
we did not use the vertices ()5.H1 and ()5.2|1 in our truncation scheme of chapter El 

5.1.2 The bare vertex 'ghost-loop only' truncation 

In ref. jHZl also a bare ghost-gluon vertex has been used. In section we have discussed 
at length why the bare vertex is capable of providing reliable results even in the infrared 
region of momentum, where one would expect effects from non-perturbative vertex dressing 
to be most pronounced. In the gluon loop the authors of [HZj use a bare three-gluon vertex 
without modifying the vertex renormalisation constant Zi. As this construction neither 
restores a correct perturbative limit of the equations (c./. our discussion in section 
nor changes the infrared behaviour of the solutions, the authors themselves omit the gluon 
loop in the main part of their investigation. 

Substituting the tree-level ghost-gluon vertex for the dressed one and neglecting the 
gluon loop the coupled system of equations ()5.7|1 reads 



msi ED 113 cniEi 113 cni 



1 



I 
I 



d'^q K{x,y,z) 



G{y)Z{z) , 



(5.10) 



G{x) 



(27r)4 xy 
d^q M{x, y, z) 
(27r)4 IS^ 



1 




G{y)G{z) . 



(5.11) 
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Again we used the abbreviations x := p^, y := g^, z := {q — p)^, s := fi"^ and L := A^. 
The kernels K and M are aheady given in eqs. ()H.2(j|l . The ghost equation ()5.1()j) and 
the ghost loop ()5.11|) are identical to the ones of our truncation in section 13.31 provided 
the projection parameter ( is set to the special case C = 4 (c./. eqs. fl3.7p . ()3.24|) . ()3.25|) ). 
As the ghost loop is the leading part of the gluon equation, both schemes share the 
same infrared behaviour for ^ = 4. Indeed, we have found the infrared solution k = 1 
in subsection 13.3.21 in accordance with the infrared analysis in jH2], where no angular 
approximation has been employed. However, we were not able to find a numerical solution 
for finite momenta connecting to k = 1 in the infrared as has been detailed in section 

Thus, similar to the truncation scheme summarised in the last subsection, one only gets 
numerical solutions for finite momenta once an angular approximation has been employed. 
The authors of ^| use the so called ymax- approximation 

Z{z),G{z) — > Z{max{x,y)),G{msix{x,y)). (5.12) 

Upon angular integration the eqs. ()5.10|) . ()5.11|) are then simplified to 

In the infrared the solutions of eqs. ()5.13I5.14|) behave power-like, c.f. eq. fl5.9p . with 
K ^ 0.77. The running coupling approaches the fix point a(0) = 11.47 in the infrared [HT] . 

The 'ghost-loop only' truncation scheme of Atkinson and Bloch successfully proved 
that bare vertices result in the same qualitative behaviour of the ghost and gluon dressing 
functions in the infrared as the dressed construction of the last subsection. This surprising 
result served as basis for recent analytical investigations in the infrared, where either a bare 
ghost-gluon vertex jHS 1^2] or a bare ghost-gluon vertex with multiplicative corrections 
[Tin ini] has been employed. The weak point of the 'ghost-loop only' truncation is its 
failure in the ultraviolet, where the absence of the gluon loop leads to a contradiction with 
perturbation theory. 

5.2 Finite volume effects on a four-torus T4 

From a technical point of view using a four-torus as the underlying manifold or choosing 
(anti-)periodic boundary conditions on a hypercube is identical. The first question to 
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answer is therefore: are periodic or antiperiodic boundary conditions adequate for the 
fields appearing in the coupled system of ghost and gluon Dyson- Schwinger equations. 
This question is easily answered for the gluon field. Due to the bosonic nature of gluonic 
excitations we have to use periodic boundary conditions for the gluons and subsequently 
for the gluon dressing function Z . Although this is not obvious, the same is true for the 
ghosts. The reason for this 'strange' behaviour of Grassmann fields is to be searched in 
the definition of the Faddeev-Popov determinant, which due to its introduction as gauge 
fixing device shares the symmetry properties of the gluon field |115t lllb] . An easy way to 
see this is to check the BRS-transformation rule of the quark field \E', which already has 
been given in eq. ()2.29j) . 



Here t°- is a generator of S\J {Nc) and s is the BRS-operator generating the transformation. 
Since s is continuous it does not change the symmetry properties of the fields. We thus 
have antiperiodic behaviour on both sides of the equation due to the fermionic nature 
of the quark field This, however, is only consistent if the ghost field c has periodic 
boundary conditions. 

With periodic boundary conditions on a hypercube with length / in every direction the 
four-dimensional momentum integrals of our DSEs have to be substituted by a sum over 
four indices, 



The quantities of interest, the gluon and ghost renormalisation functions Z{p^) and G{p'^) 
depend only on the 0(4) invariant squared momenta as all directions on the torus are 
treated on an equal footing. This suggests to relabel the points on the momentum grid 
not according to a Cartesian but a hyperspherical coordinate system. 



where the index j numbers the hyperspheres = const, i.e. the circles drawn in the sketch 
in Fig. 15.21 The index m, which numbers the grid points on each hypersphere respectively, 
will be suppressed in the following. 

In the integrals to be discretised there appear three momenta, the external momentum, 
labelled p, the loop momentum q and for the second propagator in the loop k = q—p. We 
will use the following notation: 



(5.15) 




(5.16) 




(5.17) 



X : = p^ 

y ■= (f 
z := k'^ = {q — pY 



with 



with 



with 



Xi G hypersphere i , 
Uj ^ hypersphere j , 
Zn G hypersphere n . 



(5.18) 
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Figure 5.2: Sketch of the momentum grid dual to the four-torus and the summation over 
complete hyperspheres indicated by fully drawn circles. The hyperspheres depicted by dashed 
lines are not complete due to the numerical ultraviolet cutoff in every direction of the grid. 

On the hypercubic momentum grid dual to the four-torus the momentum k = q — p is 
located on the grid for every pair of grid momenta p and q as can be seen from elementary 
vector operations. This is no longer true once we introduce a momentum cutoff corre- 
sponding to a finite extent of the momentum grid. We will detail below, how we treat the 
dressing functions at those momenta z which are larger than the cutoff. 

For the moment let us think a little on the appropriate way to introduce such a cutoff 
on a grid. Recall that an 0(4) invariant cutoff A has been introduced into all of the 
continuum DSEs eqs. dSHI) and (jOTl|l . (jEIII). In order to compare 

the continuum results with the ones on the torus a corresponding regularisation of the 
sums over grid momenta is required. A first idea might be to simply cut the sums in 
each direction according to ji,j2,j3,j4 = — A^, . . . , 0, . . . A^. Such a method, however, 
breaks 0(4) invariance and introduces numerical errors, as will be shown in section IB. 21 
(Note that in lattice Monte-Carlo simulations the analysis of the resulting data for 0(4) 
invariance necessitates special kinds of cuts through the lattice, see e.g. [12].) As can be 
seen from Fig. 15.21 an 0(4) invariant cutoff of the sums necessitates to neglect the 'edges': 
The sum extends only over the fully drawn hyperspheres, and we omit the summation 
over the dashed ones. 

The main difference between the equations on the torus and the ones in the continuum 
is the effective treatment in the infrared. The zero modes at j = 0, which are not present 
in the continuum, are neglected in all calculations presented here^. The finite volume in 
coordinate space leads to a finite value of squared momentum for the first hypersphere 
j = 1. Thus one has not to worry about possible infrared singularities. Furthermore 

^An estimate of their possible contribution is given in appendix ID. 21 
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the complicated matching procedure of the infrared integrals to the remaining integration 
region in the continuum is not necessary on a torus. On the other hand one might worry 
that the infrared part of the loop integrals are crucial to obtain solutions at all p7|. In 
appendix ID. II we show that this is not the case. We reproduce the continuum results in 
the 'ghost loop only' truncation scheme on a very coarse radial momentum grid without a 
single point in the infrared. This gives us first confidence that Dyson- Schwinger equations 
might be feasible on a torus. 

Considering all this, one anticipates already at this level some deviations in the infrared 
between the solutions obtained in these different ways. After writing down the DSEs on 
a torus in the next two subsections we discuss our numerical results. They demonstrate 
that using a torus as infrared cutoff works surprisingly well. 



5.2.1 Angular approximated DSEs on a torus 

Substituting the replacement rule / -^1^ — > j4^j into the eqs. ()5.1()j) . (jS.llj) and using 
the angular approximation (j5.12|l the DSEs in bare vertex ghost-loop only truncation read 
on a torus 



= Z,{s, L) - g'N^y^ Yl G{y,)Z{m^x{x., y,)) , (5.19) 

3 

= Us, L) + E ^^"^^ '"^ G(y,)G(max(x., y,)) . (5.20) 



All arguments Xi = (2277//)^ and yj = (j27r/Z)^ of the dressing functions G and Z are on 
the momentum grid. However, note that Zn = (q—p)^ might be larger than the ultraviolet 
cutoff even if Xj and yj are not. Nevertheless the kernels in eqs. ()5.19p . ()5.20j) can be 
calculated straightforwardly according to the expressions in eqs. flH.2(')j) . 

The corresponding equations for the dressed vertex truncation on a torus can be derived 
analogously in a straightforward manner from eqs. ()5.(j|l and ()5.7j] . However, as the 
expressions are quite lengthy and bring nothing new we do not give their explicit form. 



Chapter 5. Landau gauge Yang-Mills theory on a four-torus 



75 



5.2.2 The novel truncation scheme on a torus 

With the replacement / F discussed above, the eqs. f)3.24|) . ()3.25|) for the 

novel truncation scheme of chapter El read on the torus: 

^ = Us,L)-g^N^^^j:^^^^^^^^ (5.21) 

The kernels M and Q are given in eqs. ()3.26|) . ()3.32|) . As already stated y/z might be 
larger than the ultraviolet cutoff \/L even if ^/x and ^Jy are not. Nevertheless the kernels 
can be evaluated straightforwardly. However, if z with L < z < 4L is the argument of a 
dressing function one has the choice of two different methods. One way is to approximate 
Z{z) and G{z) by Z{L) and G{L). Another more elaborate treatment consists of matching 
the corresponding perturbative ultraviolet tail to the function under consideration. We 
have applied both methods and found only very small quantitative differences. 



5.3 Renormalisation and results 
5.3.1 The renormahsation scheme 

In order to obtain comparable results for the Dyson-Schwinger equations in the continuum 
and on the torus we have to impose the same renormalisation conditions. This can be done 
in two ways: First, one can use the solutions of the continuum equations for fixed cutoff A 
and fixed renormalisation scale ^ to calculate the corresponding renormalisation constants 
Z3(/i^, A^) and Z^^fi"^ , A"^) . These can subsequently be used in the equations on the torus. 
Second, one can subtract the torus equations at the squared momenta sq and sz and 
trade the two renormalisation constants for the values of the dressing functions at these 
momenta, namely Z(sz) and G{sg) {c.f. subsection 13.4.11) . These values are then taken 
from the continuum solution. If sz and sq are sufficiently far in the ultraviolet region of 
momentum, where finite volume effects play a minor role, the two procedures lead to the 
same results within the limits of numerical accuracy. 

In the bare vertex truncation scheme we have chosen the renormalisation condition 
= Glfi"^) = 1 for the continuum equations. For the renormalisation point /i^ we took 
the same value as for the ultraviolet cutoff: A^ = /x^ = 0.2. Of course, this choice is by 
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no means special and one is completely free to choose the renormalisation point wherever 
one likes. We subtracted the continuum ghost equation at zero momentum and the gluon 
equation at the renormalisation point /i^. This is convenient as we are then able to use the 
condition Z{sz = /i^) = 1 directly as input in the calculation. The second input is provided 
by the coefficient A of the leading order infrared expansion of the gluon dressing function, 
Zir{x) = Ax^*^. The condition G{fi'^) = 1 corresponds to A = 357.33 in our calculation. 
(Note that the coefficient A is uniquely related to the corresponding coefficient B of the 
ghost dressing function in the infrared, c.f. the discussion in section 13.3.21 ) The value of 
the coupling at the renormalisation point is taken to be = g'^/'in = 0.97. Again this 

number is completely arbitrary provided one stays in the interval a{0) > a{^i^) > a(A^). 
In the torus equations we furthermore used the values Z^lfi"^ = 0.2, = 0.2) = 0.9591 
and Z^du,"^ = 0.2, A^ = 0.2) = 1.1034, which have been determined from the continuum 
solution. 

A physical momentum scale can be determined only in truncation schemes which pro- 
vide the correct perturbative running of the coupling in the ultraviolet. With the missing 
gluon-loop this is not the case in the 'ghost-loop only' truncation scheme. We therefore 
have to stick to an internal momentum scale without physical units in this case. The sit- 
uation is different, however, in the dressed vertex truncation scheme. Here we have fixed 
the momentum scale by calculating the running coupling for the colour group SU(3) and 
using the experimental value a{x) = 0.118 at x = M| = (91.187 GeV)^ to fix a physical 
scale. 

In the numerical treatment of the continuum equations of the dressed vertex truncation 
scheme the ghost equation is subtracted at zero momentum and the gluon equation at 
the arbitrary finite momentum s = 1.048 GeV^. We solved both equations similar to 
the method described in ref. [HI], especially we introduced also the auxiliary functions 
F{x) and R{x) as defined in ref. jHE]- As input values serve the infrared expansion of 
R{x), R{x) = + . . ., and the value R{s) = 0.8 at the gluon subtraction point. For the 
calculations on the torus we determined the values Z^^jj? = M|, A^ = 1.255 GeV^) = 1.266 
and 2'3(yU^ = M|,A^ = 1.255 GeV^) = 0.966 from the continuum solution. 

For the truncation scheme of chapter El we have described in section 13.4.11 how we 
obtain the continuum results. The corresponding torus results have been calculated with 
torus equations both subtracted at the renormalisation point /i^ = 1.9 GeV^. The input 
value for the gluon dressing function at this momentum is Z{fi^) = 0.83. Requiring 
G(/i^) = 1/ \fZ{Jj^ then fixes the normalisation for G{x). For the value of the coupling 
at the renormalisation point we chose a(/i^) = 0.97 similar to the two other truncation 
schemes. 
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Barevertices,ghost-looponly 




x=p ^[internalunits] 

Figure 5.3: Shown are the ghost dressing function, G{x), the gluon dressing function, Z{x), 
and the running coupling, a{x), in the bare-vertex ghost-loop only truncation for different 
momentum grid spacings corresponding to different finite volumes of a torus. The fully drawn 
lines labelled continuum represent the respective results for continuous momenta. 

5.3.2 Numerical solutions 

Our results for the ghost dressing function, the gluon dressing function and the running 
coupling in the bare-vertex 'ghost-loop only' and the dressed vertex truncation can be seen 
in Figs. 15.31 and 15.41 In both truncation schemes we solved for three different momentum 
spacings corresponding to different volumes in coordinate space. To keep the cutoff iden- 
tical for all the spacings within each truncation scheme we have chosen three different grid 
sizes respectively. For the bare vertex truncation they are N = 17^, 31^, 51^ and for the 
dressed vertex truncation they are N = 21^, 41^, 61^. 
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Dressedvertices,ghost-andgluon-loop Dressedvertices,ghost-andgluon-loop 




Dressedvertices,ghost-andgluon-loop 

10^ 



10° 



10 ^ 10 ^ 10 ' 10° 

x=p^ [GeV^] 

Figure 5.4: The same as Fig. I5.3l for the dressed vertex truncation. 

The finite volume effects seen in both truncation schemes are very similar: Compared 
to the respective continuum solutions the ones obtained on a torus show deviations for 
the first few lattice points in the infrared. For large momenta all functions obtained on 
a torus approach the continuum ones. The biggest effect can be seen for the running 
coupling a. As aijP') is proportional to the product Z (jP')G'^ {p^) the deviations of the 
torus dressing functions from the continuum curves amplify in the infrared in a somewhat 
erratic way, such that the data points at small momenta cannot be connected by a smooth 
line. Comparing larger and smaller spacings of momentum grids one clearly sees that the 
effect is always one of the first spheres on the respective lattices and therefore moves to 
the infrared for smaller spacings. 

The most important properties of the continuum solutions can still be found in the 
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Figure 5.5: Here we compare the running coupling of the bare vertex ghost-loop only truncation 
for different treatments of the lattice cut-off, c.f. section The torus curve on the left panel 
is calculated on a grid where the edges have been cut, whereas the coupling on the right panel 
is obtained on a full hypercubic grid. The ellipses mark the region where the latter treatment 
is inferior. 



torus solutions despite some deviations in the infrared. Going from larger to smaller 
spacings a power-like behaviour of the dressing functions in the infrared with the correct 
exponents can still be inferred. For the truncation scheme with dressed vertices the gluon 
dressing function on the torus has the same shape and the same height of the bump in the 
bending region of the curve. This is a first central result of our investigation: Employing 
a torus as infrared regularisation is possible. 

Before we move on to discuss the results in our novel truncation scheme we show what 
happens, if the four dimensional momentum grid is not cut at the edges according to the 
prescription discussed in section 15.21 In Fig. 15.51 we compare the running coupling of the 
bare vertex 'ghost-loop only' truncation for different treatments of the lattice cut off. The 
torus curve in the left panel is calculated on a grid with edges cut, whereas the coupling 
on the right panel is obtained on a full hypercubic grid. Clearly one observes sizeable 
deviations in the ultraviolet behaviour of the latter curve compared to the continuum 
result. As explained in section EIH this is a result of breaking rotational symmetry on 
the lattice. On the inner, full hyperspheres we have the discrete rotational symmetry 
group Z(4) for very small discrete steps, whereas these steps become larger and larger for 
those hyperspheres which contain points on the edges of the lattice. Cutting these edges 
obviously improves the ultraviolet behaviour of our solutions. 

The numerical results for our novel truncation scheme of chapter El with a transverse 
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projector, ( = 1, are shown in Fig. 15. 61 At the time those results were found they have 
been the first numerical solutions at all in a truncation scheme without any angular ap- 
proximations In the meantime, as we have seen in chapter El solutions for continuous 
momenta are also available [H5j. Both calculations are compared in Fig. 15.61 For the 
presented torus solutions on three different volumes lattice sizes of = 13^,43^, 71'^ have 
been used^. As mentioned in the last subsection we did check cutoff effects by extrapolat- 
ing the propagator functions at z > A^ with a logarithmic tail with the correct anomalous 
dimensions. The results as compared to the ones obtained by simply setting Z{z) = Z{h?) 
and G{z) = G(A^) for all z > A^ change by less than one per mille. 

In section 13.3.21 we raised the question, which one out of two analytical solutions of 
the continuum DSEs in the infrared is connected to a full numerical solution for finite 
momenta. We have answered this question partially in section IHTil where we showed that a 
numerical solution exists which connects to the infrared solution k, = 0.595. However, in a 
numerical calculation in the continuum the analytical infrared solutions are used as input 
in the infrared and reproduced self-consistently. Thus one might argue that we showed 
the existence but not the uniqueness of the numerical solution with k = 0.595. On a torus 
we have no input in the infrared. Therefore a torus solution favours automatically the 
physical solution in the infrared, as we believe the infinite volume limit to be smooth. The 
results shown in Fig. 15.61 are clearly in agreement with the power behaviour k = 0.595 
already found in the continuum and disfavours the second solution, k ^ 1.3. The infrared 
critical exponent as calculated in refs. [0211101, = 0.595, thus has been verified. 

The gluon dressing function in Fig. 15.61 is remarkably stable against changes of the 
volume and approaches more and more the expected power solution for small momenta, 
although this process seems to be very slow. For the ghost dressing function one observes 
deviations of the first points in the infrared: An extraction of the correct infrared critical 
exponent from the numerical solution for the ghost function is hardly possible. Only some 
points come close to the analytical value of the continuum before the curve starts bending 
down again for small momenta. For the extracted value of the running coupling in the 
infrared this leads to a distinct mismatch to what one is to expect on the basis of analytical 
results. 

At first sight the fact that the power solution for the ghost dressing function could 
be not reproduced numerically to a reasonable precision may seem disappointing. Never- 
theless these numerical results themselves show that the ghost dressing function is highly 

^Note that the scale for the torus results shown in Fig. 15.61 is slightly different than in reference |83) . 
There an extrapolation to the mass of the Z-boson has been used, whereas here we simply take the scale 
which has been determined for the continuum result. 
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QQ x^=1.4x10'' 
— Continuum 
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Figure 5.6: Here we present the ghost dressing function, G{x), the gluon dressing function, 
Z{x), and the running coupling, a{x), in the truncation scheme of chapter E] for different 
volumes using a transverse projector, i.e. ( = 1. 

singular in the infrared. This reflects the long-range correlation of ghosts in Landau gauge. 
Therefore one should expect the ghost dressing function to be the one affected most by a 
finite volume. On the contrary the gluon dressing function vanishes in the infrared and 
consequently it is much less affected by a finite volume. We expect the ghost dressing 
function together with the running coupling to approach more and more the correct power 
solution in the infrared as lattice spacings are decreased and lattice sizes are increased. 

Furthermore, we add a remark on the transversality of the gluon propagator, c.f. the 
respective discussion in section 13.11 Our numerical results on a torus for different values 
of the parameter ( can be seen in Fig. 15.71 Although the solutions show the expected 
dependence on the form of the projector this dependence is not too drastic and in general 
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100 r 




Figure 5.7: The same as Fig. I5.6l for different projectors. The corresponding solutions in the 
continuum are given in Fig. 13.51 

the behaviour of these different solutions is very similar. For the gluon dressing function 
one observes that the more ( grows the greater is the deviation from the pure power 
behaviour and correspondingly from the continuum solutions, which have been shown in 
subsection 13.4.21 The points at small momenta cannot be connected by a smooth line 
any more. Based on the infrared analysis in subsection 13.3.21 one might anticipate that k 
should approach the value k, = 0.5 more and more as ( grows until there is a jump to the 
solution from k = 0.5~^ to k = 1 as the Brown-Pennington limit C = 4 is reached. We do 
not observe such a qualitative jump in our solutions on a torus. The solution shown for 
C = 4 is approached smoothly when ( approaches this limit. This clearly indicates that 
the solution n = 1 might not exist at all if one removes the torus as a regulator. 

Finally we compare our results of the new truncation scheme to recent SU(2) lattice 
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Figure 5.8: Results on the torus compared to recent lattice results (IHI EHI. As the torus 
points are very close to each other on a linear momentum scale we did not resolve the torus 
curves into single points. 

calculations jlHlllH], c./. the corresponding comparison between lattice and DSE results 
in section 13.4.21 The two graphs in Fig. ()5.8|) suggest that the differences between our 
solutions in the continuum and on a torus are much smaller than the difference to the 
lattice result, at least for intermediate momenta. We explained in chapter IHl that this 
difference can be attributed to full two-loop contributions which are missed out in the 
DSE calculations but are certainly present in lattice Monte-Carlo simulations. For small 
momenta we see that our calculations on a torus employ much larger volumes than is 
possible in lattice simulations up to now. A thorough investigation of the infinite volume 
limit of the DSE solutions on a torus including a careful analysis of possible cut-off effects 
is certainly desirable and under way. 



5.4 Summary 

In this chapter we have presented numerical solutions of truncated systems of Dyson- 
Schwinger equations for the gluon and ghost propagators in Landau gauge SU(N) Yang- 
Mills theories. We have employed a four-torus, i.e. a compact space-time manifold, as 
an infrared regulator. Apart from the infrared finite volume effects encountered on such 
a manifold we found solutions on a torus which are very close to the ones obtained in 
the continuum formulation for various truncation schemes. Thus a central result of this 
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chapter is: Dyson-Schwinger equations on a torus are feasible. 

For small momenta we found the expected finite volume effects in the solutions of the 
DSEs on a torus. These deviations from the continuum results are somewhat erratic for 
the two truncation schemes employing an angular approximation in the sense that the 
values of the dressing functions on the momentum grid cannot be connected by a straight 
line. For our new truncation scheme without angular approximations this is different: 
systematic deviations from the continuum solutions occur. We find small deviations for 
the gluon dressing function and larger differences for the ghost dressing function on a torus 
compared to the continuum one. This is what one expects, since the infrared dominant 
ghost correlation is long ranged and should therefore be more affected by the finite volume 
than the gluon correlation. 

Employing a torus as infrared regulator has certain advantages compared to the con- 
tinuum formulation. In the numerical calculation of the continuum DSEs one has to match 
the analytically determined solution in the infrared carefully to the region where the dress- 
ing functions are determined numerically. On a torus the finite volume in coordinate space 
leads to a finite value of the smallest squared momentum encountered in the calculations. 
Thus one has not to worry about infrared singularities. Furthermore it is much easier to 
obtain convergence in the numerical iteration process. Historically we solved the trun- 
cation scheme introduced in chapter IHl on the torus first and only subsequently in the 
continuum using the torus solutions as start values in the iteration process. Thus the 
torus formulation is an important technical tool. 

Furthermore we used the Dyson-Schwinger equations on a torus to corroborate a result 
of chapter 13 only one out of two analytical solutions from the infrared analysis of the 
continuum equation is connected to numerical results. For the transverse projection of 
the gluon equation this is the solution corresponding to k = 0.595. Such a statement 
is stronger if it is inferred from solutions on a torus than from corresponding continuum 
results. The reason is that we have to use the analytical infrared solution as self consistent 
input in the continuum calculations to obtain numerical stability. This is not the case on 
a torus. 



Chapter 6 



The coupled system of quark, gluon 
and ghost DSEs 

In this chapter we will enlarge our focus from pure Yang-Mills theory to Landau gauge 
QCD. We will investigate the coupled system of Dyson- Schwinger equations for the ghost, 
gluon and the quark propagators. In the quark DSE we will study the mechanism by which 
physical quark masses are generated even though the bare quark masses in the Lagrangian 
are zero. This is a genuine effect of Strong QCD. It is well known that for vanishing bare 
masses the renormalised masses remain zero at each order of perturbation theory. 

A thorough study of the infrared phenomena in the quark sector of QCD requires 
a continuum formulation. Lattice simulations of dynamical chiral symmetry breaking 
|117| I118| I119| 112m I121j have to deal with finite volume effects and in addition have to 
extrapolate from finite to zero quark mass. It is not possible to put massless quarks 
on a finite lattice. A recently performed study p."TT?j in the overlap formalism e.g. em- 
ploys masses in the range of mo = (126 — 734) MeV. These values suggest that even the 
most elaborate extrapolation method to zero quark mass needs guidance and check from 
continuum results. 

Apart from the phenomenon of mass generation we are interested in quark confinement. 
Single quark states have non-vanishing colour charge and are therefore not contained in 
the physical part of the state space of QCD. In subsection 12.3.21 we argued for a positive 
(semi-) definite metric in this physical subspace, whereas the remaining state space of QCD 
contains negative norm states as well. Consequently, negative norm contributions to the 
quark propagator are theoretical evidence for quark confinement. 

The quark propagator is an important ingredient for many phenomenological mod- 
els (see ini Hm I122j and references therein). Thus, the quark DSEs have been studied 
extensively. Various ansdtze for the gluon interaction in the quark equation have been 
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Figure 6.1: A diagrammatical representation of the quark Dyson-Schwinger equation. 

explored. The resulting quark propagators have been used in bound state calculations 
based on the Bethe-Salpeter equations for mesons (see e.g. j2Hl 11^31 1124| I125| 1126] ) or 



Faddeev equations for baryons [EZl CHDl CnH IIH3 HHHl CSSl CnZl CSHI • 



One of the central aims of this thesis is to provide a solution for the quark propagator 
which incorporates the effects from the ghost and gluon DSEs directly and not via model 
assumptions. 

This chapter is organised as follows: In the next section we will construct suitable 
ansdtze for the quark-gluon vertex such that two important properties of the full quark DSE 
are reproduced: the independence of the generated quark mass from the renormalisation 
point and the asymptotic matching of the DSE-solutions to the results of perturbation 
theory. Fortunately, the corresponding DSEs for the fermions of QED are well studied 
(a short overview is given e.g. in |139j ). We will dwell on these results and construct 
non-Abelian generalisations of Abelian vertices, which have the desired properties. 

In the second section we present solutions for the quenched system of quark, ghost 
and gluon DSEs, i.e. we neglect the quark-loop in the gluon equation. A corresponding 
approximation is frequently used in lattice simulations. Comparing our solutions with 
recent lattice results |119j we find very good agreement for the quark propagator. 

We then proceed to the unquenched case and incorporate the quark-loop into our 
truncation scheme for the ghost and gluon DSE from chapter El We present solutions 
for the full coupled system of DSEs for the quark, ghost and gluon propagators in the 
last section of this chapter. Compared to the quenched case we will find only moderate 
differences for the number of light flavours Nj < 3 jl4Uj . 

6.1 The quark Dyson-Schwinger equation 

We outlined the derivation of the Dyson-Schwinger equation for the quark propagator S{p) 
at the end of section 12.41 The renormalised equation in fiat Euclidean space-time is given 
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by 

S-\p) = Z^S^\p) + ^^Z,pCf j d\^,S{q)T,{q,k)D,,{k), (6.1) 

with the momentum routing k = q — p. The factor Cp = {N^ — l)/2Nc in front of the 
integral stems from the colour trace of the loop. The symbol r,y{q,k) denotes the full 
quark-gluon vertex. A diagrammatical representation of the equation is given in Fig. 16.11 
Suppressing colour indices the quark and gluon propagators in Landau gauge are given 

by 

S(v) = ^ = A-Uv')^*^-^^ (6 2) 

- (6-3) 

where the quark mass function M is defined as M(p^) = B{p'^)/A{p'^). For obvious reasons 
the dressing function A{q'^) is frequently called 'vector self energy' and the dressing function 
B{q^) 'scalar self energy'. Furthermore the inverse of the vector self energy, 1/A{q^), is 
denoted as 'quark wave function renormahsation'. The bare quark propagator Soip"^) 
contains the unrenormalised quark mass mo(A^) which depends on the cutoff A of the 
theory. The bare mass is related to the renormalised mass m/j(/i^) via the renormahsation 
constant Z^'- 

moiA^)=Zmifi^A^)mRifi^). (6.5) 

Here /x^ is the squared renormahsation point. 

The renormalised and unrenormalised vector self energy, A and Aq, are related by 

A,\p^A') = Z2(/x^ A^) A-\p',fi'). (6.6) 

In Landau gauge the quantum corrections to the vector self energy are finite. Correspond- 
ingly Z2{fi'^, A^) stays finite when the cutoff is sent to infinity and we have < Z2{fi'^, A^) < 
1. Furthermore the ghost-gluon vertex is not ultraviolet divergent in Landau gauge and 
we can choose Zi = 1, c.f. section The Slavnov- Taylor identity for the quark-gluon 
vertex renormahsation factor Zip thus simplifies, 

„ Zo Zi , , 

Z,p = ^ = f. (6.7) 

Previous studies of the quark equation in the so called Abelian approximation (an 
overview is given in ref . |^ ) as well as the recent investigation in ref . |141j assume implicit 
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cancellations between the full quark-gluon vertex, the dressed gluon propagator and the 
integral over the kernel of the DSE. Furthermore in the tensor structure of the quark-gluon 
vertex only a term proportional to 7^ is employed. 

In this thesis we do not have to rely on implicit cancellations since we calculated explicit 
solutions for the dressed gluon and ghost propagators, c.f. section 13.4.21 We will also 
construct explicit non-perturbative ansdtze for the quark-gluon vertex including different 
tensor structures than 7^. The advantages of such a treatment are obvious: every building 
block of the equation is explicitly given and well under control. We are able to assess the 
expedience of different vertex constructions^. Furthermore the inclusion of several tensor 
structures in the vertex is supposed to be important in bound state calculations of scalar 
meson masses |143| I144j . Finally we hope providing an explicit construction of all parts 
of the quark equation leads in turn to more thorough statements on quark confinement. 

In the following we assume an effective non-Abelian quark-gluon vertex of the form 

r,(g, k) = V:'''ip, g, k) W-^'^'ip, q, fc), (6.8) 

with p and q denoting the quark momenta and k the gluon momentum. The non-Abelian 
factor W^"'^'^^ multiplies an Abelian part V^^'^\ which carries the tensor structure of the 
vertex. This ansatz is motivated by the aim to respect gauge invariance as much as possible 
on the present level of truncation. The Slavnov- Taylor identity (STI) for the quark-gluon 
vertex is given by ^ 

(_,) G-\e) k, T^{q,p) = S'\p) H{q,p) - H{q,p) S'\q), (6.9) 

with the ghost dressing function G{k'^) and the ghost-quark scattering kernel H{q,p). At 
present the non-perturbative behaviour of the ghost-quark scattering kernel is unknown. 
Therefore we cannot solve the STI explicitly. However, comparing the structure of eq. ()6.9|) 
with the corresponding Ward identity of QED, 

i-t) k, r«^^(g,p) = S-\p) - S~\q), (6.10) 

we are able to infer some information: Whereas the ghost fields of QED decouple from 
the theory and consequently do not show up in the Ward identity, there is an explicit 
factor of G~~^{k'^) on the left hand side of eq. ()6.9|) . We therefore suspect the quark-gluon 
vertex of QCD to contain an additional factor of G{k'^) compared to the fermion-photon 
vertex of QED. Some additional ghost dependent structure seems necessary to account for 

^Results for the quark-gluon vertex on the lattice have been obtained in ref. |142) . However, at present 
the error bars from such simulations are too large to use the lattice results as guideline in the construction 
of our ansdtze. 
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the ghost-quark scattering kernel on the right hand side of eq. For simphcity we 

assume the whole ghost dependence of the vertex to be contained in a non-Abelian factor 
multiplying an Abelian tensor structure^. 

The Abelian part of the vertex, V^'^^\ can be adopted from QED. The Ward identity 
(I6.10|l has long been solved |147|I148] . Furthermore transverse parts of the fermion-photon 
vertex have been fixed by Curtis and Pennington^ to satisfy multiplicative renormalisability 
in the Abelian fermion DSE for all linear covariant gauges |155| I156j . 

The non-Abelian factor VT^'^^^' is chosen such that the complete quark equation fulfils 
two conditions: 

(i) The quark mass function M{p^) should be independent of the renormalisation point 

(ii) The anomalous dimension 7^ of the mass function known from resummed pertur- 
bation theory should be recovered in the ultraviolet. 

In the course of this section we will prove the vertex ansatz 

-2d-d/5 



_ Giz,s)Zsis,L) 
W-^'^\p,q,k) = G\z,s) Z,{s,L) , (6.11) 



[Ziz,s)Zsis,L)) 



A{x,s) + A{y,s) 



2 



A(x, s) — A(y, s) , , ... . ,B(x, s) — B(y, s) . 

+ 7^ + i){p + q)u + ^ P + q)u 

2[x-y) x-y 

+ ^^""^'^'^^^y^'^ [(x - y)7. -{i>-mp+ q)A X 

" + ^ (6.12) 



(a;-?/)2 + (M2(x) + M2(?/))2 ' 

with the new parameter d to satisfy the conditions (i) and (ii). To ease notation we used 
the abbreviations x = p'^, y = q"^ and z = {p — qY for the squared momenta, also s = /i^ 
for the squared renormalisation point and L = h? for the squared cutoff of the theory. 
The anomalous dimension 5 of the ghost propagator is 5 = —9 Nd (44 Nc — 9> Nj) at one 
loop order for Nc colours and Nf flavours. The Abelian part of the vertex is given by the 
Curtis-Pennington (CP) vertex r'^^{p,q,k). 

^For quenched calculations and in the context of angular approximated DSEs a similar strategy has 
already been adopted in refs. | 145[ I146| . 

•^This so called Curtis-Pennington vertex has been used extensively in the study of the fermion DSE 
in QED [Tl^^llWniTrniT^IT^IT^ . 
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From a systematic point of view the newly introduced parameter d in the non-Abehan 
part of the vertex is completely arbitrary. Our numerical results, however, will indicate 
that values around the somewhat natural choice ci = match best with lattice simulations, 
c.f. subsection 16.2.31 

If we would not care about (Abelian) gauge invariance we could also employ the much 
simpler vertex 

-2d-d/S 



, _ G{z,s)Z,{s,L) 
W-^''\p,q,k) = G'{z,s)Zs{s,A) 



{Z{z,s)Z,{s,L)Y 

V:''\p,q,k) = Z,{s,L)^, (6.13) 

where we have taken the bare Abelian vertex, 7,^, multiplied with an extra factor of Z2. In 
Landau gauge this construction also satisfies the conditions (i) and (ii), as will be shown 
in the next two subsections^. 

6.1.1 Multiplicative renormalisability of the quark equation 

To proceed we substitute the vertex ansatz fl6.13|) into the quark equation ()6.1|) . By taking 
the Dirac trace once with and once without multiplying the equation with ^ we project 
out the mass function M{x) and the vector self energy A{x). We arrive at 

M{x)A{x,s) = Z,{s,L) mo{L) + ^^^ J d'ql^ ^ ^^"^^'^^^^^^ Z,{s, L) A~\y,s)x 

/ _ \ -2d-d/5 

(G{z,s)Z3{s,L)) ] 

,,,d my) , (6.14) 

{Z{z,s) Z3{s,L)) J 

-2d-d/5 



Giz,s) Z3is,L)j / + y ^ (s-y)2 



{Ziz,s)Z,is,L)f V 2 2^ 

(6.15) 

where we have used the definition of the running coupling a in Landau gauge 

a{x) = — Z{x, s)G^(x, s) = a{s) Z{x, s) G^(x, s) , (6.16) 
47r 



■*In the numerical treatment we will additionally employ a vertex where the transverse part of the CP- 
vertex is left out, i.e. a generalised Ball-Chiu (BC) vertex. In Landau gauge such a vertex also satisfies 
the conditions (i) and (ii). 
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c.f. eq. (ISI2H). 

The behaviour of eqs. ()(i.l4|l . ()(i.l5|l under renormahsation can be explored by changing 
the renormahsation point s = /i^ to a new point t = v'^. We first note that the factor 
stemming from the non-Abehan part of the quark-gluon vertex is not affected by such a 
change: 

{G{z,s)Z,{s,L)) ^ (G{z,t)Zs{t,L)j 



{Z{z,s)Z,{s,L)f {Z{z,t)Z,{t,L)f 
This can be seen easily with the help of the relations 

Go{x,L) = G{x,s)Zs{s,L), (6.17) 
Zoix,L) = Z{x,s)Zsis,L), (6.18) 

between the unrenormalised and renormalised ghost and gluon dressing function, c.f. 
eq. ()3.18p . Furthermore the running coupling a{z) is independent of the renormahsation 
point, c.f. the discussion in subsection 13.2.21 From eq. ()6.6|) we infer 

Z2{t, L) A-^{x, t) = ^2(5, L) A-\x, s) . (6.19) 

With the renormahsation condition A{t,t) = 1 we have 

Z2it,L)=Z2is,L)A-\t,s), (6.20) 

and subsequently 

A{x,t) = A{x,s) A-\t,s) . (6.21) 

Substituting eqs. ()6.2()j) and ()6.2H) into the Dyson-Schwinger equations ()6.14|1 we find the 
mass function M{x) to be independent of the renormahsation point, i.e. condition (i) is 
satisfied. Note that without the extra factor of Z2 in the Abelian part of the vertex ()6.13|) 
we would violate this condition. 

Before we examine the case of the more sophisticated Curtis-Pennington type vertex 
()6.12|) . we mention two important points. First, according to perturbation theory we have 
A{x — > 00, s) —>■ 1 and ^2(5, L) 1 for large renormahsation points s. However, this is 
just a special case of the general relation 

A{x ^ oo,s) ^ Z2{s,L) (6.22) 

which can be inferred from eqs. ()6.20j) and ()6.21|) . In Fig. 16.21 we sketch the vector self 
energy renormalised at two different points s and t, with s in the perturbative and t in 
the non-perturbative region of momentum. We will see a similar picture when we present 
our numerical solutions in section 16.2.31 



92 



6.1. The quark Dyson-Schwinger equation 



A 



Z(t,L) 



1 




A(x,s) 



A(x,t) 



s 



X 



Figure 6.2: Sketch of a finite renormalisation from a perturbative point s to a non-perturbative 
point t for the vector self energy A. 

Second, the alert reader might worry about the appearance of the ghost and gluon 
renormahsation factors Z^[L) and Z^{L) in the interaction kernel of the quark equation. 
Certainly the renormalised functions M{x) and A{x) should not depend on the cutoff of 
the integral. The balance of cutoff dependent quantities in the equation is controlled by 
various factors of Z2{L) and Zm{L). We have to take care not to disturb this balance 
by the vertex ansatz. Thus the non-Abelian part of our quark-gluon vertex contains such 
powers of Z-^^L) and Z^{L) that the cutoff dependence of these quantities cancel. This 
can be easily checked using the scaling behaviour 



of the renormalisation factors for L ^ oo and the relation 7 + 25 + 1 = 0, c.f. eq. ()H.H4|1 . 

Along the same lines as for the bare vertex construction we prove condition (i) for the 
Curtis- Pennington type vertex. Plugging eqs. ()6.11|) . ()6.12|) into the quark equation ()6.1|) 
and projecting onto M{x) and A{x) we arrive at 




(6.23) 



(^) {g{z,s)Z^{s,L)^ 



2d-d/S 



M{x)A{x, s) 




z{y + M\y)) {Z{z,s)Z3{s,L)Y 




(6.24) 
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A(x,s) = Z2(s,L)H -^:r^ (Tq — } 



-2d-d/5 



xz{y + M\y)) (Z{z, s)Z;{s, L)f 

x + y jx - y y\ A{x, s) + A{y, s) 
^ 2 2z 

AA 



[X 



+ y) + ABM{y)^ (^-^ + (x + y) - ^^^^ 



+ l{A{x, s) - A{y, s))M{y)n{x, y) (^LjI _ / ^ 



(6.25) 



Here we have used the abbreviations 

A{x,s)-A{y,s) 



AA 
AB 

n{x,y) 



x-y 

B{x,s) - B{y,s) 
x-y 

X + y 



(x-|/)2 + (M2(x) +M2(y))2- 

With the help of the relations ()6.20|) and ()6.2H) we find the quark equations ()6.24|) . ()6.25|) 
to be consistently renormalised if and only if the scalar self energy behaves like 

B{x, t) = B{x, s) A-^t, s) (6.26) 

and the mass function 

M{x) = B{x, s)/A{x, s) (6.27) 

is independent of the renormalisation point. 

Summing up the results of this subsection we have shown that both vertex construc- 
tions, eqs. ()6.12|1 and ()6.13|1 . lead to a renormalisation point independent mass function 
as is required by condition (i) of section 16.11 The same is true for a Ball-Chiu type 
vertex, which is the Curtis-Pennington construction ()6.12j) without the transverse term 
proportional to Q{x, y). Note that in different gauges than Landau gauge only the Curtis- 
Pennington construction would satisfy condition (i), similar to QED |155t 1156'] . 

6.1.2 Ultraviolet analysis of the quark equation 

In this subsection we will show that the ansdtze fl6.12j) and ()6.13p for the quark-gluon 
vertex both lead to the correct perturbative limit of the quark mass function M{x). We 
first examine the case of the bare vertex construction, eq. ()6.14|1 . 
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The ghost and gluon dressing functions G and Z are slowly varying for large momenta 
according to their perturbative limit given in eqs. ()3.27j) . ()3.28j) . For loop momenta y 
larger than the external momentum x we are therefore justified to employ the angular 
approximation G{z)^Z{z) G{y), Z{y), c.f. subsection 13.3.11 Furthermore there is a 
region xq < y < x where the approximation G{z), Z{z) — >• G{x), Z{x) is adequate. We are 
then able to carry out the angular integrals in eq. ()f).14j) with the help of the formulae given 
in appendix IC.ll If we additionally take the external momentum x to be large enough, 
then all masses in the denominators become negligible since the integral is dominated by 
loop momenta y ^ x. We then obtain 



M{x)A{x,s) 



Z2{s,L)mo{L) 
Z2{s, L) a{x) 



+ 



TT 



X 



X 

j dyZ2{s,L) A-\y,s) 



G{z, s)Zsis,L) 



-2d-d/5 



xo 



[Z{z,s)Z^{s,L)Y 

-2d-d/5 



■M{y) 



Z,{s,L) f a{y) yG{z,s)Z,{s,L) 
dy— Z2{s,L) A {y,s) 



TT 



y 



{Z{z,s)Zsis,L)y 



-M{y) , 
(6.28) 



where the integral from y = to y = Xq has already been neglected. 

For large momenta y > xq the wave function renormalisation A~^ and the renor- 
malisation factor Z2 cancel each other according to eq. ()6.22j) . Furthermore we use the 
perturbative limit of the ghost and gluon dressing functions {c.f. eqs. ()3.27p . ()3.28|) ) 

r /rN iS 

G{s) 



G{z) 
Z(z) 



Z(s) 



Lulog (^-1+1 
log ) + 1 



(6.29) 



with u; = /5oa(s)/(47r) = (llA^'c — 2A^/)a;(s)/(127r). If we additionally substitute the scaling 
behaviour of the renormalisation constants Z3 and Z3, eqs. ()6.23j) . and exploit the relation 
7 + 2(5 + 1 = 0, eq. ()3.34|) . we arrive at 



M{x) = moiL) + / dyM{y) + - [ dy^M{y) . 

TT X J n J y 

Xo X 



(6.30) 



This well known equation describes the ultraviolet behaviour of the quark mass func- 
tion. The classification of its solutions has been clarified by Miransky |157|I158] . employing 
the perturbative form of the running coupling, 

-1 



a{y) = a{s) 



u log 



(6.31) 
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In the chiral limit, mo(L) = 0, we obtain the so called regular asymptotic form, 



M(x) 



X (iln(x/A|cD)) 



(6.32) 



Here {'^'$) denotes the chiral condensate which is discussed in more detail in subsection 
16.2.11 In the case of non- vanishing bare quark mass, niQ^L) ^ 0, the equation ()6.3()j) is 
solved by the irregular asymptotic form. 



M{x) 

In this case we furthermore find 



M{s) 



u; log ( - 

s 



12 



mo(L) = M(s) 



tulog ( - 



(6.33) 

(6.34) 
(6.35) 



in accordance with resummed perturbation theory. 

We thus have shown that the bare vertex construction ()6.13|1 admits a solution for the 
mass function M(x), which has the correct perturbative behaviour for large momenta. A 
similar analysis is possible for the DSE with the Curtis- Pennington type vertex, eq. ()6.24|) . 
As the vector self energy goes to a constant in the limit of large momenta, eq. ()6.22|) . all 
terms proportional to A{x) — Aiy) are suppressed in this limit. Furthermore, according to 
the perturbative expression ()6.33j) the Ai?-term contributes at most subleading logarithmic 



corrections in eq. ()6.24|) . The first term in the brackets reduces to the bare vertex form 
because A{x, s) ~ A{y, s) for large momenta x, y. Thus we obtain the same ultraviolet 
limit from eq. ()6.24|1 than for the bare vertex construction. This is certainly also the case 
if a Ball-Chiu type vertex is employed. 



6.2 The quark propagator in quenched QCD 

In the course of this section we will compare our results for three different vertex types, 
which share the non-Abelian part proposed in eq. ()6.1H) but differ in their Abelian parts. 
We will employ the bare vertex, eqs. ()6.13j) . and the Curtis- Pennington (CP) type vertex, 
eqs. ()6.12j) . Furthermore we use a Ball-Chiu (BC) type construction, which employs only 
the first three terms of the CP-vertex. In Landau gauge all these vertex ansdtze satisfy the 
conditions (i) and (ii) formulated in section 16.11 In order to compare the different vertex 
types on a quantitative level we will calculate the pion decay constant and the chiral 
condensate from the respective solutions for the quark mass function. 
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6.2.1 Pion decay constant, chiral condensate and quark masses 

The pion decay constant is calculated from the Bethe-Salpeter equation (BSE), which 
describes the pion as bound state of quark and antiquark |125j . Apart from the dressed 
quark propagator the BSE involves couplings between quarks and gluons. On the level 
of the quark DSE we have substituted the full quark-gluon vertex by an vertex ansatz. 
However, at present it is only known for certain cases how such a vertex ansatz in the 
quark DSE translates to the corresponding quark-gluon coupling in the Bethe-Salpeter 
equation |159|ll(i()] . No method is known up to know to derive the corresponding BSE for 
dressed quark-gluon vertices as the BC- or CP-vertex constructions. 
We thus have to rely on the approximation 

f2 [a M{y)A~\y) [ ^ y dM{y) \ 

= -TTi / dyy , ,^o. M{y) - -—7— , (6.36) 



47r2 7 {y^M\y)f\ 2 dy J 

which incorporates only the effects of the leading pion Bethe-Salpeter amplitude in the 
chiral limit |122j . From a comparison of the relative size of the amplitudes in model 
calculations jl24| I143j one concludes that the approximation ()6.36|1 should lead to an 
underestimation of /,r by roughly ten percent. 

The renormalisation point independent chiial condensate, i^"^), can be extracted from 
the ultraviolet behaviour of the quark mass function in the chiral limit (c./. eq. ()6.32j) ): 



Mix) . • (6.37) 



Recall X = p^ and L = A^, where A is the cutoff of our theory not to be confused with 
the scale Aqcd, which is to be taken from a fit to the running coupling, c.f. eqs. fl3.57p . 
(ITHHll . 

The renormalisation point dependent chiral condensate (^^E')/^ can be calculated via 



/74 
^^,,(g2,s), (6.38) 

where the trace is over Dirac indices, Sch is the quark propagator in the chiral limit and the 
squared renormalisation point is denoted by s = /i^. To one-loop order both expressions 
for the condensate are connected by 

= Hf^V^lco)) (^^) , (6.39) 
with the anomalous dimension 7^ of the quark mass function. 
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For the calculation of the chiral condensate we first have to determine the mass renor- 
malisation constant Zm{s,L). Recall the formal structure of the mass equation ()(i.l4j] . 
which is given as 

M{x)A{x, s) = Z^is, L) Zm{s, L) m^(s) + Z^^s, L) liM^x, s) , (6.40) 

where IIm{x, s) represents the dressing loop. In order to extract Zm{s, L) from this equa- 
tion we have to clarify the meaning of ^^(s) which is related to the unrenormalised mass 
by 

mo{L) = Zm{s,L)mii{s), (6.41) 

c.f. eq. ()6.5|1 . Evaluating eq. ()6.4()j) at the perturbative momentum x = s the matter 
seems clear. We achieve consistency with eqs. ()6.33|) and ()6.35|) . if 

mnis) = M{s), (6.42) 
Zm{s,L) = \u;\og(-] +1] , (6.43) 



s 



which is indeed the correct perturbative scaling of the renormalisation constant Zm 

Certainly one could implicitly define the finite parts of Z^ such that the relation ()6.42|) 
holds in general for all renormalisation points s. Then the parameter mji in the renor- 
malised QCD-Lagrangian would already know about dynamical symmetry breaking. How- 
ever, as the mass parameters of QCD are supposed to be generated in the electroweak 
sector of the standard model one could equally well argue that it is more systematic to 
exclude the effect of mass generation by strong interaction from mn- 

In our numerical calculations we will choose s to be sufficiently large, therefore eq. ()6.42|1 
is valid anyway. Then Z^ is determined by 



_ M(x) A{x, s) - Z2{s, L) IlAf (x, s) 
^ " Z,{s,L)M{s) 
1 YiM{s,s) 



Z2{s,L) Mis) 



(6.44) 



For the last equation we have set x = s and have used the renormalisation condition 
A{s,s) = 1. 

In the numerical calculations we have to specify the masses mR{s) as input. Choosing 
a perturbative renormalisation point s allows one to evolve the masses mR{s) to a different 
scale t by 



'■QCDJ 
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6.2. The quark propagator in quenched QCD 



For t = {2 GeV)^ typical values for the masses of the light quarks are given by the Particle 
Data Group [TTl^ : 

^(m„ + md){2 GeV) ^ 4.5MeV, m,(2 GeV) ^ lOOMeV. (6.46) 
We will use similar masses in our calculations. 



6.2.2 Renormalisation scheme and numerical method 

In the quark equation we employ a MOM regularisation scheme^ similar to the one used 
in the ghost and gluon equations in chapter El The formal structure of the quark equation 
is given by 

A{x,s) = Z2{s,L) + Z2{s,L)Ua{x,s), (6.47) 
M{x)A{x,s) = Z2{s,L)ZUs,L)mji{s) + Z2{s,L)UM{x,s). (6.48) 

We eliminate Z2 from the first equation by isolating it on the left hand side and subtracting 
the same equation for x = s. With 

^ ^ +— ^nA(x,s), (6.49) 



Z2{s,L) A{x,s) A{x,s) 

we then have 



^ 1- . Ua{x,s) + Ua{s,s), (6.50) 



A{x, s) A{x, s 

using the renormalisation condition A{s, s) = 1. In each iteration step we determine the 
vector self energy A{x) from eq. ()6.5()|1 and subsequently Z2 from eq. ()6.49|1 ^. For the mass 
function M{x) we use 

M{x)A{x,s) = Z2{s,L)Um{x,s) (6.51) 
in the chiral limit and the subtracted equation 

M{x)A{x, s) = M{s) + Z2{s, L) Hm^x, s) - ^2(5, L) 0*1(5, s) (6.52) 

if chiral symmetry is broken explicitly, i.e. itlq 7^ 0. 

For the numerical iteration we employ a Newton method and represent the dressing 
functions A{x) and M(x) with the help of Chebychev polynomials. Furthermore, we use 
a numerical infrared cutoff e, which is taken small enough for the numerical results to be 



^In quenched QED4 this technique has already been used in |154l . 

check for the numerics is to determine Z2 at different momenta x = and search for an artificial 
momentum dependence of Z2. In our calculations we find Z2 to be independent of to an excellent 
degree. 
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independent of e. Numerical difficulties arise in the case of the Curt is- Pennington type 
vertex and even more for the Ball-Chiu construction. If the external momentum x and the 
loop momentum y are both small and close to each other then the derivative-like terms 

, , A(x) - A(y) , B(x) - Biy) 

X — y X — y 

are hard to evaluate accurately. Although the functions A{x) and B{x) are constant in 
the infrared and consequently should have derivatives close to zero one encounters large 
values for AA and AB due to numerical inaccuracies in A and B. In order to evaluate 
AA and AB much more precisely at small momenta we fit the expressions 

A{x) = — , Bix) = r^772 , (6.54) 

with the parameters ai, 02, &i and 62 to the numerically evaluated functions. The scale 
Aqcd = 0.714 GeV has been determined from our fits to the running coupling in section 
13.4.21 For X — y smaller than a suitable matching point we calculate the terms AA and 
AB from the ffis. This procedure eliminates the numerical errors in the derivative terms 
and smoothes the numerical results considerably. In the case of the Ball-Chiu type vertex 
the iteration process does not converge unless we use these ffis. 

The renormalisation condition employed in the ghost-gluon system of equations is 
G'^{s)Z{s) = 1 with a{s) = 0.2 at the squared renormalisation point s = fi^. Furthermore 
we choose the transversal tensor, ( = 1, to contract the gluon equation, c.f. chapter El The 
physical scale in the quenched calculations is taken directly from the Yang-Mills results of 
section 12321 i-G. we use the experimental value a(M|) = 0.118 of the running coupling 
at the mass of the Z-boson to fix the scale. 



6.2.3 Numerical results 

In Fig. 16.31 we give our numerical solutions for the quark mass function and the inverse 
vector self energy in the chiral limit. We compare results obtained with five different 
ansdtze for the quark-gluon vertex. For the generalised CP-vertex we investigate the 
'natural' case d = 0, the value d = 6 = —9/44, already adopted in refs. |145tll46j . and the 
value d = 0.1. Furthermore we employed the bare vertex construction and a Ball-Chiu 
type vertex. The corresponding masses at the momentum = 0, the pion decay constant 

the renormalisation point independent chiral condensate and the ffi parameters for the 
functions ()6.54p are displayed in table 16.11 

The numerical results for the mass function all have a characteristic plateau in the 
infrared and show the regular asymptotic behaviour for large momenta, c.f. eq. ()6.37j) . 
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6.2. The quark propagator in quenched QCD 





M(0) 
[MeV] 


fn 

[MeV] 


(_(i|r^))l/3 

[MeV] (calc.) 


(_(i|r^))l/3 

[MeV] (fit) 




02 


bi 


b2 


bare vertex 


177 


38.5 


162 


160 


3.05 


0.99 


0.06 


1.00 


CP d=S 


150 


50.5 


223 


225 


- 


- 


- 


- 


BC- vertex 


293 


62.6 


276 


284 


1.10 


0.99 


0.29 


0.92 


CP-vertex 


369 


78.7 


303 


300 


0.83 


0.99 


0.20 


1.00 


CP d=0.1 


464 


87.5 


334 


330 


0.79 


0.99 


0.34 


0.95 



Table 6.1: The mass M(0), the pion decay constant f^, calculated with eq. (I6.36j) . the 
renormalisation point independent chiral condensate calculated with eqs. (I6.38p and ()6.39|) . 
and the condensate obtained by fitting the expression (I6.37j) to the chiral mass function in 
the ultraviolet for all four vertex types. Recall S = —9/44 in quenched approximation. If not 
stated otherwise the parameter d in the vertex construction is taken to be = 0. For the case 
of the CP-vertex with = 5 we did not get good fits in the infrared. 

The bare vertex construction and the CP type vertex with d = S both generate masses 
much smaller than typical phenomenological values of 300 — 400 MeV. The BC- and the 
CP-type construction with d = provide good results, whereas the choice d = 0.1 leads 
to a somewhat large mass. The lattice calculations taken from ref. jll9j favour masses 
around 300 MeV with the caveat that they are obtained by an extrapolation from sizeable 
bare quark masses to the chiral limit''. The numerical solutions for the wave function 
renormalisation 1/A can be seen in the right diagram of Fig. 16.31 Whereas the ultraviolet 
asymptotic behaviour of all vertex constructions is similar we observe sizeable differences 
for small momenta. Again the bare vertex construction and the CP-vertex with d = 6 are 
clearly disfavoured by the lattice data. 

Our approximate calculation of the pion decay constant should underestimate the ex- 
perimental value = 93 MeV by about ten percent, c.f. the discussion below eq. ()6.36p . 
We thus have best results for the CP- vertex construction with d = and d = 0.1. Fur- 
thermore we obtain very good agreement between the two different methods to extract the 
chiral condensate, c.f. subsection 16.2.11 Compared to the value (— (^\E'))^/'^ = 227MeV 
from the phenomenological study summarised in ref. |162j most of our results are larger. 

Our favourite vertex for further investigations will be the CP-type construction. Al- 
though the bare vertex construction is by no means capable to reproduce phenomenological 
values of M(0) and f.^ we will keep this vertex for the sake of comparison. However, we 

^This extrapolation is intricate, as can be inferred from the data points in the ultraviolet region of 
momentum, where the expected regular asymptotic behaviour in the chiral limit is not reproduced by the 
lattice data. 



Chapter 6. The coupled system of quark, gluon and ghost DSEs 



101 



10 N I I I I I I I n 




10^ 10' Id' 10° lo' 10" 10^ 10* lO' 

x=p^[GeV^] 




' latticedata 



Q^ F I I I I I I I i 

10"' 10" 10"' 10° lo' 10" lo' lO'' lo' 

x=p^[GeV^] 



Figure 6.3: The mass function M(x) and the inverse vector self energy 1/A{x) of a chiral 
quark are shown. We compare the results for five different vertices with lattice data taken from 
ref. \m\. 

discard the BC-type vertex due to numerical problems in the infrared. 

Apart from the case d = 6 we obtain very good fits for the scalar and vector self energy, 
A{x) and B{x), for small momenta. The results for the fit parameters can be found in 
table If). 11 It is interesting to note that the exponents 02 and 62 in the fit functions of 
eq. ()6.54|) are all fitted very close to one. Such a behaviour might turn out to be crucial 
for the continuation of the quark propagator to negative p^, i.e. timelike momenta, in 
future work^. 

A technical point is illustrated in Fig. 16.41 Here we display the dressing functions 
calculated with the bare vertex and the CP-vertex construction {d = 0) for two different 
renormalisation points, s = 8100 GeV^ and t = 1.9 GeV^. Although these two points differ 
by three orders of magnitude the resulting mass functions for each vertex are indistinguish- 
able in the plot. The wave function renormalisation 1/A{x,s) is multiplied by a constant 
factor when we change from s to t, c.f. eq. ()6.2H1 . We thus find the behaviour expected 

^There are several possibilities how the quark propagator might look like for negative |163[ 1164] , 
Confinement seems to require that there are no poles on the negative p^-axis. This would prevent a 
Kallen-Lchmann spectral representation |lf)5| which is mandatory for a propagator representing a physical 
particle. In phenomenological models two different ansdtze for the quark propagator have been explored: 
an exponential form with a pole at — > —00 (see e.g. jl24[ 11251 1166j ) and a form proposed by Stingl 
|167j with complex conjugate poles in the negative half of the p^-plane |168[ 1135] , It has been shown, 
however, that the exponential form is ruled out by scattering processes |135) . 
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6.2. The quark propagator in quenched QCD 




Figure 6.4: The behaviour of M(x) and 1/A{x,s) under a change of the renormalisation 
point from /i^ = s to /i^ = t. 
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Figure 6.5: In these diagrams we show that a change of arguments, C^^z) G{z)G{y), in 
the non-Abelian part of the quark-gluon vertex is disfavoured by lattice data. 



from the discussion in subsection 16.1.11 

In Fig. 16.51 we compare the results for two different momentum assignments in the 
non-Abehan part of the vertex. We employ 

W-'''"'\x, y, z) = G\z, s) Z3(s, L) , (6.55) 

in accordance with the d = case of eq. 1)6.111) and the modified assignment 

W-'^'^^ix, y, z) = G{z, s) G{y, s) Z,{s, L). (6.56) 
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Figure 6.6: Results for three different forms of the running coupling in the quark equation: 
The running coupling calculated in chapter Eland the two fits given in eqs. ()3.57|) . (I3.58|) . 



Both expressions behave similar under finite renormahsation and become equal in the 
limit of large momenta. Therefore both expressions lead to a mass function satisfying the 
conditions (i) and (ii) formulated in section IU?T1 However, as can be seen from the plot, 
the lattice data for the wave function renormahsation 1/A clearly disfavour the second 
momentum assignment. This will be important in the next section, where we investigate 
the quark-loop in the gluon equation. 

Fig. 16.61 compares results for the bare vertex and the CP-type construction for three 
different forms of the running coupling in the interaction kernel of the quark equation. The 
two fit-functions, 'Fit A' and 'Fit B', have been given in eqs. ()3.57|) . ()3.58|) Furthermore we 
used the running coupling calculated from the quenched ghost and gluon DSEs in section 
13.4.21 Although there is the (presumably) artificial bump at ~ 0.1 GeV^ in the running 
coupling, the mass functions obtained from the DSE-result and from 'Fit A' are virtually 
indistinguishable. 'Fit B', however leads to somewhat smaller masses. This observation 
suggests that nearly all the dynamically generated mass is produced from the integration 
strength above p = 500 MeV, as is indicated by the vertical line in the plot of the running 
coupling. This is a favourable result as it would have been very unsatisfying if the artificial 
bump contributed a considerable amount to the quark mass function. 

Finally we observe the effects of explicit chiral symmetry breaking in the plots of 
Fig. 16. 71 We give results for three different quark masses, m(2 GeV) = 5 MeV, m(2 GeV) = 
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6.3. Incorporating the quark-loop in the gluon equation 
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Figure 6.7: These diagrams show our results when three different bare quark masses are 
employed. In the diagram on the right small quark masses correspond to small values for 1/A 
in the infrared. 



100 MeV and m(l GeV) = 1000 MeV. These values correspond roughly to the ones given by 
the Particle Data Group for the up/down-quark, the strange-quark and the charm-quark 
|102j . For small momenta we note again that the dressed vertex generates more mass in the 
quark equation than the bare vertex construction. This effect becomes much less dominant 
for the heavy quarks, where more and more of the infrared mass stems from explicit chiral 
symmetry breaking and not from dynamical mass generation. Furthermore in accordance 
with our analysis in subsection Ifj . 1 . 21 we observe the same ultraviolet behaviour of the mass 
function for both vertex constructions. 



6.3 Incorporating the quark-loop in the gluon equa- 
tion 

In section 13.4.21 and in the last section we have discussed results for the ghost, gluon 
and quark propagators in quenched approximation. Now we go one step further and 
investigate the unquenched equations, i.e. we include the back-reaction of the quarks 
on the ghost-gluon system. To this end we incorporate the quark-loop in the truncation 
scheme developed in chapter IHl Due to our experience with the ghost- and gluon-loop 
this might seem to be a straightforward task. However, we will encounter some additional 
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Figure 6.8: Diagrammatical representation of the quark-loop in the gluon equation. 

problems which are related to the quark-gluon vertex. 

The contribution of the quark-loop to the gluon equation is given by 

nrr' = / A TT{j,S{q)r^{q,k)S{k)}, (6.57) 

where x = p'^ is the external gluon momentum and y = and z = k'^ = [q — pY the 
squared momenta of the two quarks running in the loop. The trace is over Dirac indices. 

In section lUTTl we have proposed an effective quark-gluon vertex Ty{q, k) with Abelian 
and non- Abelian parts such that the quark equation is multiplicatively renormalisable and 
one-loop perturbation theory is recovered for large momenta. However, this construction 
is not capable to account as well for the one-loop behaviour of the unquenched gluon 
equation unless we switch the arguments of the non- Abelian part eq. (jFnT|l . to 

different momenta. In the quark equation such a change of momentum leads to results 
clearly disfavoured by the lattice data, c.f. Fig. 16.51 We therefore have to use different 
momentum assignments for the quark-loop and the quark equation. Certainly, this is 
a deficiency which has to be resolved by a more elaborate vertex construction in future 
work. The aim of the present study, however, is to present an effective construction which 
captures essential properties of the theory. 

Taking care of symmetries we propose the following ansatz for the non- Abelian part of 
the quark-gluon vertex in the quark-loop: 

G{y)Z,{L)) [G{z)Z,{L)j 
{Z{y)Z,{L)f' {Z{z)Z,{L)f' 

(6.58) 



W;:^:Lioo,{^.y.z) = G{y)G{z)Z,{L) 



Here x = is the squared gluon momentum, y = q^ and z = k"^ = {q — pY are the squared 
quark momenta, and L = A^ is the squared cutoff. The Abelian part of the vertex is 
already symmetric with respect to the quark momenta. 
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6.3. Incorporating the quark-loop in the gluon equation 



Plugging the Curtis-Pennington type vertex into the quark-loop and contracting the 
free Lorenz-indices with the tensor [c.f. eq. ()H.7j) ) 

V;<J{p) = 6,,~CP^, (6.59) 

we obtain 

n ^ / ,4 G{y) G{z) {G{y) G{z) Zl{L))-'^-'/(''\ 

(27r)4^^ i % + M^iy) z + M\z) {Z{y) Z{z) Zl{L)YI^ 

xA~\y)A-\z) [M±A^ (w,{x, y, z)A{y)A{z) + W^{x, z)B{y)B{z)) 

{W^{x, y, z)A{y)A{z) + W^{x, y, z)B{y)B{z)) 
{W^{x, y, z)A{y)B{z) + We{x, y, z)B{y)A{z)) 



\ 2 




1 Aiy) - 


A{z) 


2{y- 




, B{y) - 


Biz) 



y- z 

(6.60) 



with the kernels 



Wi{x,y,z) = ir^ + z{— ttt ~ o n ^ ' (6-61) 



\^ 3^; 3a;2/ 3 3x 3^^ 

W^2(x,2/,;^) = (6.62) 

M/r , C^' + /I , (2C-6)y Cy' 

W^3(a;, y,z} = TT^-Z — + 77^ + ^ 77 



3a;2 \ 3x 3x2/ ^^.2 

^y__{C±iWM^ (6.63) 

3 3x 3x^ 

-2(2^ , / 4 , 4C!/\ 2 , 4a 2iy' 



W,{x,y,z) = _ + + + (6,64) 

m.y.) ^ g-.(l±i.|^|Vi.^.g. (6.65) 



33/^ Y 3^ 3^^ J 3 3^ 3^^^ 



3x2 \^ 3x2 y 3 3^, 3^,2 
2 2 

W-r{x,y,z) = -- + z + y--y, (6.67) 

X X 

iy8(:^;,?/,^) = 2(-^ + ^) . (6.68) 
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Note that the symmetry factor 1/2 and a factor 1 / (3a;) from the left hand side of the gluon 
equation have already been absorbed in the kernels. From this expression the correspond- 
ing one for the bare vertex construction can be read off easily by setting = and 
replacing the remaining factor {A{y) + A{z))/2 in eq. ()6.60|) by unity. 



6.3.1 Ultraviolet analysis of the quark- loop 

In section 13.3.11 we encountered quadratic divergences in the gluon equation for C 7^ 4. 
These quadratic divergences occur in the ghost- and gluon-loop of the gluon equation and 
show up in the quark-loop as well. To identify such terms we expand the dressing functions 
in the integrand of the quark-loop around large loop momenta y with the difference {z — 
y) still larger than any quark mass. To leading order this expansion amounts in the 
replacements 

Giz) ^ Giy), 
A{z) A{y), 



Ajy) - A{z) 

y- z 
B{y) - B{z) 



B'iy), 



y- z 

(Aiy) - Aiz))iy + z) A'iy)iy + z) 



2{{y - zf + {M\y) + M\z)f) ' 2{y - z) ' ^^'^^^ 

with the derivatives A' and B' . Note that the first two equations are identical to the 
angular approximation employed previously in subsection 13.3.11 For large momenta x and 
z the denominators in eq. (jfj.fjOj) simplify and the angular integrals are trivially performed 
using the integrals given in appendix 1(111 We arrive at 



3 3a; \ 3?/ 3x 



Keeping in mind a factor (1/y) hiding in the derivatives we are now able to identify three 
quadratically divergent terms: (4 — C)/3a; in the second line, —2(4 — ()y/3x in the third 



108 



6.3. Incorporating the quark-loop in the gluon equation 



line and Ay/x in the last line. The first two of them are proportional to (4 — Q and 
reminds us of similar terms occurring in eq. These terms are artefacts of the regu- 

larisation and will be subtracted from the kernels. However, we encounter the additional 
^-independent quadratic divergent term Ay/x originating from the transverse part of the 
Curtis-Pennington vertex. Such a term is already known from corresponding studies in 
QED |169j . Although first suggestions have been made how the Curtis-Pennington vertex 
should be modified to avoid this problem |139j . a convincing solution has not been found 
yet. In this thesis we therefore choose the pragmatic strategy of subtracting this term by 
hand together with the other quadratically divergent parts. 

Moreover we subtract all further terms proportional to {A — Q. Although these terms 
are not quadratically divergent they are artefacts of the regularisation. We then obtain 
a ^-independent expression for the quark loop at large momenta. Together with the 
corresponding expressions for the ghost-loop and the gluon-loop, eq. we obtain a 

transversal gluon propagator in the ultraviolet as it should be the case in Landau gauge. 

Collecting all modifications together we arrive at the new kernels 



Wi{x,y,z) = Wi{x,y,z) , (6.71) 

bx 

W2{x,y,z) = 0, (6.72) 

Ws{x,y,z) = W,{x,y,z) + '^^y^^~^\ (6.73) 

W^{x,y,z) = W^{x,y,z)-^ ^ ^, (6.74) 

6X 

W,{x,y,z) = W,{x,y,z)- ^y^^^J'^~^\ (6.75) 

W,{x,y,z) = W,{x,y,z)-^ ^ ^, (6.76) 

bx 

W,{x,y,z) = Wr{x,y,z)- ^y^^y^'\ (6.77) 

X 

Ws{x,y,z) = Ws{x,y,z). (6.78) 



Note that the subtracted terms are chosen to preserve the symmetry of the kernels with 
respect to the squared quark momenta y and z. 

Without quadratic divergences we are in a position to extract the leading logarithmic 
divergence of the quark-loop. With modified kernels the ultraviolet limit of the quark-loop 
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is given by 



-2 , A'{y) fl , -2^.^ 



+B'{y)M{y)- +^ (^-1 + -M'{y)j j . (6.79) 

Similar to the situation in the DSE for the quark mass function, c.f. subsection 16.1.21 the 
leading ultraviolet term is the first term in the curly brackets. Substituting the ultraviolet 
limit of the vector self energy, eq. ()6.22j) . the ghost and gluon dressing functions, eqs. ()6.29j) . 
and choosing the perturbative renormalisation condition G{s) = Z{s) = 1 we arrive at 



n 



uv _ 2iV^ g' 



quark 3(2^+1)^16^2 



cjlog(-j+l - ^tulog (^-j + 1 



(6.80) 



Here we have u = /?o«('S)/(47r) = (llA'c — 2Nf)a{s)/{127T). Moreover 6 is the anomalous 
dimension of the ghost propagator which is related to the corresponding anomalous dimen- 
sion of the gluon by 7 + 25 + 1 = 0, c.f. eq. ()3.34p . Combining the expression ()6.80|) with 
the results for the ghost and gluon loop, eq. ()3.33p we obtain the anomalous dimensions 

(6.82) 



UNc - 8Nf 

which are in accordance to one-loop perturbation theory for arbitrary numbers of colours 
Nc and flavours Nf. 

6.3.2 Infrared analysis of the quark- loop 
In section 13.3.21 we employed the power law ansatz 

Z{x) = Ax^^, G{x) = Bx-'', (6.83) 

for the ghost and gluon dressing functions at small momenta x. Substituting this ansatz 
into the gluon equation we have found the ghost loop to be proportional to x~'^'^, dominat- 
ing the gluon-loop in the infrared. When we compare our expression for the quark-loop, 
eq. ()6.6()|1 . with the ghost loop in eq. ()3.25|1 we find two ghost dressing functions, G{y) and 
G{z), in each loop respectively. However, as the momenta y and z in the denominators of 
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6.4. Numerical results 



the quark loop are negligible compared to the quark masses in the infrared we anticipate 
that the quark-loop is less divergent than the ghost-loop for small momenta, provided the 
parameter d is small. 

An explicit calculation along the lines of our analysis in subsection 13.3.2! and appendix 
1(^.3! shows that the quark loop is proportional to 3;-2K+2+K(i/<5 -^^ j^^ie infrared. Therefore the 
quark loop is suppressed for small momenta provided the parameter d fulfils the condition 

d < (6.84) 

K 

As we have k ^ 0.5953 and 6 = —1/4 for Nc = 3 and Nf = 3 we find the condition 
d < 0.84, which is satisfied for all quark-gluon vertices employed in our calculation. From 
a numerical point of view we encounter serious instabilities in the quark and the gluon 
equation once d is taken to be larger than d ^ 0.2. 

We conclude that the quark-loop does not change the infrared behaviour of the ghost 
and gluon dressing functions found in section in.3.2[ In pure Yang-Mills theory as well as in 
QCD we thus have an infrared finite or vanishing gluon propagator and a ghost propagator 
which is more divergent than a simple pole. The Kugo-Ojima confinement criterion and 
Zwanziger's horizon condition are both fulfilled not only in pure Yang-Mills theory but 
also in QCD. This is a central result of this chapter. 



6.4 Numerical results 

The numerical treatment of the integrals in the quark, ghost and gluon equations has been 
described in subsections 13 . 4 . 1! and 16 . 2 . 21 The iteration process is done for the ghost-gluon 
system and the quark equations separately: we first iterate the Nf mutually uncoupled 
quark systems until convergence is achieved, feed the output into the ghost and gluon 
system, iterate until the ghost-gluon system converges, feed the output back into the quark 
equations and so on, until complete convergence of all equations is achieved. Similar to 
the quenched calculations we used a{s) = 0.2 at the renormalisation point s = yU^, and a 
transverse tensor to contract the gluon equation, C = 1- 

In contrast to section 16.21 we fix the physical scale of the system not by the condition 
q;(M|) = 0.118 but by adjusting the pion decay constant to the experimental value. This 
choice has an important advantage: Whereas the asymptotic behaviour of the running 
coupling depends strongly on Nf via the coefficient Pq {c.f. appendix IA.2|) . the pion decay 
constant turns out to be almost independent of Nf. In order to be able to compare different 
constructions for the quark-gluon vertex we fix the scale by requiring = 85 MeV for the 



Chapter 6. The coupled system of quark, gluon and ghost DSEs 



111 





M(0) 




(_(^^))l/3 


a{Mz) 


A MOM 




[MeV] 


[MeV] 


[MeV] 






[MeV] 


vertex 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


qu. 


unqu. 


bare d=0 


172 


169 


37.4 


37.2 


155 


162 


0.116 


0.144 


694 


718 


CP d=6 


146 


153 


49.1 


49.8 


219 


215 


0.116 


0.138 


694 


578 


CP d=0 


358 


346 


76.5 


76.2 


291 


295 


0.116 


0.141 


694 


644 


CP d=0.1 


451 


425 


85.0 


85.0 


321 


325 


0.116 


0.143 


694 


669 



Table 6.2: A comparison between the quenched (qu.) and unquenched (unqu.) results for 
the quark mass M(0), the pion decay constant the renormalisation point independent chiral 
condensate, the running coupling at the mass of the Z-boson and ^^qSd different vertices 
and values of the parameter d. The unquenched calculations are done for Nj = 3 chiral quarks. 
Furthermore we have S = —9Nc/ {AANc — 8Nf) = —0.25 in the present case. 



CP-construction with d = 0.1 and keep this scale when employing other vertices^. 

In table 16.21 we compare results for the quenched and unquenched system of equations^°, 
c.f. section 16.21 The quark mass, the pion decay constant and the chiral condensate differ 
only slightly for each vertex construction respectively. The only sizeable difference occurs 
in the running coupling. As expected from perturbation theory the unquenched running 
for Nf = 3 results in larger values of the running coupling at = (M^)^ compared to the 
quenched case Nf = 0. We obtain a{Mz) ~ 0.140, which is somewhat larger than usually 
quoted values from experiments^. If we increase the number of flavours in our calculation 
we encounter large numerical uncertainties and do not obtain convergence for Nf > 5. 

All employed vertex constructions allow for nontrivial solutions of the quark equation 
corresponding to dynamical chiral symmetry breaking. However, similar to the quenched 
case the bare vertex construction and the CP-type vertex with d = 6 generate much 
too small quark masses compared with typical phenomenological values^^. For d = we 
obtain good results for the quark mass, the pion decay constant and the chiral condensate, 
whereas the choice d = 0.1 leads to somewhat large values. It is interesting to note that 

^Recall from subsection 16.2.11 that the approximation H6.36|l for the pion decay constant does not 
include the effects of subdominant pion Bethe-Salpeter ampHtudes. Therefore the value = 85 MeV is 
chosen to allow for corrections compared to the experimental value /jr = 93 MeV. 

^°The scale of the quenched results has been transformed to our modified scale. Both scales are very 
close to each other as can be seen by comparison with table 16. l1 

^^However, such large values are not yet excluded by experiment. A recent analysis of experimental 
data from r-decay suggests a{Mz) ~ 0.129 |ll)6| . 

^^This is in accordance with previous quenched |145| and unquenched jl46j investigations of these cases 
in a truncation scheme with angular approximations in the Yang-Mills sector. 
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6.4. Numerical results 




x=p^[GeV^] 




x=p^[GeV^] 




Figure 6.9: Displayed are the ghost and gluon dressing function, Z and G, the running coupling 
a, the quark mass function M and the inverse vector self energy XjA. The calculations are 
done quenched and unquenched with Nj = 3 quarks in the chiral limit. 



d = of all values is preferred as in this case the quark equation resembles most the 
fermion equation of QED. 

In Fig. 16.9! we display the ghost, gluon and quark dressing functions corresponding 
to the unquenched and two representative quenched cases in table 16.21 We find different 
anomalous dimensions in the ultraviolet corresponding due to the change from Nf = to 
Nf = 3, c.f. eqs. fl6.29|l . ()6.32|1 . ()6.82j) . As expected from the infrared analysis in subsection 
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Figure 6.10: Results from the unquenched calculation with Nf = 3 massive quarks. We used 
the renormalised masses m„/rf(2GeV) ^ 4.5 MeV and ms(2GeV) ^ 95 MeV. 

If). 3. 21 the back- reaction of the quark- loop in the gluon equation does not affect the infrared 
behaviour of the ghost and gluon dressing functions. Consequently the infrared fixed point 
of the running coupling is the same as in pure Yang-Mills theory. Thus the Kugo-Ojima 
confinement criterion and Zwanziger's horizon condition, c.f. section are satisfied in 
quenched and unquenched Landau gauge QCD. 

Our results for the case of explicitly broken chiral symmetry are shown in Fig. 16.101 We 
choose Nf = 3 with renormalised quark masses corresponding to rriu/di!^ GeV) = (4.4—4.6) 
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6.4. Numerical results 



MeV and ms(2GeV) = (90 — 95) MeV. These masses are well in the range suggested by 
the Particle Data Group |l()2j . Compared to the chiral case the behaviour of the ghost 
and gluon dressing functions hardly changes. For the quark mass function we obtain the 
irregular asymptotic solution in the ultraviolet as expected from the analysis in subsection 

EI2 

Furthermore we provide fits to our results for the quark propagator employing the fit 
functions^^ 



m 



x/Alcnr , (6.85) 



[AiO)]-' + h (x/AIcd) + h2 (x/A 



2 



m 



with X = and the six parameters gi, g2, hi, h2, h^, h^. We used the renormalisation point 
independent current-quark mass rh, which is related to the renormalised mass M{s) by 

Mis) (^IHs/AIcd])'^ , (6.87) 

to one loop order. For the running coupling, the ghost and the gluon dressing function we 
use the form 'Fit B', given in eq. ()3.58|) and the fit functions from eqs. ()3.59p . In table IU7^ 
we give our values for all parameters as well as the numerical results for M(0) and [^4(0)]"^. 
Note that the scale Aq§^^ is different to the corresponding scale in the chiral limit due 
to the different ultraviolet behaviour of the quark-loop when quarks with non-vanishing 
bare masses are employed. When plotted the fits are virtually indistinguishable from our 
results in Fig. I6.1U1 

Unquenched lattice calculations employing dynamical quarks are complex and time 
consuming |17()j . To our knowledge such simulations for the propagators of QCD have not 
yet been performed. From our results in the Dyson-Schwinger approach we do not expect 
drastic differences between quenched and unquenched propagators on the lattice. 

Finally, we investigate positivity violations in the gluon and quark propagators. Recall 
from our discussion in subsection 12.3.21 that the condition 

oo 

j dtdt' fit',p} Si-it + t'),p} fit,p} <0 (6.88) 







^•^For simplicity we did not include the regular asymptotic term in the fit function for the quark mass 
and we did not care about logarithmic corrections to the ultraviolet behaviour of the wave function 
renormalisation. Nevertheless both fit functions work very well. 
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A M OM 
^^QCD 

[MeV] 


a 


b 


c 


d 


mR 
[MeV] 


rh 
[MeV] 


M(0) 
[MeV] 


9i 


92 


CP d=0 


600 


1.35 


1.04 


1.36 


1.81 


4.4 


4.6 


354 


2.30 


1.29 














90 


96 


507 


3.42 


1.09 


CP d=0.1 


620 


1.32 


1.06 


1.42 


1.78 


4.5 


4.6 


429 


2.22 


1.28 














95 


96 


564 


3.80 


1.14 





[MeV] 




h 


h2 


h3 


/i4 


CP d=0 


4.4 


0.638 


0.503 


5.263x10- 


-3 


0.544 


5.202x10-3 




90 


0.671 


0.305 


1.420x10- 


-3 


0.321 


1.402x10-3 


CP d=0.1 


4.5 


0.538 


0.448 


2.740x10- 


-3 


0.477 


2.706x10-3 




95 


0.590 


0.303 


1.051x10- 


-3 


0.317 


1.037x10-3 



Table 6.3: Parameters for the fits to the unquenched results with Nj = 3, a{s) = 0.2, 
S = -0.25, 7™ = 12/27 and Po = 27/3. 

is sufficient for a given propagator S to violate the reflection positivity axiom of Euclidean 
quantum field theory. Here / are complex valued test functions, c.f. subsection 12.3.21 The 
one-dimensional Fourier transform S{t,'p) of the propagator S{po,p) is given by 

Sit,p):= I ^^(po,p-)e^^«*. (6.89) 

Provided there is a region around to where S'(— to, p) < one can choose a real test function 
f{t) which peaks strongly at to to show positivity violation. In Fig. 16.111 we display the 
Fourier transform of the nontrivial part D{p^) = Z{p^)/p'^ of the gluon propagator and the 
vector part <Jv{p'^) = A{p'^) / {p"^ A{p'^) + Bij)^) of the quark propagator (we chose = 0). 
Clearly the Fourier transform of the gluon propagator is negative on a large interval. The 
resulting positivity violation for the transverse gluon propagator in Landau gauge can be 
interpreted as a signal for gluon confinement, c.f. subsection 12.3.21 Note that the ghost 
propagator violates positivity trivially for G{p'^) > as can be seen from the definition 
in eq. ()3.3|1 '^^. From the Fourier transforms of ayiv^) displayed in the second diagram 
of Fig. 16.111 we cannot conclude positivity violation for the quark propagator. As we 
have an indefinite metric in the unphysical part of the state space of QCD this result 
is not in contradiction with quark confinement. In addition we have to keep in mind 

Already the bare ghost propagator violates positivity. This is not surprising due to the unphysical 
spin-statistic relation of the ghost field. 
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gluonpropagator 



quarkpropagator: 
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Figure 6.11: Here we display the one dimensional Fourier transforms of the gluon propagator, 
D{—t,p^), and the vector part of the quark propagator, cTy(— We observe violation of 
reflection positivity for the gluon propagator but not for the quark propagator. 



that eq. 



is a sufficient but not a necessary condition for the violation of positivity. 



Further investigations of the quark propagator will be done in future work. 



6.5 Summary 

In this chapter we have presented solutions of the (truncated) Dyson-Schwinger equations 
for the propagators of Landau gauge QCD. We first concentrated on the Dyson-Schwinger 
equation for the quark propagator. We proposed several ansdtze for the quark-gluon 
vertex which consist of an Abelian part carrying the tensor structure of the vertex and 
a non-Abelian multiplicative correction. Our guiding principles for the construction of 
these vertices have been two important conditions on the truncated quark equation: it 
should be multiplicatively renormalisable and recover perturbation theory for large exter- 
nal momenta. We showed that the resulting quark mass function is independent of the 
renormalisation point and has the correct asymptotic behaviour for large momenta. 

In the quark equation both the ghost and gluon dressing function show up at least 
implicitly. In the quenched approximation we employ our solutions of the ghost and gluon 
Dyson-Schwinger equations in the truncation scheme of chapter El In the unquenched 
calculations we include the back-reaction of the quarks on the ghost and gluon system and 
solve the quark, gluon and ghost Dyson-Schwinger equations self-consistently. 
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All our solutions show dynamical chiral symmetry breaking. However, only carefully 
constructed vertex ansdtze have been able to generate masses in the typical phenomeno- 
logical range of 300 — 400 MeV. Constructions with an Abelian part satisfying the Abelian 
Ward-identity are superior to other vertex ansdtze. We obtained very good results for 
the quark mass, the pion decay constant and the chiral condensate by employing a gener- 
alised Curtis-Pennington vertex. In the chiral limit both, the quark mass function and the 
vector self energy are close to recently obtained lattice results. This agreement confirms 
the quality of our truncation and in turn it shows that chiral extrapolation on the lattice 
works well. 

In the unquenched case including the quark-loop in the gluon equation with Nf = 3 
light quarks we obtain only small corrections compared to the quenched calculations. 
In particular the quark-loop turns out to be suppressed in the gluon equation for small 
momenta. We thus showed on the level of our truncation that the Kugo-Ojima confinement 
criterion and Zwanziger's horizon condition are satisfied in Landau gauge QCD. 

Furthermore we searched for positivity violations in the gluon and quark propagators. 
We confirmed previous findings that the gluon propagator shows violation of refiection 
positivity. Thus the gluon is not contained in the physical state space of QCD. We did 
not find similar violations for the quark propagator. It is an open question whether these 
violations exist and if not, how else confinement shows up in the quark propagator. 



Chapter 7 

Conclusions and outlook 



The central objects of interest in this thesis have been the two-point Green's functions of 
QCD, i.e. the ghost, gluon and quark propagators. We have presented solutions for the 
corresponding set of coupled Dyson-Schwinger equations, employing ansdtze for the ghost- 
gluon vertex, the three-gluon vertex and the quark-gluon vertex. These ansdtze have been 
constructed such that important constraints from general principles are satisfied: both, 
the running coupling and the quark mass function are independent of the renormalisa- 
tion point. Furthermore, we obtained the correct one- loop anomalous dimensions for all 
propagators. 

We have been able to achieve considerable progress by overcoming the angular approx- 
imations used in previous studies. These approximations have proven to be quantitatively 
unreliable for small momenta. Our solutions corroborate previous findings of a vanish- 
ing gluon propagator in the infrared, whereas the ghost propagator is highly singular at 

= 0. In the momentum range covered by lattice simulations the lattice data and our 
results agree very well even on a quantitative level. Differences occur for intermediate 
momenta where the Dyson-Schwinger equations suffer from neglecting genuine two-loop 
diagrams. An important aim in future studies is the inclusion of these diagrams. 

For the running coupling we obtain a fixed point in the infrared. Together, an infrared 
finite coupling and an infrared vanishing gluon propagator disprove the old idea of infrared 
slavery. Instead the ghost propagator turns out to be the dominant degree of freedom in 
the infrared, at least in Landau gauge. This is in accordance with Zwanziger's horizon 
condition and the Kugo-Ojima confinement scenario. A strongly infrared diverging ghost 
propagator signals an unbroken global colour charge. In the Kugo-Ojima scenario this 
charge is used to show that the physical states of QCD are colourless. However, as two 
other central assumptions of the scenario, the existence of a well defined BRS charge and 
the violation of the cluster decomposition property, are currently unproven a complete 
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verification of tlie Kugo-Ojima scenario is still lacking. 

First steps towards a solution of the Dyson- Schwinger equations for the ghost and 
gluon propagators in general gauges have been performed in this thesis. We investigated 
the infrared behaviour of the ghost and gluon dressing functions in these gauges employing 
a bare vertex truncation. In all linear covariant gauges we find identical results to Landau 
gauge. For general ghost-antighost symmetric gauges, however, we do not find power 
solutions in the infrared when bare vertices are employed. 

Another piece of progress achieved in this thesis is the solution of Dyson-Schwinger 
equations on a torus, i.e. for periodic boundary conditions for the fields. For various 
truncation schemes these solutions suffer only mildly from finite volume effects in the in- 
frared and arc very close to the continuum results for intermediate and large momenta. 
From a numerical point of view we found it to be easier to obtain solutions on a torus 
than for continuous momenta. The reason is the property of the finite volume torus to 
act as regulator in the infrared. Thus one avoids all problems with infrared singularities 
encountered in the continuum formulation of the Dyson-Schwinger equations. Our treat- 
ment of the Dyson-Schwinger equations on a torus is an ideal starting point for further 
investigations, provided the role of the zero modes can be clarified further. First, we are in 
a position to study finite volume effects in more detail by comparing solutions on different 
volumes to the continuum results. Second, by varying the extension of one direction in 
space-time the inclusion of finite temperature effects might be accomplished in a relatively 
easy way. Here the main qualitative question arises about the fate of the Kugo-Ojima 
confinement criterion at the deconfinement transition. Finally, changing from periodic 
to twisted boundary conditions we hope to be able to include topological effects in the 
Dyson-Schwinger equation approach. 

The quark propagator is the basic input in many phenomenological models which de- 
scribe mesons and baryons as bound states of quarks and gluons. In this context the 
Dyson-Schwinger equation for the quark propagator is the link between underlying QCD 
and the model frameworks. By solving the unquenched set of quark, gluon and ghost 
Dyson-Schwinger equations in different truncation schemes we have provided an impor- 
tant step in connecting these models with the fundamental theory. We obtained dynamical 
chiral symmetry breaking in the quark equation and found masses in the order of phe- 
nomenological values. Similar to previous studies in quenched QED we had to build the 
quark-gluon vertex carefully along general principles to obtain satisfying results. A central 
aim in future work will be to extend these calculations to the corresponding Bethe-Salpeter 
equation for mesons and thus make further contact with experiment. 
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Appendix A 

Notations, conventions and 
decompositions 



A.l Euclidean space conventions 

In Euclidean space-time we use the metric g^u = 5^^ and Hermitian Dirac matrices 7^ 
related to the standard Minkowski ones PP by 

1' = -^I'u^ J = 1-3, (A.l) 

1' = ill- (A.2) 

The Dirac matrices satisfy the Clifford algebra 

{7^,74 = 25^,. (A.3) 



A.2 Perturbative running coupling 

The non-perturbative definition of the running coupling in Landau gauge is given by (c./. 
eq. (IT^ in section ITT^ : 

a{p') = a{^^'')Z{p\^^'')G\p^^^^). (A.4) 

Here Z is the dressing function of the full gluon propagator, G the corresponding dressing 
function for the ghost propagator and /i^ is the squared renormalisation point. Since we 
have the renormalisation condition 

Z{^\^^)G\^\^^) = l, (A.5) 
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A. 3. Definitions and decompositions of correlation functions 



the momentum dependence of the running couphng is the same as the dependence of the 
renormahsed couphng a{fi^) = g'^{fi'^)/47T on the renormahsation point /i^. 

For large renormahsation points /i^ the behaviour of can be calculated from 

perturbation theory. To three-loop order the Particle Data Group gives the expression 

rrn2i 



1 - 



2/5i In [HfiyAlcD 



PI 



X 



In [ln(/iVA^cD)] - 



+ 



8P! 



(A.6) 



in the MS renormahsation scheme. The coefficients of the /^-function are defined by 



da 


= 2(3{a 


N Po 2 


47r2 


P2 4 


(A.7) 




= 11 - 


3 ^' 






(A.8) 




= 51- 


—Nf , 
3 ^ ' 






(A.9) 






5033 ,, 


325 




(A.IO) 


P2 


= 2857 


Nf + 

9 * 


27 ^ ' 





where N^. is the number of colours and Nf is the number of fiavours. 



A. 3 Definitions and decompositions of correlation func- 
tions 

In this section we give definitions and conventions for some correlation functions needed 
in this thesis. Further definitions are given at the appropriate places in the main body of 
this thesis. 

A. 3.1 Ghost, gluon and quark propagators 

The full ghost, gluon and quark propagators in coordinate space are defined as 
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The inverse bare propagators in coordinate space are easily derived from the quadratic 
part of the action S = J dzC with the Lagrangian C given in eq. ()2.1()|1 . One obtains 

6^S 



Df''\x-y)_ 
[S^'^'^'ix - y)Y' 



5c''{x)d'{y) 
6^S 



(A.14) 



5^S 



1 



SK(^)A<-M = ' \ + ( 1 - ^ ) 1 - . (A.15) 



+ mQ)5{x-y) 



(A.16) 



6^!'^{x)^!^{y) 

with the gauge parameter A and the bare quark mass itiq. After Fourier transformation 
the corresponding expressions in momentum space are given by 

-1 



Ay 



-1 



(A.17) 
(A.18) 
(A.19) 



A. 3. 2 The ghost-gluon vertex 

The tree level ghost-gluon vertex, y^I^^""^^ , is derived from the ghost-gluon part of the action, 

S,host-,iuon = I d'x' f ) gr'' {d^-n A^c' + i^g r'^c'' A^d^ c'] . (A.20) 

For general values of the gauge parameters a and A the vertex is given by 

■p(0)a6c/ \ _ " ^ghost-gluon 

^ ^ ' 5Al{x)5c\y)5c'^{z) 

= -gr'' [^ (l - f ) {d',5\z - y)) 5\z - x) + i^d^ [b\z - y)5\z - x)) . 

(A.21) 

Using the momentum conventions of Fig. lA.ll the Fourier transformed bare ghost-gluon 
vertex reads 



r^;^'''%k,p,q) 



(f[xyz\ Tf'{x, y, z) e^(*^-+5-?'-P-) 
gr^^{27,)U\k + q-p) 



a\ a ' 



(A.22) 



where the abbreviation d^x d^y d^z =: d^[xyz] has been introduced. Note the symmetry of 
the vertex in the ghost momenta p^ and if a = 1. For convenience we define a reduced 



124 



A. 3. Definitions and decompositions of correlation functions 



V 



V / ki \ 

b c / ^ \ 

a I 

Figure A.l: Momentum routing for the tree level ghost-gluon and four-ghost vertices. 

ghost-gluon vertex function T^/^\p,q) by 

rW'^^^(A:,p,g) = gr'^i2n)U\k + q-p)T^;\p,q), 

Furthermore we will need the full one-particle irreducible ghost-gluon vertex in coordinate 
space, given by 

T"!''{x, y, z) = ,,,,,,,,, , (A.24) 
A. 3. 3 The four-ghost vertex 

The bare four-ghost vertex Tf^l^'"^^ is derived from the four-ghost part of the action, 

54,. = / d'x' { f (l - f ) ^g^f-f'^-c^cW] . (A.25) 



One obtains 

-r^(0)abcd / 

I \ T II 7 111 

With the momentum conventions given in Fig. lA.ll the Fourier transformed bare four-ghost 
vertex is given by 

^f^t'ih, k2, ks, k,) = f (l - f ) Xg^r''r'%27:f5\k, + k,-ks- k,). (A.27) 
Again we define a reduced vertex function F^^^^ for convenience: 

Tf^t'\hMMM) = g^r''r''\2nY6\k, + k^-k,-k,)Tf 



.(0) ^ ^ [t^ 
Agh 2 V 2 



rSl. = Sfl-S)A. (A.28) 
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The full four-ghost vertex in coordinate space is formally given by 
A. 3. 4 The other vertices 

The quark-gluon vertex, the three-gluon vertex and the four-gluon vertex do not show 
up in the derivation of the ghost Dyson-Schwinger equation in appendix iBl We therefore 
refrain from giving the expressions for these vertices in coordinate space but merely state 
the corresponding expressions for the bare vertices in momentum space. All vertices given 
in this subsection are reduced vertices i.e. we have suppressed a factor {2-nY and a delta 
function for momentum conservation. 

The bare quark-gluon vertex in momentum space is given by 

Vf\k,q)=ige^,, (A.30) 

where is a generator of SU (Nc) gauge transformations. 

For the bare three-gluon vertex in momentum space one finds 

^^}f%P, Q, k) = gr'^ {5,, {p - q), + 5., {q - k), + 5,, {k - p).) , (A.31) 

where /"^'^ is the structure constant of the gauge group. The three-gluon vertex is com- 
pletely symmetric with respect to the three gluon momenta. 
The four-gluon vertex is given by 

ri.lT(^i' h, h) = g' {rT%e. - + r'j'^^. - 5,^ 5,p] 

j^jebdjeac]^^^^ 5a/3 - 5ua ^^if}] } , (A. 32) 

and is symmetric in the Lorentz indices of the four attached gluons. 



A. 3. 5 Decomposition of the connected ghost-gluon correlation 
function 

For the derivation of the ghost Dyson-Schwinger equation from the generalised Lagrangian 
(I2.11|) performed in appendix^lwe will need some decompositions of connected correlation 
functions. These decompositions are derived in this and the next subsection. We will use 
the generating functional of connected Green's functions, W , and the effective action, F, 
see appendix IE) for details. Furthermore we will exploit the fact that the fields and sources 
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can be written as functional derivatives of W and F: 



m _ _ SW_ 

Sa ' 5a 

sr _ _ sr 

5c ' 5c 



= c, 



(7, 



5W 

5r 



— 



J,. 



(A.33) 



The sign conventions have been chosen such that derivatives with respect to c and a are 
left derivatives whereas the ones with respect to c and a are right derivatives. 
With the help of the matrix relation 



5(t> ^ 5(t>^ 



and the identity 



I 



4 5a\y) Sc'^iz) 
a z- 



(fz 



5^W 



5c'^{z)5c^{y) 5a"'{x)5a'^{z) 



(A.34) 



(A.35) 



5c'^{z) 5a°'{x) 

we decompose the connected ghost-gluon correlation function, {A'^^{x)'^{y)c'^{z)) , in the 
following way: 

5^W 



{Al{x)c\y)c\z)) 
5 



5J^{x)5a^{y)5a''{z) 



5J-{x) 



5''V 



5&{y)6c'iz) 



-1 



,4 . SAiim) 6 



d Uy 



5J-{x) 6Ai{ui) 
5^W 



6^r 



d [^11*21*3] 



5c''{y)6c''{z) 
6^W 



-1 



5^T 



-I 

~ I ^"''"'^"''''5J^{x)5J;^{ui) 5a^{y)5a''{u2) 5A'^{ui)5c''{u2)5cf{u3) 5af{u3)5a''{z) 

= J d^[uiU2U3]D;t(x - Ui)D'^(u2 - y)rt'^(ui, U2, Us)D'J(us - z). (AM) 

Here we used the abbreviation d'^[uiU2U3] := d'^ui d'^U2 d'^u^ and the definitions of the gluon 
propagator D^^, the ghost propagator Dq and the ghost-gluon vertex given in previous 
subsections. 



A. 3. 6 Decomposition of connected four-ghost correlation func- 
tion 

Furthermore we need the decomposition of the four-ghost correlation function into one- 
particle irreducible parts. We start at a stage where the sources are still present and set 
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them to zero at the end of the derivation. We first give the decomposition of the connected 
ghost-antighost-ghost three-point function 



5 



5^W 



5a\y) 



2-r n-i 



6'T 



4 SAKm 
a ui 



6 



6a\y) 6Al{ui] 
5^W 



5^V 



d [uiU2Uz] — 



6^W 



-1 



5^T 



6^W 



(A.37) 



Then we decompose the connected four-ghost Green's function: 



c''{x)(f{yy{z)(^\w)) 



5^W 



6a''{x)5a^{y)6a''{z)6a'^{w) 

/ ?/,r,7/,0 I 

5a 



5^W 



X 



\y)6J;{ui) 6a'{z)6af{u2) 
5^T 5^W 



(A.38) 



Carrying out the remaining derivative gives four terms. The two terms where the derivative 



acts on the second and on the last propagator vanish, because the term 



5W 



vanishes 



5a>'{y)SJt{ui) 

when the sources are set to zero. The contribution where the derivative acts on the first 



propagator can be treated using eq. ()A.36|) . In the expression with the derivative acting 
on the vertex we use 



5^W 



5^V 



6a^{y)6J^{ui) Sa''{x)6Al{ui)6cf{u2)6c9{u3) 6a<'{x)6a^{y)6cf{u2)6c9{us) 



5^W 5^V 



5^W 5^W 5^T 



- j d'^[uiU^] 



5a''{x)a''{u^) 5a^{u5)5a^{y) 6c9{u3)6c^{u5)6c''{u4)6cf{u2)' 

(A.39) 
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Collecting all this together we arrive at 



6^r 5^W 5^W 



X 



6At{ui)6c9{u^)6c''{uQ) 5a^iuQ)5a\y) 5a^{z)5a'{u2) 
5^W 5^W 



j d!^[uiU2U:>,UiU^] i — 



5^W 5^T 



X T-^7-Y7^ ■ (A.40) 



Interchanging some Grassmann fields in the correlations and using the definitions for the 
propagators and vertices given in the previous subsections we arrive at 

{c\y)c''{w)c''{xy{z)) = J d*[uiU2U3U4U,ue]{D;i{ui-U4)D^J{x-U5) 

xTf^^u,, u,, U5)D^\ue - y) D'^{z - U2) 
xrf{uuU3,U2)Di,\m-w)} 

+ J d^[uiU2U3U4U5] {D^%x - U4) D^\u5 - y) 

xD'Jiz - U2) r1lt/{u5, Us, M4, U2) Dg{u3 -w)}, (A.41) 
which is the decomposition of the four-ghost correlation used in appendix El 
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The derivation of the ghost 
Dyson-Schwinger equation in general 
ghost-antighost symmetric gauges 

In this appendix we will derive the ghost Dyson-Schwinger equation from the generalised 
Lagrangian ()2.11|) . As there are no direct couplings between quarks and ghosts contained 
in the Lagrangian it is sufficient to employ only the ghost-gluon part. We start by trans- 
forming the ghost-gluon part of the Lagrangian into a more suitable form using partial 
integrations under the assumption of vanishing fields at infinity. We obtain 

+i (l - I) gr'^c^d, {A^c") + i'^gr''^eAld,c\ (B.l) 
The partition function of the theory is given by 

Z[J, a, a] = j V[Acc] exp |- ^ c/S £ + j d^z (AV" + ac + ca) | (B.2) 

with the sources J, a and a of the gluon, antighost and ghost fields, respectively. The 
action is given by iS[J, c, c] = j d'^z C. The generating functional of connected Green's 
functions, W[J, a, a], is defined as the logarithm of the partition function. The functional 
Legendre transform of W is the effective action 

T[A, c, c] = -W[J, a,a]+ I d'^z {A^r + ac + ca) , (B.3) 



129 



130 



which is the generating functional of one-particle irreducible vertex functions. The fields 
and sources can be written as functional derivatives of the respective generating functionals 
in the following way 



6(7 ^ 6a ^ 6J^ 



The sign conventions have been chosen such that derivatives with respect to c and a are 
left derivatives whereas the ones with respect to c and a are right derivatives, 

5 6 
. , := left derivative — := right derivative. (B.5) 

6 [a, c) 6 [a, c) 

Given that the functional integral is well-defined, the Dyson-Schwinger equation for 
the ghost propagator is derived from the observation that the integral of a total derivative 
vanishes. We take the derivative with respect to the antighost field and obtain 

= y V[Acc] ^ exp I - y rfS £ + j (fz (A + ac + ca] 

V[Acc] (^J_ll^^ + exp |- y + j d^z (AV" + ac + ca)| 

J_^Uj1S^ + a Z[J, a, a]. (B.6) 
6c I 

Now we use the relations ()B.4|) and apply a further functional derivative with respect to 
the source cr^iy). We arrive at 

= {^-^f\y) + ^\^)^\y) + - y)5c!}j Z[J, a, a] (B.7) 

with explicit colour indices and space-time arguments. Setting the sources equal to zero 
we obtain the ghost Dyson-Schwinger equation 

The derivative is easily calculated 

^ - SVW + f(l-|)^//-//M.)c/(.)c»(.) 

(l - f ) gf-^a, (A;{zy{z)) + /jgf'^Al(z)d,c''{z). (B.9) 
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Whereas in the covariant formahsm full and connected three-point functions are the same, 
the four-point correlations have to be decomposed into disconnected and connected parts. 
For the four-ghost correlation function this results in 

{c\y)c''izy{zyiz)) = {c\yy{z)) {c\z)cf{z)) - {c\yy{z)) {c\zy{z)) 

+ {c\y)c\zyizyiz))^^,n.. (B.IO) 

Keeping in mind the Grassmann nature of the ghost and antighost fields we then obtain 

- Siz - y)5,, = d^c^{yY{z)) + f (l " f ) ^g' f^"^ f^'^ {{c^{y)A^)cf {z)c^{^)) 

+ {{c^{yy{z)) (c^(z)c^(z)) - {<^{y)cf{z)) 

+ (l - f ) 9f'^'{<^{y)d, {Al{z)c\z))) + ^9f'''^{Ay)Al{z)d,c\z)) , 

(B.ll) 

where all correlations are connected Green's functions. We now use the relation 

^ ^ J 6c'^{z) 5a-"(x) J 6c^{z)6c^{y) 6a''{x)6a'^{z) 

=: j d^z [D^{z - y)] D^^{x - z) 

(B.12) 

and multiply eq. (jRTT]) with -[D'}^{x - z)]-^ = [{c%z)c''{x))]~\ We arrive at 
[D-^ix-y)]-' = d'5{x-y)S-' 

- f (l - f ) ^9'r''f^'' J d'z [D^Jix - z)]-' {{c\y)c\zyizyiz)) 

+ {c^\yy{z)) {Az)cf{z)) - {c\y)c^iz)) {c\zyiz))} 

- ^ (l - f ) ar'-J d'z - z)]-\<^{y)d, {A^,{z)c\z))) 

- ^^gr'' [ d'z[D^a%x - z)]-'{c\y)Al{z)d,c''{z)). 



2' 



(B.13) 



Before we decompose the connected Green's functions into one particle irreducible ones 
we have to take care of the space-time derivatives. Noting that 



I d%v]d; (Sin - z)) 5{u - v)jj^^^ (B.14) 
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with the abbreviation d ud v =: c? [wf], and 



:9u-rr^ = Auv]S{u- z)6{u-v)——^r^- 



d^iuv] d; {5{u - z)5{u - V)) jj^^^ (B.15) 



we can replace the derivative terms by the bare ghost-gluon vertex defined in appendix 
IA.3I The tadpole term can be treated in the following way: 



= 2 I d'z[D^a'ix - zr'r'jf^^ {{c\yyiz)) {c'izyiz))} 
= 2 J d^[zuv] [D^%x - z)]-^6{z - u) 6{u - v) f^'^jfa^ D§{z - y) Dg{v 
= 2 J d^[uv] 5ix - y) 5{z - u) 5{u - v) f'^jf'^^ D{f{v - u). (B.16) 



u 



Plugging the expressions for the ghost-gluon loop and the one for the tadpole into eq. ()B.13|) 
and using the expression for the bare four-ghost vertex given in appendix IA.3I we obtain 



[^G - VT' = d'6{x - y)5'^' - I d'[uv] rgf ^"(x, u, v, y) D{f{v - u) 

d^uv] 5{z - u) 5{u - v) [Z^^"(x - z)]-^ {c\y)c\z)cf{uy{v)) 
d'lzuv] rW^'^^(^,«,t;) [D^a'ix - z)]-' {c\y)Al{vy{u)) . 

(B.17) 



To decompose the connected Green's functions into one-particle irreducible ones we 
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use the relations 



{Al{v)c^{yy{u)) = J d'[z,z^z,] Dliiv - z,) D'^{y - z,) Ti^^{z,, z,, z^) D, 

(B.18) 

{^(y)c''{z)cf{uy{v)) = J d^[u,U2U3U^u,ue] {dII{ui - u^) D{^{u - u^) 

X D^(^{V - U2)Tl^gl\u3,U5,U^,U2) D^q{u3 - z)^ , (B.19) 



which have been derived in appendix IA.3I 

Substituting these expressions into eq. ()B.17|) we arrive at the final expression for the 
ghost Dyson-Schwinger equation in coordinate space: 

- j d^[uv]Tf^!^'^^''{x,u,v,y) D{f{v - u) 

- i y d'^[zUVUiU2U3UiUr^] T^^gl^^^y, Z, V, u) Df^{ui - U^) D{j{u - Wg) 



X 



Tt''\u,, X, u,) Dfj{v - U2) rf (Ml, Us, U2) D^^{U3 



i j d'^[zuvuiU2U3U4]Tf^lf'^\y,z,v,u) D(^{u-Ui) 

X D^q{v - U2) r^^r(^3, X, M4, U2) (^3 - Z) 
d^[zUVZ,Z2Zs] r(°)''^^(2/, U, V) D%{V - Z,) rf '^(Zi, ^3, X)DI^{U - Z3) 

(B.20) 

where an additional minus signs arises from the interchange of the colour indices / and 
g in the bare four-ghost vertices and from the interchange of j and i in the ghost-gluon 
vertex. 

After performing a Fourier transformation we obtain the respective expression in mo- 
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mentum space 



G 

2 



+ (^) lj£ys J D,^ip - q,) Daiql) 

X r^(p, qi) DG{q2) r^(-P + qi + ^2, 92) ^g(p - 9i - (I2) 

- (-^c ) 2 y f^^big2]r2), i^Glgi) ^g(p - gi - 92) r43h(p, gi, ^2) ^0(52) 

+ (-^^)(|^ j d'qTf{p,q)D,,{p-q)TM,P)DG{q) 

(B.21) 

where the colour traces have been carried out and the reduced vertices defined in appendix 
lA. 31 have been used. 
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Methods to solve DSEs in flat 
Euclidean space-time 

C.l Angular integrals 

Working in Euclidean space-time the four dimensional integrals in the loops of Dyson- 
Schwinger equations can be transformed to hyperspherical coordinates as 

J d'^q ...^^J dq^q^ J d<^> J d^ sm{^) J ddsin^{d) .... (C.l) 



We choose the external momentum to point in the four-direction, enclosing the angle 
9 with the loop momentum q^. All integrands are then functions of the squared external 
momentum, the squared loop momentum and the angle 9 only. The other two angular 
integrals can then be performed trivially yielding a factor of 47r. We thus have the relation 

OO TT 

j d'qfip\ q\ 9)^2nJ dqV J d9sm\9) f(p', q\ 9), (C.2) 



which is used frequently throughout this thesis. 

Sometimes the remaining angular integral can be performed as well using the integra- 
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C.2. Tensor integrals 



tion formulae 



d9 



TT 



x{x-y) y{y-x)_ 



(C.3) 



TT 

/ 



de 



TT 



Q(a^ - y) ^ 0(1/ - x) 



X 



y 



(C.4) 



d9sm^{9) 



TT 

2 



(C.5) 





7r 



TT 



d9sm\9)z = -{x + y) 



d9sm^{9)z^ = ^{{x + yY + xy) 



(C.6) 



(C.7) 



where we have used the abbreviations x := p"^, y := q"^ and the squared momentum z is 
defined a.s z = {p — qY = x + y — 2^/xy cos{9). A derivation of eqs. ()C.3|) . ()C.4|) can be 
found e.g. in |141j . 



C.2 Tensor integrals 

The exphcit expression for the scalar bubble integral /, defined in eq. ()C.8|) . can be easily 
evaluated in Euclidean space-time using the Feynman-parameterisation (see e.g. ref. jH]). 
With the squared momenta x = p^, y = q^ and z = {p — qY the result is given by 



/(a,6,p) := / 

J y^z" 

2 .,,_, r(2-a)r(2-fe)r(a + &-2) 

T{a)T{b)T{4-a-b) ' ^ ^ ^ 



The corresponding tensor integrals can be reduced to scalar integrals by extracting combi- 
nations of momenta and the symmetric tensor 5^,^ according to the symmetry properties 
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of the integrand: 

Jf,{a,b,p) :-- 



4„ ?M 
ya^h 



Ji{a,b,p) Pf,. 



^^9%^ = Ki{a,b,p) Pf,p^ + K2{a,b,p) x5^u, 



■ya^b 
ya^b 



(C.IO) 
(C.ll) 



Li{a,b,p) p^PuPp 



ya^b 



+ L2{a,b,p) X {Pp5^p+p^5pp+pp5py) , (C.12) 
^i^Ml^Ml. = Mi{a,b,p) p^PuPpPa 

+ M2(a, b, p) X {5f,y PpPa + 6pp p^Pa + V PpPm+ 

^lypPfiPa + SuaPpP/i + 5p„p^Py) 

+ M3(a, 6, p) x^ {5^u 5pa + 5^p 5„a + V "^p^) • (CIS) 



The scalar integrals in these expressions are calculated by contracting them with appro- 
priate tensors, writing all scalar products in terms of squared momenta x, y and z and 
applying eq. ()C.9|) . One arrives at 
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T{b)T{8-a 


-b) 




2r(5 — a] 


r(3 -6)r(a 


+ b- 


2) 


r(a) 


r(6)r(8-a 


-b) 




21(4 -a) 


r(4-6)r(a 


+ b- 


2) 


r(a) 


r(6)r(8-a 


-b) 
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„2-a- 



2(-3 + a + 6) 



2-a- 



X 
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X 
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(C.14) 
(C.15) 
(C.16) 
(C.17) 
(C.18) 
(C.19) 
(C.20) 
(C.21) 
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C.3. Analytic expressions for some diagrams in bare vertex approximation 



32 A^^,'-' 



Figure C.l: Momentum routing for the sunset and for the dressing diagram in the ghost 
equation. 

C.3 Analytic expressions for some diagrams in bare 
vertex approximation 

In this appendix we give explicitly the expressions for some diagrams needed for our 
investigation in section |^21 AH algebraic manipulations have been done using the program 
FORM jnni • Our ansdtze for the small momentum behaviour of the ghost dressing function 
G, the transversal gluon dressing function Z and the longitudinal gluon dressing function 
L are the power laws 

G{x) = Bxf^, Z{x) = Ax", L{x) = Cx\ (C.22) 

where we have used the abbreviation x = p^. 

We first evaluate the sunset diagram in the ghost equation given diagrammatically in 
Fig. IC.ll With the bare four-ghost vertex given in eq. ()A.28|) and the abbreviations for 
the squared momenta x = p^, yi = (gi)^, 1/2 = ('?2)^, = {p^ QiY and Z2 = {p — qi — ^2)^ 
the sunset diagram reads 

2(27r)« V2V 2J J J xyi J 1/222 

The factor 1/x in the first integral stems from the left hand side of the ghost equation. 
We now integrate the inner loop with the help of formula (jdQj) and obtain 

^ - 512 vrs 1 2 V ~ 2) V P(l - p) r(2 + 2/3) 7 ^ ^''^ ' ^^'^^^ 

where Zi is the total squared momentum flowing through the integrated loop. The second 
integration is done in the same way. We arrive at 

^ ^3, g'Z.B' faf^ a^ 2 r3(i + /j) r(-3/3 - 1) 



5127r4 V2 V 2/ / r3(l -/3)r(3 + 3/3) 

^3(}^jjryun_ ^(^25) 
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As each integration step eats up the two squared momenta in the denominators of the 
integral kernels only powers of x to the anomalous dimensions of the dressing functions 
in the loop (here 3/5 from three ghost propagators) survive. This mechanism works in the 
same way for all diagrams and explains the pattern in the eqs. ()4.4p . ()4.5|) and ()4.6|) of 
section 14.21 

Next we evaluate the two contributions in the gluon equation needed for the argument 
below eq. ()4.9|) . The explicit expressions for the kernels of two- loop gluon diagrams are 
rather lengthy but the calculation is done along the same lines as in the ghost sunset 
diagram above. Therefore we just give the final results: 

^ -27 g'N^Z.A^ r(-l - 4a) r(l/2 - a) r(3a) r^(l + a) 
^^^^ ^ 4096 7r4 r(4-3(T)r2(2-a)r(3/2-a)r(4 + 4(T) 

2-4<^ (-1 + 3a) (10 + a - 66a^ + 63a^) (5 + 43a + A7a^) 
:= x^'^ (\/')t4"^ , (C.26) 

1536 7r4 (1 + 35)13(2-5)1(3 + 35) 
:= x^'ilVYLTL- (C.27) 

Finally we calculate that part in the dressing diagram of the ghost equation which 
contains the longitudinal part of the gluon propagator for the special case a = 0, 2, where 
L{x) = 1. These are the linear covariant gauges. With the momentum assignments x = p^, 
y = and z = k"^ = {p — qY the longitudinal part of the diagram is given by 

jjdress _ Zi r 4 ^y'^ 98^ 

Ul - -^^Jd<l<l,^^Pu^ (C.28) 

where again the factor 1/x stems from the left hand side of the equation. Writing the 
kernel in terms of squared momenta we obtain 

Ut- = /AA5/ fA-^-^ + ^V (C.29) 

^ (27r)4 J ^ ^ \2z^ Ayz^ Axz^ Axy J ^ ' 

After integration we obtain 



= _^y£^Ai?x^ f----- + oV(0) 



(C.30) 



0. (C.31) 



Although r(0) is formally divergent this contribution vanishes due to vanishing coefficients. 
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C.4 Numerical methods for flat Euclidean space-time 



Here we detail the numerical method we employed to solve the coupled system of ghost 
and gluon Dyson-Schwinger equations, ()3.24|) and ()3.25|) . For the convenience of the reader 
we display the equations again: 

^ - 4\l^^^G{y)Z{z) , (C.32) 



G{x) " J (27r)4 xy 

- Z-i + g — — — G{y)G{z) 



Z{x) 3 i (27r)4 xy 



^N, f A Q{x,y,z) G{yY'~''^^-^''W{zY'-'/^-^'^ 



+^ 3 J (27r)4 xy Z{yY Z{zf ' ^^'^^^ 

Recall X = , y = and z = {q — pY = x + y — 2^J\xy) cos 6. The integral over the 
loop momentum q is transformed to four- dimensional hyperspherical coordinates. Two of 
the four integrals are then trivial and yield a factor of 47r. The remaining angular integral 
and the radial integral have to be performed with the help of numerical routines. We use 
the Gauss-Legendre quadrature rule described in |171j . To achieve high accuracy we split 
the radial loop integral into three parts, y G [0,e^], y G (e^,x] and y G {x^Xjjv]- The 
second split is necessary as the integrands are not smooth at the boundary x and too 
much accuracy would be lost, if one uses a quadrature rule that spans the whole region 
(e^a;c/y]. 

According to the value of their argument the dressing functions Z as well as G have to 
be handled differently. In the infrared region, y, z & [0, e^], Z and G behave like powers and 
are replaced according to eq. ()3.44j) . Recall that the approximation by leading powers in 
the infrared is justified by the analysis of subleading contributions at the end of subsection 
13.3.21 The infrared matching point is chosen sufficiently low. In physical units we 
have = (0.55 MeV)^ in our calculations. In the high momentum regime, y G {x,xuv]j 
arguments z occur which are larger than the numerical cutoff xuv- There we approximate 
the respective dressing functions by the expressions 

Ziz) = Z(/)[u;log(y)+l]' , (C.34) 

G{z) = G(/) [c^log(y) + (C.35) 

according to the one loop behaviour of the solutions as has been detailed in subsection 
13.3.11 Here u = llNca{l)/12TT and the squared momentum / is a perturbative scale. We 
chose / = (174 GeV)^ to be slightly lower than the numerical cutoff xuv = (177 GeV)^. To 
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Figure C.2: The y-dependence of the function z) for different values of the parameter 

a. Only the choice a = 35 leads to momentum independence in the infrared. Due to the 
symmetry of the ansatz for Ziiy, z) the 2;-dependence is the same for b = a. 

be able to perform the angular integrations for momenta [e^jXt/y] we expand the dressing 
functions in Chebychev polynomials and solve the coupled system of equations for the 
expansion coefficients using a Newton iteration method. Details of this technique can be 
found in appendix B of ref. [FTTj . 



C.5 One-loop scaling 

In the framework of the truncation scheme presented in section 13.31 we have shown that 
the substitution 

Zi^Ziix,y,z;s,L)= ^^^y,^,^ (^-36) 

for the gluon vertex renormalisation constant Zi together with a bare three-gluon vertex 
yields the correct one-loop scaling of the gluon loop in the gluon Dyson- Schwinger equation. 
This is true for any values a and b. Of course, in a full treatment of the coupled ghost- 
gluon system Zi{s, L) would be independent of momentum. Therefore a choice of a and b 
which keeps Zi as weakly varying as possible seems the most reasonable one. 

This choice can be inferred using the scaling of the dressing functions extracted from 
the renormalisation group equation, see ref. jHEl for details. The dressing functions can be 
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expressed as 

1+2(5 

-5 

R-\x) , (C.37) 

where the running couphng provides the correct one loop scahng in the ultraviolet. Con- 
sequently the function R{x) approaches unity for high momenta. Furthermore, from the 
known infrared behaviour of Z{x), G{x) and a{x) (c.f. subsection I3.3.2|) one infers that 
R{x) is proportional to x'^ in the infrared. Writing Zi in terms of a{x) and R{x) yields 

. {^^f''„,, (c,38) 

In the perturbative region -R(?/), R{z) 1 and the function Zi is therefore slowly varying 
for any a and h according to the logarithmic behaviour of the running coupling a. In the 
infrared, however, a approaches its fixed point while the functions R behaves like a power. 
Consequently the choice a = h = ?>5 guarantees the weakest momentum dependence of Zi 
which is illustrated in Fig. IC.2I Shown is the y-dependence of the function Zi{x^y^z = 
y;s,L). (Note also that due to the symmetry Zi{x,y, z; s, L) = Zi{x, z,y; s, L) and the 
absence of an explicit x-dependence this is sufficient to demonstrate its momentum de- 
pendence.) In the perturbative momentum regime the function Zi does not vary with 
the parameter a. So all three choices give the same logarithmic running in momentum as 
required to give the correct one loop scaling behaviour of the integral. In the infrared, 
however, a change in a gives rise to substantial changes in the behaviour of Zi, with only 
the choice a = 36 leading to a constant. 



Z(x) 



G{x) 



a{x) 

a{x) 
a(s) 



Appendix D 

Discretisation and finite volume 
effects 

D.l Radial discretised DSEs 

In this appendix we study possible discretisation errors in the Dyson- Schwinger equations^. 
We investigate a discretised version of the angular approximated DSEs in the ghost-loop 
only truncation scheme discussed in subsection 15.1. 21 From the infrared analysis of these 
DSEs one finds power laws for the ghost and gluon dressing functions. These power laws 
are used in the numerical treatment of eqs. ()5.13|) . ()5.14|) to solve the integrals from zero 
momentum to an infrared matching point e^. It has been claimed that these integrals 
are crucial to find numerical solutions for the ghost and gluon propagators j^. In this 
section we will show that they are not. We will introduce a prescription to discretise the 
radial integrals in the DSEs such that the smallest momentum showing up in the equations 
is nonzero. Despite the infrared divergence of the ghost dressing function we reproduce 
the solutions from the formulation with continuous momenta without having to include 
the region in the very infrared. This is very important as otherwise the attempt to solve 
Dyson- Schwinger equations on a four-torus, performed in chapter El might be hopeless 
from the very beginning. 

For the discretisation of eqs. ()5.1Hj) . ()5.14|1 we use the same prescription for the radial 
integral as one gets for all four Cartesian momentum directions if employing a finite volume. 
Note that contrary to the finite volume case there is no easy geometrical interpretation 
for a discretised radial integral. 



I owe Jacques Bloch the idea to the study performed in this section. 
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With / denoting the inverse spacing of the radial momentum g, we thus substitute 

2 



dq 2q 



3 



2ti 
-J 



2vr \ \ ^ 



:d.i) 



into eqs. (jOT^ . ^X^. We obtain 



1 



1 

Zixi 



9 g'^Nc /27r 



4 487r2 V / 



i-l 



487r2 



i=i 

i-l 



;d.2) 



N 



2j 



2%- 



(D.3) 



The squared momenta Xi are given as Xj = (?27r//)^. The largest momentum xat 
corresponds to the numerical cutoff L = in the continuous version, eqs. ()5.13|) . ()5.14|) . 
The momentum xq = is discarded in the calculation. We solve the discretised eqs. ()D.2j) . 
()D.3j) by iteration employing the Newton method. To directly compare the continuous 
solution with the discretised one we use the same cut-off and renormalise at the same 
momentum. This is most easily done by taking Z3(/i^,A^) and Z3(/i^,A^) in eqs. ()D.2|1 . 
()D.3j) from the continuum solution. 

It is quite amusing that there is a second way to solve eqs. ()D.2|) . ()D.3|) . which cor- 
responds to the conversion of the continuous integral equations ()5.13|) and ()5.14|) into 
differential equations^, see ref. ^J\: Subtract eqs. ()D.2j) . ()D.3|1 from themselves for two 
different arbitrary momenta Xi and Xg- One obtains 



1 



9 g'^N^ (2tx 



G{x, 



4 48772 V / 



min{i—l,s—l) 



j=mm(t,s) 



(D.4) 



Z{xi) 



Zixs) 487r2 



2j f y] , 3y, 



Xe 



2x, 



mm(j— l,s — 1) 



G{y,) 



j=min{i,s) 



G\y,] 



(D.5) 



^Note that such a conversion is only possible for angular approximated DSEs. 
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Figure D.l: Comparison of the gluon dressing function Z{x), the ghost dressing function 
G{x) and the running coupling a{x) obtained from the DSEs with continuous momenta and 
from the discretised version. For the radial grid of the discretised version we used three different 
momentum spacings. 



If we now set s = 1 the sums in the equations run from j = Given the values for 
Z{xi) and G{xi) the equations for Z{x2) and G{x2) are then simple quadratic equations 
which can be solved easily. Given these solutions we can solve the equations for i = 3 
and so on. The complete solution is thereby built point for point from the infrared to the 
ultraviolet region of momentum. The advantage of this method is that one does not have 
to use a numerical iteration procedure to solve the equations. Of course, the specification 
of the input values Z{xi) and G{xi) corresponds to a certain choice of the renormalisation 
constants Z3(/x^, A^) and Z3(/x^, A^), i.e. to the boundary restriction from a certain choice 
of renormalisation point. The disadvantage of this method is thus that one has to improve 
iteratively on the input values Z{xi) and G{xi) until the same Z3(/i^,A^) and Z3(/i^,A^) 
are obtained as for the continuous integral equations. Both ways of solving eqs. ()D.2jl . 
()D.3jl certainly yield the same results. 
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Our solutions for the dressing functions on the radial momentum grid are compared 
with the solutions from the continuum in Fig. ID. II For the dressing functions and the 
running coupling we find nearly the same results with both methods. There are only small 
discretisation errors for our largest momentum spacing. For smaller spacings these errors 
become more and more irrelevant and we thus conclude that there is a smooth transition 
to the continuum as the spacing more and more decreases. 



D.2 The influence of zero modes on the solutions of 
DSEs on a torus 

An important point when formulating the Dyson-Schwinger equations on the torus could 
be the treatment of the zero modes. In addition, on the torus an infrared analysis like 
the one in flat Euclidean space-time is not possible, and one is left with the problem 
how the dressing functions behave at vanishing momenta. Guided by the intuition that 
especially the long ranged modes should be affected by the finite volume we assume in 
the following Z{x — ^ 0) = just like in the continuum and G{x — > 0) = const if zero 
modes are neglected. Phrased otherwise we assume that the zero modes are the missing 
ingredient to ensure the correct infinite volume limit for the torus results. Therefore, if 
on tori of different volumes G{x = 0) shows no sign of becoming divergent, the infrared 
enhancement seen in G{x —>■ 0) or in the flat space-time results has to be due to the torus 
zero modes of gluons and ghosts. 

Therefore, in this appendix, we will show that the assumption G{x = 0) < oo does not 
lead to a contradiction in the equations on the torus if zero modes are neglected. To this 
end we focus on the truncation scheme without angular approximations. First we rewrite 
eqs. dSm, (HT^ as 



r A Q{x,y,z)G{y)-'-<^'Giz)-'-<^' 
^ 3 J (27r)4 xy Z{yf^ Z{zf^ ' ^ ' ' 



According to Appendix 1(151 we have chosen a = b = 36, where 6 = 



—9/44, the anomalous 
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dimension of the ghost. The kernels have the form: 

If {x-yf\ lfx + y\ 1 sin^e , 



M(.,!,,.) = ifiz^. + I-il^Vi + l^-l-, (D.9) 



z \ A 2 A X J 2 2x Ax 

I flx"^ 2 19-C 5-C 2 Cl/^ 

Q{x,y,z) = + x — ^y + —^y+Q — 

z^ \8 y 8 4 8 X 



+ -[ — X — y + C— 

z \y 4 4 X 

19- Cx 3-4C 9C?/ 
8 y 2 ^ 

C , 5-C^ , .2 C 



~'~ \a; ~'~ 4y / ~'~ 8xy ^ ^ 

We first analyse the behaviour of the integrands in the limit y ^ for finite momenta 
X. Then Z{z) — > Z{x) and ^(2;) — » (^(a;) and the kernels times the respective dressing 
functions are to appropriate order in momentum y: 

9MmK(.,y,z) G(0)Z(.f^, (D.ll) 

xy x^ 

^M^M{x,y,z) ^ G'(0)G'(x)^(l + (C-2)cos2e) , (D.12) 
xy 

Gjyy'-'^'Giyr'^'^' G{oy^-^'G{x)-^-^^ / Ccos^e \ 

Z(y)3'5Z(^)35xi/ ^^-^'^'^^ Z{0)^^Z{xy'xy V ^•••y' • ^'^•'''^ 

Furthermore, z = x + y — 2^fxy cosQ has been used and terms proportional to cosB 
have been dropped, as they either integrate to zero in the continuum or cancel each other 
in the sums on the torus. Each of the expressions ()D.11|) . ()D.12|) . ()D.13|) is then the 
appropriate term for j = on the right hand side of the Dyson-Schwinger equations on 
the torus. Clearly one observes that only a finite ghost mode G(0) avoids trouble with 
divergences. This is especially true for the kernel Q' of the gluon loop, as Z~^^{y — > 0) 
is more singular than the simple pole, so this kernel vanishes for small momenta y. The 
other two expressions ()D.11|) and ()D.12|) are finite. One is then left with the ambiguous 
quantities sin^ B and cos^ B which will be replaced by their integrals from zero to 27r in the 
calculation at the end of this section. The arbitrariness of this procedure is considerably 
moderated by the observation that any number plugged in for the trigonometric functions 
yields the same qualitative result at the end of this section. 

Second, we take the limit z — > 0, which on the torus is identical to B ^ 0. The ghost 
kernel sin^B/2;^ alone would certainly diverge as B ^ 0, but taking into account the 
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Figure D.2: G(0) compared with the results for finite momentum x on the torus. For 
convenience we have kept the logarithmic momentum scale and plotted the zero modes on the 
left border of the figure. 



power law behaviour Z{z) ~ z'^'^ for the gluon dressing function the integrand is zero in 
this hmit. This is vahd for k > 0.5, which is in agreement with the infrared analysis in the 
continuum. We therefore may omit the points z = in the ghost equation. The situation 
is different in the gluon equation where the kernel of the ghost loop has a finite limit z —>■ 0: 
M{x,x,0)/xy = (C + l)/(2x^). Therefore with a finite ghost dressing function G{0) the 
points z = in the ghost loop contribute but no divergences occur. In the gluon loop the 
kernel Q' multiplied by the dressing functions approaches zero for vanishing momentum z 
due to the power law behaviour of Z~^^{z 0). 

To obtain a definite value for G{0) we now investigate the behaviour of the eqs. ()D.6|) . 
()D.7j) in the limit x — ^ 0. The integrands are then given by 

^M^A-(x.,..) ^ G(y)Z(yf^, (D,14) 

xy y^ 

9mAM(.,y,.) - G(,)G(,)(i^i^2!!^ + .„) . (D,15) 

xy \ xy J 

where the kernel Q' is of no interest, as we know the gluon loop to be subleading in the 
infrared. Clearly, the kernel M of the ghost loop in the gluon equation is now singular for 
X — i> 0, corresponding to a vanishing gluon dressing function in the infrared. This result 
confirms our working hypothesis that the gluon mode Z{Q) is not affected by the finite 
volume of the torus. The integrand of the ghost equation is finite up to the point y = 0. 
There the pole in the kernel is cancelled by the behaviour of the gluon dressing function 
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Z[y) ~ y'^'^ resulting in a zero for vanishing momentum x and y. 

We therefore arrive at a consistent set of equations for the ghost propagator at = 
and a vanishing gluon Z{0). In Fig. ID.2l we show the resuhs for the ghost dressing function 
gained on two different volumes on the torus. Within numerical accuracy the values of 
G{0) are the same for the two volumes. Obviously terms with high loop momentum y 
contribute most to the right hand side of the ghost equation for vanishing momentum 
X. Furthermore, one observes that the actual value of G{0) is not in accordance with an 
extrapolation of the ghost curves to the infrared. There is also no sizeable change of the 
gluon and the ghost dressing function when G{0) is set to zero by hand. This has been 
done in all calculations performed in chapter El 

Having shown that G{0) < cxo on the torus even in the infinite-volume limit and 
assuming that the torus should provide a reasonable infrared regularisation of physics in 
fiat space-time we conclude that the divergence of G{0) is very probably due to the torus 
zero modes of gluons and ghosts. Noting furthermore that a diverging G{0) is related to 
Zwanziger's horizon condition and the Kugo-Ojima confinement criterion this indicates a 
direct relation between zero modes, the Gribov horizon and confinement. 
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